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ABSTRACT 


The impingement of turbulent jets is studied both experiment- 
ally and analytically. This work involves mainly impingement upon 
solid, smooth boundaries, however, a simple case of a jet impinging 
against a uniform stream is also considered. Three distinct flow re- 
gions exist in this phenomenon; the free jet region, the impingement 
region and the wall jet region. Emphasis has been placed in studying 
the impingement region which is the most interesting for practical 
applications and the least studied owing to its complexity. The rela- 
tively simple case of axisymmetric impingement is studied first. De- 
tailed explorations of the time-average velocity and pressure fields 
are performed and semi~empirical methods are proposed to predict the 
velocity and pressure fields and the wall shear stress in the impinge- 
ment region. Experimental data of previous works in the wall jet 
region are re-analyzed and shown to comply with a simple similarity 
analysis. The results are summarized and average values of pertinent 
coefficients are obtained. The more complex cases of plane and circu- 
lar jets impinging obliquely are studied next and semi-empirical methods 
are proposed to predict the wall pressure in the impingement region. 
The wall pheag stress in this region is predicted by extending the 
previous method to account for the effects of obliqueness. Relatively 
simple methods to predict the velocity field and wall shear stress in 


the wall jet region are developed and are shown to agree with previous 


and present experimental results. These studies refer to cases where 
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the height of impingement is sufficiently large to allow the assumption 
that the jet is a point source of momentum. In order to understand the 
physics of the problem when the impingement height is "smali", axisym- 
metric impingement at heights comparable to the jet diameter is also 
studied in detail. It is shown that significant differences exist in 
the impingement region whereas the wall jet region remains unaffected 
of "small" or "large" height considerations. Further, the development 
of a circular jet in a uniform, opposing stream is studied as an im- 
pingement problem. A semi-empirical method is developed to predict the 
flow field, based on turbulent jet theory and potential flow considera- 
tions. Finally, as an application of knowledge gained thus far, a 
criterion for scour from jets is developed and is shown to agree with 


previous experimental results. 
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CHAPTER I - INTRODUCTION 


The impingement of a turbulent jet on a solid surface is of 
interest in many engineering problems, such as paint spraying, jet blast 
drying, shielded arc welding and in connection with the operation of 
VTOL aircraft. In the field of hydraulic engineering the impingement of 
jets is of particular interest in the design of hydrauiic outlet works 
and in other problems related to scour from jets. The phenomenon of 
scour from jets has attracted the attention of hydraulicians (Refs. 1 - 
8) in the last 35 years, however, no satisfactory solution has been 
obtained so far. The reason for this, appears to be a general lack of 
appreciation for the mechanics of impingement. It was felt, therefore, 
that the first step in a logical development of this subject ought to be 
the systematic study of the impingement problem. The present study is 
mainly an effort toward systematic evaluation and solution of the prob- 
lem of turbulent jets impinging on smooth, plane, walls. This problem 
has been attacked in the past (Refs. 9 - 27) by workers of various 
engineering backgrounds. The great majority of these works deals with 
plane or circular jets impinging normally, while very little is known 
about the oblique impingement which is naturally more interesting from 
the practical point of view. For purposes of introduction, some very 
general considerations will be outlined in this chapter whereas a com- 


prehensive review of literature will be presented in the next chapter. 


Consider a submerged, turbulent jet impinging normally ona 


smooth wall, as shown in Figure 1.1. The outline of the jet shown by a 


solid line, is to be interpreted as the locus of points where velocity 
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becomes negligibly small. At present, it has been established that there 
exist three distinct flow regions (I, II, III in Figure 1.1). In region 
I, the flow develops as if the wall did not exist, because wall effects 
which are transmitted upwards in the form of increased static pressure 
become negligibly small beyond a certain distance from the wall. This 
implies that in region I, the flow is identical to the well-known free 
jet (28),(29) and hence this region is appropriately called the free-jet 
region. In region II, the effects of the wall are significant. The 
static pressure assumes values higher than the ambient pressure and in- 
tense pressure gradients are set-up in both directions, parallel and 
normal to the wall. These pressure gradients, in turn, cause the flow 
to deflect and become eventually parallel to the wall. Region II is, 
therefore, called the "deflection" or, more frequently, the "impingement" 
region. After turning, the flow becomes almost parallel to the wall. 

At the same time, the static pressure drops to ambient values after some 
distance from the stagnation point and a flow pattern is established 


which resembles in many respects the classical wall jet. Region III is, 


therefore, appropriately called the wall-jet region. 


A complete study of this problem requires investigation of all 
three regions. - However, the flow characteristics in region I are well 


known already (28),(29) and all that is needed to know is how far from 


the outlet this region extends. Region III is a wall-jet type of flow. 


Considerable theoretical and experimental information is already avail- 


able for the classical wall jet (Refs. 30 - 35. Reference No. 35 is an 


extensive bibliography on the subject and contains many references which 


are not listed here). The objective, therefore, in reg 
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establish the similarities with the classical wall jet and to account 
for differences that are likely to exist owing to the different geo- 
metries that produce the wall jet. Coming to region II, one could see 
that this is the most complex region of the impinging jet. Here, both 
components of velocity are of the same order of magnitude and pressure 
gradients are significant in both directions, hence boundary layer sim- 
plifications do not apply. Furthermore, since the direction of the 
velocity vector and the static pressure at any arbitrary point are not 
known beforehand, it is obvious that equipment more sophisticated than 
the ordinary Pitot-tube is necessary. This calls for rather laborious 
measurements and probably accounts for a general lack of detailed ex- 
perimental data in this region. At the same time, the impingement 
region is the most important from the practical point of view, because 
here the hydrodynamic action on the boundary is most severe, owing to 
simultaneous existence of large shearing stresses and significant 
pressure gradients. Hence it is only natural that in the present work 
emphasis was placed on the study of region II. Region I was studied 
only in so far as was necessary to verify its free-jet character and to 
establish its boundary with region II. Intensity of study in region III 
was dictated mainly by the amount of information already available. It 
was found that ny few experimental data would suffice to complement 
those already available, hence the main task in region III was to improve 


upon the analytical solutions proposed earlier, and to extend the theory 


to more advanced problems, such as oblique impingement. 


Given the complexity of the problem, it was considered impera- 


; Sees ‘ 
tive to proceed in an inductive manner. The simplest case of impingement 
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is the normal impingement of a plane jet. This was studied by the author 
as part of his M.Sc. program (25) and the results were published recently 
(26). The corresponding axisymmetric case, i.e., a circular jet imping- 
ing normally, was studied next and the results are presented in Chapter 
IV. Experimental and analytical results of these works were obtained 
for impingement heights larger than about ten nozzle diameters (or nozzle 
widths) which means that the jet was fully developed when approaching 

the wall (28),(29). It was thought, however, that these results would 
not apply to smaller heights where the jet is not fully developed, due 

to the changes in the physics of the problem. In order to explore 
impingement at small heights, a study of axisymmetric impingement was 
carried out keeping the impingement height less than 6 nozzle diameters. 
The results, which verified the expectation of significant differences, 
are presented in Chapter V. The findings of Chapters IV and V consti- 
tuded the foundation for additional work on more complex situations. 
Owing to the relatively minor practical importance of impingement at 


small heights, the remaining studies were mainly for large impingement 


heights. 


The case of plane jets impinging obliquely was studied next 


(Chapter VI). Here the number of parameters defining the problem is 
increased by one, namely the angle of impingement. A further step in 
complexity is the case of circular jets impinging obliquely. Here, in 
regions II and III, the flow ceases to be axisymmetric and, therefore, 


in addition to the angle of impingement, a new parameter must be intro- 


duced to account for variations in the direction normal to the plane of 
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A practical application of the knowledge gained thus far, is 
given in Appendix A, where a criterion for scour from circular jets 
impinging normally, is developed. Finally, another problem of imping- 
ment is considered (Chapter VIII), namely the development of a circular 
turbulent jet in a uniform opposing stream. This study is of interest 
with respect to efficient disposal of effluent discharges. A critical 
review of literature and some fundamental theoretical considerations 
are presented in the next chapter, whereas descriptions of experimental 


set-ups and experimental techniques are given in Chapter III. 
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CHAPTER II - REVIEW OF LITERATURE 


In this chapter, an effort is made to review critically the 
existing material on turbulent impinging jets. The literature cited 
herein contains accounts of previous contributions to the author's best 


knowledge, however, no claim to absolute completeness will be made. 


2.1 Stagnation Flow Problem 


Consider first a very simple case of impingement, i.e., the 
case of an infinitely wide, uniform stream, impinging normally on a 
smooth, plane, wall. This is the well-known stagnation flow problem 
(36) which affords an exact solution of the Navier-Stokes equations of 
motion. For simplicity, consider a two-dimensional (plane) flow, the 
axisymmetric case being analogous to this. With reference to Figure 2.1, 
the velacity distribution in the frictionless, potential flow outside 


the boundary layer of thickness 6, is given by (36): 
Us ax $V Uaeeenay @s1) 


where a is aconstant. If Ps is the stagnation pressure and p is 


the static pressure at an arbitrary point (x, y), then: 


1 2. xe, yr 
x Ss + (252) 
Pp Pp 7) a (x y ) 


Ss 


For the laminar flow within the boundary layer, it is assumed that: 


Ele 6a (0) Sah i ola SO) (2.3) 


p.-p = xpalix’ + Fy) (2.4) 
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FIGURE 2.1 PLANE STAGNATION FLOW 
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FIGURE 2.2 CLASSICAL WALL JET 


In this way, the equation of continuity is satisfied identically and the 
functions f and F are determined by substitution into the Navier- 
Stokes equations and subsequent reduction to ordinary differential 
equations. It is noteworthy that the assumed functional form of the 
u~component (Equation 2.3) implies that the boundary layer thickness, 

6, is independent of x. If the co-ordinate y is normalized by the 
transformation Hn = 2 y, then it is shown that the normalized 
velocity component, u/U, is a function only of n, i.e., the velocity 
profiles are similar. If 6 is defined to be the value of y for which 


u/U equals 0.99, then 6 is given by: 


2.40 (plane flow) 
yi 5 - (2.5) 


1.98 (axisymmetric flow) 
The wall shear stress, Ty» can also be calculated and is given by: 


1.2326 (p a Vav) x (plane flow) 
2:0 


1.312 (9 a Vav) r (axisymmetric flow) 


where r is-radial distance from the stagnation point. It is seen that 
te increases linearly with distance from the stagnation point. Natural- 
ly the above considerations are meaningful so long as the boundary layer 
remains laminar. Stagnation flow has been found to occur not only near 
plane walls, a also near any cylindrical body (36). Considering a 
turbulent mee ke jet, it may be noted that in a small neighborhood of 


the stagnation point, the approaching stream will have an almost uniform 


downward velocity (time-averaged). If at the same time the intensity of 


turbulence is sufficiently reduced, then stagnation flow is likely to 


occur within this small region. 
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2.2> ‘The Classical Wall Jet 


The problem of the wall jet produced by a two-dimensional or 
radial nozzle and growing on a smooth wall has received considerable 
attention (Refs. 30 - 35) and can be considered solved. The now classic 
work of Glauert (30) appears to be the first theoretical treatment of 
this problem, including solutions for plane and radial wall jets and 
for both laminar and turbulent states. In the case of turbulent wall 
jets, the theory was based on matching separate solutions for the inner 
and outer layers (see Figure 2.2) using the eddy viscosity model of 
turbulence. Later experimental investigations have shown that such 
assumptions are not realistic and, therefore, it was thought instructive 
to base the following presentation on the work of Schwarz and Cosart 


(31) which is most appealing for its logical and deductive development. 


Consider a plane jet issuing from a two-dimensional nozzle of 
thickness d, with a uniform initial velocity Uo: growing on a smooth 


wall as shown in Figure 2.2. In this case the static pressure is 


constant; the transverse component of velocity, v, is much smaller 


than the forward component, u, and the transverse extent of the flow 


is relatively small. Hence the boundary-layer approximation can be 


applied to the Reynolds' equations of motion, which after simplifications 


become: 
OT Z 
du ou 1 i du (2.7) 
uatbVv a FOS Qh ae . 
Ox oy OQ oy he 
and Fi (2.8) 
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where u, v are the time-average components of velocity, is is the 
turbulent shear stress and v, 9 are the kinematic viscosity and 


density of the fluid respectively. Experiments have confirmed that, 


much like in the case of a free jet, the flow is self-preserving, i.e.: 
Beers | Sy) S, (2.9) 


where 6, is the larger value of y at which u = u_/2. Assuming that 


the turbulent shear stress obeys a similar relation, if scaled with the 


maximum local dynamic pressure, then: 
2 * 
05 = pu g(n) (210) 


Gc m 


Using Equations 2.9 and 2.10 and 2.8, Equation 2.7 can be reduced to: 
n if] 
papas f° ae Ldn |= 2 ce Peas 
ie che dn dx dn 
m 
fe) (@) 


$B 4 v £ (2511) 


In order that Equation 2.11 be satisfied for all n, the coefficients 
of functions of nm must be independent of x. Therefore: 
dé. 
_ dx 


a ee ae 


* Strictly speaking, similarity of shear stress can only be assumed in 


(212) 


the form tT. = k(x)g(n), where k(x) is a function to be determined. 
c 


It is, however, possible to show [see for example Poreh, Tsuei and 
Cermak (17)] that k(x) can be taken equal to the wall shear stress, 
t . Further, by assuming that the defect law holds in the inner layer 
nee 2.2) and comparing to Equation 2.9, one can show that T/pu_, 


2 
is independent of x, i.e., k(x) can be replaced by pus 
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6 du 
L ae apni f A 
< dx = const = Se (2543) 
m 
ar oe Bes c 
eer coee en ey (2.14) 


A further relation is obtained by integrating Equation 2.7 from y =o 


to © , which yields: 


2 d Z 


uy = hicuees (S,u, ) C205) 
7 oe 
where F =f fad @@y constisitand aw, wiisethe shear weilldcity VT /p. 
e) 
Relations 2.12 and 2.13 imply that: 
‘peihax pegs s and a = C ee (2.16) 


2, 2. m u 


where Cc. and §$ are constants and x is a corrected value of x, to 


account for virtual origin effects. Using the latter and letting Ce 
ou 
to be the local skin friction factor Tf A » Equation 2.15 becomes: 


Sil ls . (clan ai (ais) 
Z 
which also shows that coe is independent of x. Strictly 
speaking, the exponent f is always less than - $. However, the term 
om is generally a small quantity, so that as a first approximation: 


1 
2 == 2.18 
B ; (2.18) 
It is obvious now that with this value of 8, relation 2.14 cannot be 


satisfied simultaneously with 2.12 and 2.13. Yet for values eft 


greater than about 0.005, the viscous term in Equation 2.7 can be neglec- 


ted and hence the condition of Equation 2.14 will vanish. For n < 0.005, 
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viscous effects will be significant with a resulting loss of self- 
preservation. The above quoted limiting value of n was proposed by 
Schwarz and Cosart (31) on the basis of data in the inner, or boundary, 
layer. Since 5/6, ~ 0.150, it may be stated that self-preservation is 
lost only in a layer of thickness equal to about 3% of the boundary 
layer thickness. It is noteworthy that calculations and measurements 
of the turbulent shear stress (31),(17) have shown that it does not 
vanish at the locus of maximum velocity (n = 6/6,), a fact that pre- 
cludes use of any eddy diffusivity or mixing length model. This is 


believed to be due to the asymmetric character of the profile. 


For the case of a radial wall jet, analogous considerations 
are valid, however, here un, will vary approximately as et The 
focal” shin friction factor’, Ces can easily be shown to depend upon 
the nozzle Reynolds number, Ro = Ud/y, by means of dimensicnal con- 
siderations. However, when more complex geometrical configurations are 
atest ect ty such as jets impinging normally or at an angle, dimensional 
analysis will not be very fruitful due to the increased number of para- 
meters. Insight as to dependence of Ce upon initial parameters of 


the problem can be gained from the following reasoning. 


It is assumed that the velocity u obeys a defect-type law 


in the boundary layer, i.e.: 


u-u 
iy es na ROR e 
uy £65) sale) a= 


This law is normally valid for turbulent boundary layers on smooth and 
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rough walls alike. Schwarz and Cosart (31) have indirectly verified its 
validity in the inner layer of the wall jet on smooth walls, whereas 
Rajaratnam (34) has found it to hold in the case of rough walls. Apply- 
ing Equation 2.9 at y= 5, it can be shown that 6 = const. 6, : 


therefore, the above equation can be re-written as: 
7. 6 
= £,(n) 3 my 5. 
* 2 
Solving this for u/u, and comparing to Equation 2.9, it is easy to 


show that: 


Me - AMe-i 


“in £,() 


= const. 


which implies that the friction factor, Ce = 2() ,» depends only upon 
0 

the shape of the functions f and £7 These apes pel are not likely 
to change with the geometry producing the wall jet, but they can change 
slightly with Ro the nozzle Reynolds number. Hence it is postulated 
that Ce depends only upon Ro? being independent of the geometry which 
produces che wall jet. 
2.3 Turbulent Impinging Jets 

Consider a turbulent jet, plane or circular, impinging on a 
smooth wall, as shown in Figure 2.3, where for simplicity the impinge- 
ment is ried sic be normal. In this figure, d is the width (or 
diameter) of the nozzle, Us is the uniform velocity at the nozzle, 
and H is the impingement height. For oblique impingement, H is 
defined to be the distance between nozzle and wall, measured along the 
centerline of the jet. The co-ordinate system x, y (or r) is used 
for regions I and II, whereas the system x (or r), z is used for region 


III. The difference hetween the static pressure at any 
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FIGURE 2.3 TURBULENT IMPINGING JET 
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and the ambient pressure is denoted simply by p, and will be referred 


to as "pressure". The wall pressure is’ p_, Py is the stagnation 


w 
pressure and To is the wall shear stress. A short paper by Murray 
et al, Rouse (9) in 1956 contains the results of a study to determine 
the flow pattern and some gross characteristics of a circular jet 
impinging normally. (Values of the parameters defining the problen, 
such as the angle of impingement 4, the dimensionless height H/d 


and the nozzle Reynolds number Ro = Uid/y, are listed, for this, and 


for the following references in Table 2.1, at the end of this chapter.) 


The impingement height H was made so large compared to the 
jet diameter that the jet could be considered as a point source of 
momentum. A fundamental contribution of this work (9) is the proposition 
that when H/d is very large, the two parameters ie and d do not 
influence flow properties separately. Rather, they can be combined 
into a single parameter, known as the kinematic momentum and equal to 
ndU,/4 ae reasoning borrowed from free-jet theory). This immediately 
leads to the realization that there must be a limiting value of H/d, 
below which the physics of the problem will change. Two questions are, 
therefore, raised: (i) how does the flow behave when H/d is "small" 
(i.e., smaller than the limiting value), and (ii) what is the limiting 
value of H/d. Very generally, it may be stated that the answer to 
(ii) will be related to the length of the potential core of the jet, 
whereas the answer to (i) will be related to the distinct physics of 
an impinging jet and an impinging shear layer for "large" and "small" 
values of H/d respectively. (In the remaining part of this thesis, 


the term "impinging jet" without further qualification will denote 
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large impingement heights.) 


It was reported in this paper (9) that the impingement region 
begins at a distance Xo from the nozzle equal to about 0.8 H, and the 
well-known irrotational flow solution for this region was presented. 
Measurements of wall pressure yielded the surprising result that the 
integral of the wall pressure was about 30% more than the value of the 
initial momentum of the jet. This was attributed to the occurrence of 


negative pressures on the wall. 


In 1959, Bradshaw and Love (10) presented the results of an 
experimental study for a circular jet impinging normally. These include 
data on velocity, static pressure and shear stress in the impingement 
region. For example, it was shown that the wall shear increases from 
zero at the stagnation point to a maximum value occurring within the 
impingement region, decreasing rapidly afterwards. Similarity of radial 
velocity profiles in the wall-jet region was verified. The length scale 
6, increased linearly with r whereas the velocity a (Figure 2.3) 


varied approximately in inverse proportion to r. Self preservation in 


the wall jet begun at r/H®0.45. 


In the same year, Poreh and Cermak (11) presented the results 
of an experimental and analytical study on the same problem. The well- 
known flow development in the free-jet region was verified. Some data 
are presented in the wall-jet region, however, the choice of velocity 
scale was made in a manner that would force the profile to assume the 
shape corresponding to the free jet. Hence, it was necessary to esti- 


mate the velocity scale by extrapolation, a procedure that is inaccurate 
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and similarity of the profiles was not established conclusively. In the 
impingement region, the wall pressure and radial velocity were measured 
and the results were presented in a generalized form by means of dimen- 
sional analysis which is essentially identical to that of (9). For 
example, pH /M, (OM = initial momentum of jet = ond” /4) is plotted 
against r/H and a single curve results, showing that pH /M, = Ge 
The irrotational flow solution was compared against data on wall pressure 
and was found to be valid within a radial distance of r/H < 0.05 from 


the stagnation point. 


A technical report by Schauer and Eustis (14) in 1963 contains 
the results of an experimental and analytical study of plane turbulent 
impinging jets. This includes measurements of wall pressure and one 
measurement of wall shear in the impingement region. Some data were 
reported for angles of impingement other than 90° (oblique impingement). 
The analysis was based on an integral approach which led to rather com- 
plex equations that can be solved only implicitly. Furthermore, a large 
number of rather strong assumptions were introduced, a fact that probably 
accounts for a large portion of the analysis failing to describe the ex- 
perimental results. It is worthwhile to examine these assumptions. For 
the impingement region it was assumed that: 

(i) The engin extends from the stagnation point equally along the 
wall and along the centerline of the jet. There is no obvious 
reason why this should be so and none is given. A weak justi- 
fication can be furnished only in the case of normal impingement 
if one refers to the corresponding flow of a potential jet, which 


again is related to the potential flow of two equal jets imping- 
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ing against each other. In this case the central streamline be- 
tween the jets, after impingement, is straight and can be replaced 
by a wall. Because of the reversibility of potential flows, the 
pressure distribution on the centerline of the original jets is 
identical with that along the wall. If the impingement region is 
defined so that the pressure is greater than, say, 1% of the 
stagnation value, then the impingement region will extend equally 
in both directions. 

The stagnation pressure is assumed equal to the centerline dynamic 
pressure at the boundary between regions I and II. For normal 
impingement, where the centerline is also a streamline, all the 
way to the stagnation point, this assumption implies inviscid 

flow at least near the centerline. For oblique impingement, where 
the stagnating streamline does not coincide with the centerline, 
the flow needs to be, in addition, irrotational. These constitute 
strong simplifications and they do not correspond to reality. For 
example, for normal impingement it was found (26) that the dynamic 
pressure at the boundary of regions I and II is about 30% larger 
than the stagnation pressure. 

A similarity function was used to correlate wall pressures, which 
is srimierntcs about the stagnation point. The data show, however, 
that for @¢ # 90°, the wall pressure distribution is skew and 
deviation from the assumed curve becomes significant for ¢ < 70°. 
For © # 90°, the stagnation point is removed from the intersection 
of jet centerline and wall by a certain distance, s. In order to 


predict this quantity, an equation was borrowed from the corres- 
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ponding potential flow solution which, however, gives the distance 
to the center of gravity of wall pressure. This is consistent 
with the previous assumption of symmetry, but in view of the actual 
skewness of the wall pressure, it does not correspond to reality. 
(v) In order to obtain an extra integral equation, entrainment in the 
impingement region was neglected. Careful examination of later 
data (26) shows that significant entrainment continues to take 


place well into the impingement region. 


In the wall-jet region the analysis was based on the following 


assumptions: 
(i) The local skin friction factor Ce = 2(uy/u_)” was assumed to 
u_6 
vary with the local Reynolds number R = a according to 


a Blasius type relationship. Since in the wall jet R is approx- 


imately proportional to mn it is in fact assumed that 
Ce oc x te This contradicts previous considerations that 
showed Ce ec ty but in view of the small power of xy 


involved, the discrepancy is not important. Hence, this assump- 

tion could be viewed as a source of unnecessary complication. 
(ii) It was further assumed that the shear stress T is given by 

Prandtl."s hypothesis, i.e., Tr © uié6,——. This, however, is 


contradicted by reality as was mentioned before (Sec. 2.2). 


In 1964, Tani and Komatsu (15) presented the results of an 
experimental and analytical study of a circular jet impinging normally. 
Based on experiments with values of H/d equal to 4, 8 and 12, they 


reported that the impingement region extends about 1.6d to 3d radially 
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from the stagnation point and 1.6d to 2.2d above the wall. No distinc- 
tion was made between "large" and "small" impingement heights. In this 
work, this constitutes an ommission, because the value H = 4d is 
smaller than the potential core length (about 6d) and, therefore, one 
could expect. it to be classified as a "small" impingement height. On the 
other hand, the value H = 12d is about twice the length of the poten- 
tial core and, therefore, it would rather belong to the "large" impinge- 
ment height problem. An interesting theoretical contribution is included 
in this paper and it deserves closer examination. It was assumed that 
the flow in the impingement region is inviscid but not irrotational. In 
this case the equations of motion will not contain viscous or turbulence 
terms, which makes it possible to derive a differential equation for the 
streamfunction (Stokes) W that does not contain other unknowns. The 
streamfunction was then written as an infinite series expansion (using 


the co-ordinate system r, z as shown in Figure 2.3): 
Yo = $x? [£() + reg (2) t roh(z) + 0, 1 (2.19) 


where f, g, h, .... are as yet unknown functions of z. Substituting 
Equation 2.19 into the differential equation for it is possible to 
obtain e¢, h, .... in terms of only the function f(z). The z-component 


of velocity, u, is then given by: 


2 4 
wbtaleat pnpictea ty ee = | 2? [ae (ett 24) eno) 
4 2 2 8 So oe 
dz dz 
where A, bo» . are integration constants. The meaning of f(z) is 
made obvious by letting r= 0. Then, Mes Gar Ele) tl Cut eis) MepGes 


sents the magnitude of the centerline velocity. The function f(z) has 
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to be evaluated empirically for any given set of initial conditions and 
this constitutes a weakness of this approach. Furthermore, for large 
impingement heights the assumption of inviscid flow is justified only in 
a small neighbourhood of the stagnation point rather than throughout the 
impingement region, as was pointed out earlier. In Chapter V, however, 
it is shown that this theory can be employed to advantage for the cases 
of small impingement heights, where the total pressure head remains 
constant on the centerline. The results of an experimental study of 
plane jets impinging normally, were presented by Cola (16) in 1965. 

One measurement of vertical velocity in the impingement region is inclu- 
ded and the results agree with the experimental and theoretical results 


of Beltaos and Rajaratnam (26). 


A paper by Poreh, Tsuei and Cermak (17) in 1967 is a theore- 
tical and experimental study of the wall jet region in the case of 
axisymmetric impingement. The theoretical development is similar in 
many respects to that of Schwarz and Cosart (31) which was outlined in 


ection 2.2. Firstly, it was shown that: 
ge coer ; On. “SF : Ce weeaer (25219 


It was not realized that Cy is independent of r (Section 2.2), which 
would immediately give n=l. The values of the exponents m and n 


were obtained empirically as: 
m = -1.1 and n = 0.9 (2,22) 


However, a re-analysis of their data (Chapter IV) showed that 5, varies 


linearly with r and if a small virtual origin correction is introduced, 
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then b, oc tr; Ls@eg oD From the integral momentum equation, the 
following relation can be derived (which is analogous to Equation 2.17 


for the plane case): 


& Ce 
m= - (1 Any (2323) 


where F is defined as in Section 2.2, and Cy = dé,/dr. It follows 
that, strictly speaking, m< - 1. However, from the available data the 
term C,/4EC, can be estimated to be very small in comparison to unity. 
From dimensional considerations, and using n= 1, Ce can be shown to 
depend only upon the nozzle Reynolds number, Ro» so that it decreases 
with increasing Ro: This means that for very small values of Ro» 
negligence of the above term may not be realistic. Assuming then that 
Ro is not very small, the value of m can be taken as -1l. Again, a 
re-evaluation of the data presented in (17) showed that Un, x“ j/r 
(Chapter IV). The value of this approximation is not merely academic. 


From dimensional considerations, the authors of (17) have shown that: 
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If m=-l and ne=1 the above expressions will become free of the 


parameter H. This constitutes a major simplification, especially when 
considerations of oblique impingement come into play. Then the com- 
plexity of the problem will be considerably reduced by knowing that H 


is irrelevant. 
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Measurements of turbulent stresses confirmed the conventional 
boundary layer approximations with respect to normal stresses, and at 
the same time demonstrated that root-mean square values are also self- 
preserving with scale ut However, due to scatter it was not possible 
to establish conclusively the similarity of turbulent shear stress 
variations across the jet, even though a trend toward similarity is 
clearly demonstrated. It is interesting to note that the turbulent 


shear stress did not vanish at the point of maximum velocity. 


In 1968, Cartwright and Russel (18) presented an experimental 
and analytical study of a plane jet impinging normally. Measurements 
were reported mainly for the wall jet region. The law of the wall was 
verified for the boundary layer. The linearity of 6, was also obtained 


but the exponent of x, describing the variation of u_ was given as 


a 
-0.39, a rather surprising value, since as has already been pointed out 
(Section 2.2), this exponent cannot be greater than -0.5. Yet, when 
the data were re-analyzed (Chapter VI) it was shown that an exponent 
equal to -0.5 is reasonable, if a small virtual origin correction is 


allowed. In the impingement region only the wall pressure was measured, 


and a similarity property was obtained by plotting P/ De vs x, /H. 


An interesting numerical study of the plane turbulent jet 
impinging ee was presented by Wolfshtein (19) in 1970. It was 
assumed that the Reynolds stresses are related to the mean strain by a 
scalar eddy viscosity, which in turn was assumed to depend on the level 
of energy fluctuations and a length scale. Thus, it was possible to 


obtain iterative, finite-difference solutions which were compared with 
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the data of Schauer and Eustis (14). No claim to general validity of 
the theoretical model was made owing to the arbitrary selection of the 
length scale. However, this model predicts the correct trends with the 
exception of the wall shear in the impingement region. This was pre- 


Big betes whereas later work (26) showed 


dicted to be proportional to 
that the wall shear stress in region II is practically independent of 
Ro° It is noteworthy that the model predicts strong generation of 
turbulence near the stagnation point, a phenomenon that was also obser- 
ved experimentally (Ref. 10 of Wolfshtein (19)). Owing to the very 
small thickness of the wall boundary layer in the impingement region, 

it is extremely difficult to explore this boundary layer experimentally. 
Yet, knowledge of boundary layer properties is important, not only in 
an academic sense, but also from the practical point of view for both 
scour and heat transfer considerations. In order to estimate boundary 
layer thickness from data on wall pressure and shear stress, knowledge 
of the state (laminar or turbulent) of the boundary layer is imperative. 
The above considerations point toward a turbulent boundary layer, a 
possibility that is supported by the irrelevance of Ro with respect 
to wall shear stress. The turbulent counterpart of the classical 


stagnation flow problem (Section 2.1) has been treated by the same 


author, Wolfshtein (37). 


A systematic experimental study of a circular jet impinging 
normally was published in 1970 by Hrycak et al (22). The variation of 
potential core length with R, was studied extensively. It was demon- 
strated that for Ro < 4,000 this quantity normalized with the nozzle 


diameter, d, depends stron 
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tial core length varies only between 6.8d and 6.5d. Next, the free jet 
laws were verified in the free-jet region. In the impingement region 
measurements of wall pressure and maximum radial velocity ve (maximum 
in a vertical, occurring at the edge of the boundary layer) were pre- 
sented. No distinction was made between small and large impingement 
heights, hence no generalization of the results could be achieved. 
Valuable data on boundary layer thickness were also presented. These 
were for H/d = 3, 7, 10 and 20 but for one only value of Reynolds 
number (Ro = 54,000) and were estimated to be accurate to within 15% 

and 30% for H/d = 3 and 20 respectively. It was established that near 
the stagnation point the boundary layer thickness 6 is almost constant. 
Using the available data on vO to evaluate the constant a of Equa- 
tion 2.5, the thickness can be predicted from stagnation flow consider- 
ations (Section 2.1). However, the predicted values are about one-half 
of the measured ones (owing to a numerical oversight, Hrycak et al (22) 
predicted values that are only 30% less than the actual). In the wall- 
jet region, u, was found to vary as ¢ by'ty yet a re-analysis 

showed that the relation un « ]1/r gives also a good fit (Chapters IV 
and V). It is noteworthy that these data were for values of H/d rang- 
ing between 4 and 20, i.e., for both small and large impingement heights. 
Yet, when these data were plotted in the form u/U, vs r/d_ they 
collapsed in a single curve, thus leading to the expectation that the 
wall jet region is independent of small or large height considerations. 
This in turn would imply that the wall-jet region depends only upon the 
initial momentum of the jet for all impingement heights. Further 


support to this is furnished by the data of Yakovlevskii and Krashen- 
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ninikov (27), not only for normal impingement but also for oblique as 
will be seen in Chapter VII. The length scale, Oo» was given as 
6,/d = const Gade. the constant being a weak function of H/d. 
Again, re-analysis showed the linearity of So» Sey 6,/d = C, (r/d) 
with Cy taking values between 0.071 and 0.082 as H/d changes from 
2 to 20. This, however, does not necessarily imply dependence upon 
H/d, because C, has been found to vary even for classical wall jets. 


2 


For example, C, varied from 0.056 to 0.085 in the work of Schwarz 
and Cosart (31) for a plane wall jet. It is believed that this varia- 


tion is due to varying nozzle conditions which cannot be accounted for 


at present. A survey of literature by the same authors, Gauntner et al 


(21) was also published in 1970 and presumably preceeded the above work. 


A brief paper by Bradbury (23) in 1972 is an attempt to cor- 
relate wall pressure data by means of a rather weak dimensional argument. 
This correlation is essentially identical to that of Poreh and Cermak 
in 1959 (11) and can be derived formally by means of dimensional ana- 
lysis using the assumption of a point source of momentum as set forth 
by Murray et al, Rouse (9) in 1956. This paper seems to ignore most of 


the significant research prior to 1972. 


A paper by Kamoi and Tanaka (24) in 1972, is mainly an investi- 
gation of turbulence properties near the stagnation point for various 
angles of plane jet impingement. It is again demonstrated that strong 
turbulence exists near the stagnation point. Measurements of r.m.s. 


(root-mean-square) values of wall pressure and shear stress showed that 
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for large impingement heights the r.m.s. is maximum at the stagnation 
point itself, whereas for small impingement heights the maximum is re- 
moved somewhat from the stagnation point. In fact, the authors state 
that coincidence of r.m.s. maximum with the stagnation point indicates 
that the boundary layer is turbulent. The state of the boundary layer 
depends upon the parameter H/d, the impingement angle and the nozzle 
turbulence level. For normal impingement the critical values of H/d 
are 8 and 4 for turbulence levels of 0.5% and 3.7% respectively. Time- 
average values of wall pressure and shear stress were also measured for 
some runs, however, the shear stress data appear to be suspect because 


of non-zero values of Es at the stagnation point. 


In 1973, Beltaos and Rajaratnam (26) presented the results of 
a systematic experimental and analytical study of plane jets impinging 
normally. The impingement region was studied in detail, including ex- 
plorations of the velocity and pressure fields. It was shown that the 
impingement region begins at x/H = 0.7. Flow properties were predic-— 
ted semi-empirically and the analysis will be illustrated in Chapter IV 
for the axisymmetric case. The wall-jet region was analyzed according 
to considerations outlined in Section 2.2. It was established that 
for small values of Ro the exponent of xy describing the variation 
of us cannot be taken as -0.5, being definitely smaller than this 


value. 


Finally, a paper by Yakovlevskii and Krashenninikov (27) deals 
with the wall-jet region for oblique impingement of circular turbulent 


jets. This will be discussed in greater detail in Chapter VII. 


cot DN Pe ey 
Pe yen 
ae: . ¥ 
S , Vy ; 
* { 
X it 
' ras 4 
10 Pe “s 
' fe Set Rel iy 


pobkyeagure sig #4 pen ter by ¥ Os + ot 
meek tiene. a eadghed snteight COR, 4 i ae) 


pies syataun ans (apet a .anied oo 


royce) vxebenod ae i¢ asase ety es td ive it ie 
on oid his ohanr ivonitnye Fah Jetta Oni sacha aa ie 


(h Io eauhaybepds ip ett -aeiomens bak Cibiees! sakes ; 
=uctd F fay py oSgeon We UE Site SEO 6 aiaited ois tytn) «re 


2 - been sem “ORRe AnH fab ti Sa Liew ¢iwaes sy hit Yo" 


— 


ainkan psaqens ‘ef en aia ad in rH daereIe ated ” sows 
ae sittin suk aokanga afr hall i ut ‘ne sntean onan 
baits ele 

Yo-qsfubey oy noraieeg (as) ‘ence ‘Shs ueoaiad ERRE * 


a a 


aqehigrebeat gat welq “pu ae barca 7. ua tng 3 
“tn, on bivient fies a es yshiaw wel HOT 29%, sosmognceat our ! oe: 


id a 


ait ints wae: tala # sootett ‘pahacig brie. vation pee Je 


-Stieae gyay pa Ranmifinty ‘eutt- t6 = Whe 92 anrgrd olgas +8 


/t wetaeHo aa Buidpanantlt a8 Apis disease “Ai ais rt nel 
Te 


netheooon, pein, env ngtaos ee-Hiey oar one. 2 a eaemenya tae Ls ie: 
a 1 
7 ; 


nce Ai lider. no baw di° ~ SR nepense mt ‘Seats navoid surat 


abil tua oA Gniddtrseah th Fo Foose qns is . ‘9 eoute I 
a : a , 


Life: 


aba nett) qothma Claslatioh Sica ae O- i quan 3d toate ‘2 a 


ba os a ‘ . ps 


vale paeY 
acid 5 t 
EA 


30 


2.4 Summary 


A critical review of literature pertinent to the problem of 
turbulent jet impingement, has been presented in the previous sections. 
It was shown that previous work was mostly concerned with regions I and 
fat. in region I, flow sWavearenierics have been conclusively verified 
to be identical with those of a free jet. In region III, an improvement 
upon the existing analytical approaches has made it possible to effect 
considerable simplification so long as the Reynolds number, Ro» is 
larger than a certain limiting value which will be discussed in later 
chapters. This simplification consists of showing the wall jet develop- 
ment to be independent of the height of impingement. Hence, the velocity 


and length scales are given by: 


for plane jets 
6, = (2.24) 


Cy $) for circume jets 


» for plane jets 


E 


(2-25) 


(=) 


> tor circular jets 


where Cos C. are constants (different for plane and circular jets), 

and the overbar denotes that distance along the wall must be measured 

from a virtual origin. Regarding the wall shear stress, it was shown 

that it is proportional to the local value of the maximum dynamic head 
i.e., the local skin friction factor is independent of xy (Or Lys 


This factor depends only upon R, i.e., 
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(G = Ce (RD (2526) 
So that |-dC2/aR' << 0. 
£ fe) 


Furthermore, it was pointed out that if H is made very small, 
the physics of the problem is likely to change, thus bringing out the 
distinction between "large" and "small" impingement heights. The above 
considerations are ordinarily valid for "large" heights, however, based 
on existing data, it was tentatively shown that the wall jet region de- 
velops in a unique way, independently of large or small height consi- 


derations. 


Finally, turbulence measurements (24) show that the state of 
the boundary layer in the impingement region depends upon the relative 
height H/d, the angle of impingement, $, and the nozzle turbulence 
level. Even for the smaller turbulence level investigated (0.5%) the 


boundary layer will be turbulent for all 6 so long as H/d > 8. 
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CHAPTER III - EXPERIMENTAL EQUIPMENT AND TECHNIQUES 


Sr Experimental Set-Up 


Experiments were conducted with air jets in air using speeds 
less than 300 fps so that compressibility effects were negligible. Use 
of air instead of water has the advantage that the necessary experi- 


mental set-up is less bulky and much easier to handle. 


Circular and plane jets were produced by means of suitably 
designed nozzles and plenum chambers. The experimental set-up for plane 
jets is shown in Figure 3.1. The aspect ratio of the nozzle is about 
65, so that the flow could be assumed to be two-dimensional in the cen- 
tral portion of the nozzle in the region studied. In the chamber, a 
series of circular tubes were placed near the entrance and a sequence of 
five screens was placed further downstream, as shown in Figure 3.1(a). 
The tubes serve the purpose of "straightening" flow streamlines which 
tend to become curvilinear after the sudden expansion at the entrance. 
The screens were provided in order to ensure mixing and thus uniformity 
of the flow. The flow rate was adjusted by means of a pressure regula- 


tory valve. 


To produce impingement, a 1/4 inch thick plexiglass plate, 
6 inches wide and 5.75 feet long, equipped with 1 foot high sidewalls, 
was used. Adjustment of the angle of impingement was made possible by 
means of a simple suspension system which is outlined in Figure 3.2, 
whereas Figure 3.3 shows photographs of the plane jet set-up. A circu- 


lar plenum chamber 3.5 inches in diameter was used for circular jets 
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(a) Set-Up 
Air supply 


Plenum chamber 


! 
4 plexiglass 
, sidewalls 


Static outlet 


Air supply 


ot 


Tubes 


(b) Nozzle and Chamber Details 
5 screens 


16 thick galvonized iron 


H 0.088" 


FIGURE 3.1 SET-UP FOR PLANE JETS 
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FIGURE 3.2 CHAMBER SUSPENSION-PLANE JET 
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PLANE JET SET-UP 
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and the details of construction are much the same as those of plane jets. 
This was made out of plexiglass and was possible to equip with nozzles 
of different outlet diameters. In the present experiments two outlet 
sizes were used, 0.923 and 0.253 inches. Impingement was produced by 
means of plexiglass plates. The set-up was somewhat different for nor- 


mal and oblique impingement, as shown in Figure 3.4(a) and (b). 


The temperature of the supply air was somewhat less than the 
room temperature (= 68°F) and for this purpose a heater was installed 
between supply and chamber, so that the air issuing from the nozzle was 


kept close to room temperatures. 


Finally, for studying the circular counterjet, a smaller plenum 
chamber was built (plexiglass) and was placed along the centerline of a 
2 foot square wind tunnel. This chamber was constructed in the same 
fashion as the previous two but its diameter was only 1.5 inches so as 
to minimize disturbance to the approaching stream. A nozzle with outlet 
diameter of 0.20 inches was used so that the ratio of stream dimension 
to nozzle diameter was made very large (120). To minimize boundary layer 
effects, the chamber was placed near the entrance of the wind tunnel. 


This set-up is shown in Figure 3.5. 


3.2 Experimental Techniques 


3.2.1 Nozzle Velocity 


Referring to Figure 3.6, the Bernoulli theorem applied between 


sections 0- 0 and 1-1, gives: 
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(a) NORMAL IMPINGEMENT 


(b) OBLIQUE IMPINGEMENT 


FIGURE 3.4 CIRCULAR JET SET-UP 
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FIGURE 3.5 


SET-UP FOR COUNTERJET EXPERIMENTS 
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FIGURE 3.6 APPLICATION OF BERNOULLI THEOREM 


(b) TOP VIEW 


FIGURE 3.7 THE PITOT CYLINDER 
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where p is the mass density of the fluid, Pp, U, and Eno, are 
values of static pressure and velocity at sections 1-1 and 0-0 
respectively. sig 8 A sAy are the flow areas at sections 0-0 and 
1-1, then (u,/u,)° = (a,/A,)°. This quantity is always less than 


0.0067 for the present experimental equipment, and therefore, the pres- 


sure difference, Ap = Pree is given by: 


oo (3.2) 


The validity of Equation 3.2 was checked by comparing manometer heights 
with nozzle velocity heads at different velocities bec No difference 
was recorded between corresponding readings. A vertical water manometer 
was used for heads larger than 3 inches. Assuming that reading errors 
were at most +0.05 inches, the maximum head recording error would be 
about: 1.62. Since os « YAp , the maximum error in obtaining uy would 
be less than 1%. Thus, the nozzle velocity was measured indirectly by 
recording the pressure Ap. This procedure has the advantage that ue 
may be recorded within seconds at any instant, so that corrections for 


small supply fluctuations can be applied. 


3.2.2 Velocity and Pressure Measurements in the 
Impingement Region - Pitot Cylinder 


In the impingement region the direction of the velocity vector 
is not known before hand. In addition, the static pressure is higher 
than the ambient (atmospheric). As ordinary Prandtl-Pitot tubes cannot 
handle this situation, a Pitot cylinder was used, employing a technique 
that does not require rotation of the instrument. This technique was 


developed by the author (25),(38) and is summarized below. The Pitot 


ots es 4% tates 9. q or or 2 
if : 7 

0-0 be t+ Ft "aksobe 9 wrtasioe tan pega ots 

bape 0 - 9 enoksoe. 26 enbes molt i sm, 

ned evel eynwee’ pt into waae ya 


-soxy ada _,seitemed? hie, sma 


ad hg tend ereeveers andres et soci ine $2 oc to t 
snomsiigh o 0 geRaootay rearsad tay te sbsod meet 
nonam masa Leokeaew A, Te cmconiresiiiyl oct 
ett, galls send gamma “agndaad € cat aigin stad Yok a oH 

Oe bass wc ‘a, «abcou ea.0t eon a ia 
nf op 9 eed kb os79 awed oom ails - ‘a wo wo uae 
| efsqw wed Semana cies egnpsitn ‘ebaton sont “SE ned» vt Le 


Ft) 


dels 1 agama: oad, on grubeoayy: eat as eresnore oda. oie 


” 


py? all shY¥a orig’ ja, as eine ws ss ssanee shite be poksebas od 
Wakes tie se ite saapilies ‘ela, 


Tatoo eibonlay anys) Bee pels aa iota ets at 
‘ai Sra eat ah 169% ait noe ok basil avohed cial 30H at 


YDase Pel at Buus “! 
rongeo ets Yoo LM “-cang euand be eA - Cobpmiqvoms) , dnatden sit madd 


supindomde & h, gat yoga. Late, ar sebntitto sore & vaplsauiin atid etas 


gee suplatoss ekat - ‘mesa yd i aphsador erlupos Jou vache: set 
JOSLT gle “nfad Sas ‘ofan a Spares. (oR, (2s) ern od, yd, beqaloweb 


d 


“~~ 


41. 


cylinder is a cylindrical instrument of small diameter having three 
orifices near its edge which is rounded (Figure 3.7). For two-dimensional 
flow, the instrument is placed with its axis perpendicular to the plane 
of flow so that no gradients exist along the axis. If the orifices l, 
2, 3 (Figure 3.7) are connected to a manometer, three pressure readings, 
Py> Po» Px» will be obtained Peeom Lei AG ch are related to the local 
values of 9, V, and Py (static pressure). On this principle the 
instrument can be calibrated so that measurement of Py> Po» and Py 
enables calculation of 96, Vand Po? and the procedure is: 

1. Measure Py> Po» Ps 

2. Compute K = (PD, = P,)/(P, 7 Po)- 

Ss. From a calibration plot of K Vsid,. obtain >. 

4, From a second calibration plot of K, vs o, 

G01, 2) 3). obtain Ky. 


5. Compute velocity from: 


= an ae = j = ue 2 (3753) 


DeSoto : 5 hoe A Ae pe (3.4) 


Calibration curves for the factors kK, Ki> Kos Ky and ss - K, are 
shown in Figures 3.8, 3.9 and 3.10 for an instrument 1/8 inch in diameter 
with a polar angle 8 between orifices of 45°. This method was found to 
be satisfactory when applied to plane shear flows (25). It was also 


found to give satisfactory results in the case of axisymmetric shear 
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FIGURE 3.10 CALIBRATION FACTORS K-K, AND K3-K, (1/8 in. Probe) 
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flows, because here, even though the flow is not two-dimensional, local 
gradients along the axis of the probe vanish if the probe is oriented 
with its axis perpendicular to the radial line passing through orifice 

2. For the case of large impingement heights, where the width of the 
jet is much larger than the probe diameter, the 1/8 in. probe yielded 
satisfactory results. However, when the same probe was used in cases of 
small impingement heights, erroneous results were obtained, which were 
attributed to the relatively large size of the probe. In order to over- 
come this difficulty, a new probe of diameter 0.04 in. was built. Owing 
to the small size of the probe it was not possible to construct separate 
orifices. Hence, only one orifice was drilled and each measurement 
involved rotation by 45° twice so that readings corresponding to 1, 2, 3 
could be obtained. This, however, did not slow down the work, because 
the inside diameter of this probe was considerably larger than the inside 
diameter of the fine tubing leading to the orifices of the 1/8 in. probe, 
so that response times were much smaller in this case. As the orifice 


size can affect the calibration curves, a new calibration was performed, 


the particulars of which are shown in Figure 3.1l. For the factor 
K, in the Figure 3.11, the corresponding curve of the 1/8 in. probe 
is also shown for comparison. It is seen that the curve is shifted 


somewhat due to differing relative size of the orifice. For measure- 
ments in the wind-tunnel, associated with the counterjet problem 
(Chapter VIII), a third probe of diameter 0.04 in. and having only 
one orifice, was built. In this case, however, the relative size of 
the orifice was made to be close to that of the large probe. Cali- 


bration curves for the wind tunnel probe are shown in Figure 3.12. 
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FIGURE 3.11 CALIBRATION FACTORS Ky > Kos K, (0.04 in. .Probe) 
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FIGURE 3.12 CALIBRATION FACTORS Ky> Kos K, (Wind Tunnel Probe) 
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It is seen that in this case (Figure 3.12) the calibration curve for Ky 
is almost identical to that of the 1/8 in. probe. In order to estimate 
effects of probe size to the measurements, an error analysis was perform- 
ed (Appendix B) using the experimental results for small impingement 
heights. The percent error depended upon the local transverse gradients 
of pressure and velocity. The maximum error was 2% for velocity, 3.5% 
for the static pressure and 1° for the angle 6. The percent error 

in velocity and pressure increases generally for increasing distance 


from the centerline. 


To measure wall pressure, a large number of 1/16 in. diameter 
static taps were installed on the impinging plates. Details of the 
spacing for plane jets are shown in Figure 3.13. For normal impingement 
of circular jets, pressure taps were located along a straight line as 
shown in Figure 3.14(a). For inclined impingement of circular jets a 
circular plate 2 ft. in diameter was used, with taps along a radial line. 
This plate was able to rotate about its center. Radial distributions of 
pressure were measured at 15° increments of polar angle. These were 
plotted separately and the co-ordinates of points where p/P. = 0563955 
.9, .8, etc. were determined. These were afterwards plotted on the wall 
plane and contours of equal pressure were drawn. The taps were spaced as 
shown in Figure 3.14(b) and a close view of this plate is shown in 
Figure 3.15. Wall pressures yielded vertical water columns generally 
less than 3 inches, hence, a precision manometer was used which can 
measure alcohol columns at slopes as small as 1:25. To facilitate the 
work, a valve was placed between taps and manometer so that it was pos- 


Sible to take readings from different taps by simply turning a knob. 
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PLATE AND TAPS-OBLIQUE IMPINGEMENT 


FIGURE 3.16 VALVE AND MANOMETER 
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D24 


Figure 3.16 shows this arrangement. 


In the impingement region, a thin boundary layer grows on the 
wall beginning at the stagnation point. The maximum velocity in the 
vertical, which occurs at the edge of the boundary layer, was measured 
using ordinary total head tubes and assuming that the static pressure 


at this location was equal to the wall pressure in the same vertical. 


3.2.3 Measurement of Wall Shear Stress 

Preston's technique (39) was employed to measure wali shear 
stress with the aid of the calibration curves given by Patel (40). This 
technique is briefly discussed below. For a Preston tube of external 
diameter qd resting on a smooth boundary, if Ap is the difference 
between the total pressure indicated by the tube and the static pressure 


on the boundary, it could be shown that: 


Apd T Wd 
Q bi (22) (3.5) 
2 
; 2 2 2 2 . ‘ 
Setting -Ap..= Apd /4ov"—s and Ae ad Tat. /4ov', Equation 3.5 becomes: 
Be let Daal) ew) , (3.6) 


Using the equations given by Patel (40) for different ranges of Ap, 


and T Equation 3.6 is plotted in Figure 3.17. This curve is normal- 


ox? 
ly valid for zero pressure gradient boundary layers. Certain limitations 
exist, therefore,when the pressure gradient is not zero, as is the case 


in the impingement region. Patel (40) suggests that in favourable pres- 


sure gradients use of the Preston technique will involve errors less than 
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6% if the following conditions are satisfied: 


(i) 0>A>-0.007 or |A] < 0.007 
(ii) u,d /v < 200 (3.7) 
(iii) a eae 


where x is distance along the wall, and A = x 3 “2. ; u, = /t /p 


pu, 
These criteria are rough guides rather than absolute restrictions. Con- 
dition (i) was derived empirically with the aid of fence readings. Con- 
dition (ii) ensures that the tube is located in the region of wall simi- 
larity. Condition (iii) is imposed to entre that the flow is not sub- 
jected to re-laminarization. In the present case, the only physically 
plausible sequence of different flow states is firstly laminar and later 


turbulent. 


Two total head tubes with outside diameters of 0.049 and 0.046 
inches were used as Preston tubes in the present set of experiments. 
These tubes were also used for measuring velocities in the wall jet 
region, as the static pressure becomes equal to the ambient (atmospheric) 


in the wall jet. 
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CHAPTER IV -—- AXISYMMETRIC IMPINGEMENT 
4.1 Introduction 


A relatively simple case of impingement is studied in this 
chapter. Emphasis has been placed in exploring time-average flow charac- 
teristics in the impingement region. Data of previous investigators in 
the wall jet region are re-analyzed according to theoretical considera- 


tions outlined in Chapter II and the results are summarized. 


Lee Experiments and Experimental Results 


The experimental set-up has been described in Chapter III. 
Two nozzles of internal diameters of 0.923 and 0.253 inches were used. 
The velocity at the nozzle was varied between 153 and 293 fps and in 
this range compressibility effects are known to be negligible. The 
nozzle Reynolds number Ro was in the range 10° to 100. The height H 
was varied from about 15 to 20 inches, thereby varying the ratio H/d 
Eremsz liz to5./. On the whole, five experiments were conducted (see 
Table 4.1) and for description of various symbols see the definition 
sketch, Figure 4.1. The variation of the u-component of velocity with 
r at different x-stations, in regions I and II, is shown in Figure 4.2. 
These profiles are replotted in a dimensionless form in Figure 4.3. 
From Figure 4.3 it was found that for = less tian about 0.95, all.pro- 


files are described by one curve. (These profiles are said to be similar.) 


Further, this general curve could be given by the equation: 


=“ =. exp (- 0.693 my) (4.1) 


; Kee 
yo ee Sa) # ! 
a ae 


Fle win: eg 
sit nh SSP Roe at seine te quan. nei, 
_sanratig wilt AGRE aE ge ORaRNS ob bapsta non 
mt axqugars cov eva to nod a the 


+ 


Saeco ea8 stoner i fe eee 
YG ( te | fl' 


tT weagnld jet fodProel it os 
fan, ey. andbith avid. as eo) : 

wk Bins ait set. bivs wel, cattaed bobsge 

ad tq bigen, o. oe otom vad $300 ‘iiaaenenien 
a snd stot | 
b\ taer pl i ee -ausioak 08 ast stings govt _ a 
mee), ‘vayoboes on aanamasene> aveth pated ats Ps 0,802) pee 5 

7 — ay giao: ‘98 sa8), ates ‘aunts do wok ixoumb od bec ED 
| dphy estsaLsy id Aaa day Qe dolswkyey oft” 1-8 aright dosail 


2 


. 


nee bicneitl mt cies ‘i il trent ‘ones mk: vwaotanaenn jnez0ikh a6 ) 


ro) esting Fst wis ae cpeihesotis 5 of beasolqet ot jsitierd ouodt 
“org. DRG aoe ‘i duos mated seek ; 40% eile neue saw rs fut omupht oot 


e +5 ¢dimde oe od ODES ase eth c— saath) satan aero “ bsdtooneh id goit? 


nahi mit ed govig sd Dives ova trian elds rodst 


, : 7 7 : BS g 
(a) | ay oon oS. ca00 (oa 


7; 


305 


SS eee 


SINAN YEdXa AHL dO STIVIYG 1) Pa Ley 


III uofsey | II uopTsey |} [ uorey 


Set {TIueny peansesy 


vi 


% 


Ot = 1e°9 €6¢ 


OT 4: 70-8 991. 


LIC 


96T 


occ 0 


€c6°0O 


CSc 0 


€76°0 


=09 OT 


£5 61 


29S 


G61 


mS 


Sire Oh kag See 


« * : ae 
‘e* 


=. 


DP. 


HOLAAS NOILINIAYG T°? dans 


NOMA YLSIO 
“adg | U7 ASNSS3ad TivVM 


T1IVM soe"d 
d Sq 
aes LLL eg I le LM 
| ° i Prise: NOILVNOVIS 
_ 
NOoIO3y | ra ye 
NOILNGIYLSIG 3YNSS3yd 
ae wae Bed Filvisess4aDxa 
NOIS3¥ 13f 33a4 | 


A NOLAgINLSIC 
ALIDOIZA WWDILYBA 


\ para. 
4 ; n s : Y 
a1ZZON—_S 


ire (q) 


| 
3 


° - 


7 3 Pe =) 
oye a eee 


VOITUsIEe ee 


arama ogni | Lf NTE j yeh 


aXe 


58 


FIGURE 4.2 VELOCITY PROFILES IN REGIONS I AND II 


120 
ts 
H 
100 © 0.463 
40632 
00.752 
a0 ¢ 0800 
4 0.900 
— 60% ® 0.952 
7 0976 
i © 0.986 
eo DATA FROM 
RUN No. 1 
20 ¥ 
< +, *—s—___s 
0) ——O--—-——+» *——___, 
1 3 
: r (in.) oY 
yg Po seep 7 


4 0.900 


: * 0632 ¥ 0925 
le as 20752 #0952 
* 0800 «0976 


© 0.986 


f(y )= 


¢ 0.928 

¢ 0.955 
0974 
RUN No.3 


FIGURE 4.3 DIMENSIONLESS VELOCITY PROFILES IN REGIONS I AND II 


i a eee 5 Ws ar a 
; a 
4 
t 
t 


Sabie BS Swit : 2 


= ie +: 


_—. PF tne Waivers cone 
he SS Sh 80, Ba aaa 30 ar An £o 


zhi i ce j 
; j a 


‘or 


hl Wie PSVOLOIR AT BROOM RS GMey ees MoTeMaMIG” , £., SAVOTS 


_) P. FI 


ag. 


where un is the maximum value of u at any x _ station, occurring at 
the centerline, n = r/b with bi being the value of r where 

u = Fu, and e is the Napierian base. The quantities u. and db. 
are generally known as the velocity and length scales. The data of 
Figure 4.3 are also well described by the Tollmien solution (28),(29). 
For values of x/H greater than 0.95, the velocity profiles are no 


longer similar. 


The variation of the dimensionless scale u/U, with x/d is 
shown in Figure 4.4. For x/H up to about 0.86, Uo/u, increases 
linearly with x/d, (the virtual origin appears to be located at the 


nozzle itself) thus verifying the equation of the free jet: 


ae a 
UO - » ed (4.2) 
oO 
where Cy = 6.8 and 6.3 for runs 1 and 3 respectively. This type of 
variation for C has been noticed before (29). The variation of the 


ih 


dimensionless length scale b /# with x/H is shown in Figure 4.5. 

For x/H up to about 0.86, Bi increases linearly with x, witha 
slope of 0.093, which is close to the generally accepted value of 0.096 
for free circular jets (29). ~ For x/H > 0.36; be increases more 
rapidly, assuming a maximum value at x/H So Die For x/H > 0.95, db. 
decreases, but in this region the velocity profiles are not similar and 


bo is of little practical interest. 


In the impingement region, the v component was also measured 


and Figure 4.6 shows that = (uv) is approximately independent of x 


for small values of r. This property is of importance in predicting 
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u, in region II and has been observed earlier by the author (25), (26) 


for plane impinging jets. 


As mentioned before, the static pressure in the impingement 
region is greater than the ambient pressure and Figure 4.7(a) shows 
typical pressure profiles in region II. It is seen that for any x 
station, the pressure decreases continuously from a maximum value of 
Po. 8 the axis as r increases. The pressure profiles are replotted 
in a dimensionless form in Figure 4.8 with p/P, versus m0 
nes being the value of r where p= > Po: It is seen from Figure 4.8 


that the pressure profiles are similar and well described by the equation: 


P = exp (- 0.693 n = (4.3) 
Be P 
where nN = mb Figure 4.7(b) shows two typical pressure profiles on 


the wall. It was found that the wall pressures are also well described 


by Equation 4.3. 


In runs 4 and 5, the wall shear stress was measured in the 
AROS 2 
impingement region and Figure 4.9 shows typical variations of u, , where 


u, is the shear velocity, Vt Ip. From Figure 4.9 it is seen that i 


* 
increases with r up to a maximum value of ee and then decreases with 
any further increase in r. When T, was replotted with ae versus 
r/H as shown in Figure 4.10, it was found that the results of runs 4 and 
5 and those of Bradshaw and Love (10) are described by one curve for 


r/H less than about 0.2. This aspect is discussed in greater detail 


later. 
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FIGURE 4.9 WALL SHEAR STRESS VARIATIONS 
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4.3 Analysis of Results 


4.3.1 Impingement Region 


For the free circular jet it is well known (28), (29),(36) that 


one could write: 


u £, (,.p.x) (4.4.a) 


cf 


I 


b 


ite £5. Q1.b0,%) (4.4.b) 


where M, is the momentum flux from the nozzle, equal to O7 ac uy 5 
Using the principles of 'Dimensional Analysis’ one could reduce Equation 


4.4.a to the form of Equation 4.2: 


For the impinging circular jet, one could write: 

u = f., (MP >%>H) C45 Fay) 
b = fy, (M,»0>%5H) (4.55b) 
Equation 4.5.a could be reduced to the form: 


€ 


U 
fo) 


=u f (>) (4.6) 


Au) 
On 
m 


Combining Equations 4.2 and 4.6: 


£. 'G/H) 
oe igs a ee (4.7) 
Ung C, G7 


The experimental results shown in Figure 4.11 agree with the formulation 


in Equation 4.7. The beginning of the impingement region could be 
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arbitrarily fixed at that value of x/H where u/uo¢ is ‘equal to; say, 
0.98 (or some other convenient value). The start of region II could in 
principle be fixed from pressure considerations also; but it was found 
that, since the pressure data scatter more than the velocity data, the 


velocity criterion is superior. 


The functional form of Equation 4.7 is evaluated below. In the 
impingement region, from Figure 4.7, it was found that ~. (uv) does 
not depend upon x for small values of r. That is, in the neighbour- 


hood of the axis of the jet: 
(uv) = By @) (4.8) 


which could be reduced to the form: 


I! 
g, (x) = = & (rw) (4.9) 
and 85 (o) = E 2 (wv) = constant (4.10) 
Yr f=" 0 


Using the equation of continuity, Equation 4.10 can be re-written as: 


pote dibint ics Vite eile ie uk (4.11) 
dx 2 2 
where K is a constant. Integrating: 
2 
it ae Sek + (4.10) 
e 1 


where KR is an integration constant. Using the formulation of Equa- 


tion 4.6, Equation 4.12 could be written as: 


Cc jie oe Nw Casa 
T oe ty Nia =F &) 8 Ky (4.13) 
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[2KH a x 
where K, = ae ; Ky = oH 7 1 and € = tl Combining Equations 


fo2 and 4,132 


(4.14) 


The present experimental results and those of Tani and Komatsu 
(15) were used to evaluate K, and Ky. Teewas founds that Ky = 3580 


and KS = 10. With these values, Equation 4.14 becomes: 


ob aBS31T0 all = 6) (4.15) 


This is shown plotted in Figure 4.11, along with the experimental obser- 
vations. This curve intersects the US line at x/H = 0.86 and 

this value of x/H is conveniently taken to denote the end of region I 
and the beginning of region II. Recalling Equation 4.13, it is easy to 


verify that with a value of K, = 0 the gradient of u. at the wall wili 


3 
be infinite. This is physically unacceptable and, therefore, it must be 
understood that Equation 4.15 cannot apply very near the wall. Rather, 
one could expect that the irrotational variation (a, « (H - x)) would be 
more reliable when x approaches H. The available data are not suf- 
ficiently close to the wall to permit verification of these considera- 


tions, however, they do indicate that an irrotational range would not 


extend further than about 2% of H above the wall. 


For the impingement region, considering the axial pressure 
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excess: p = 8, (MOH, x) (4.16) 


Cc 


Equation 4.16 could be reduced to the form: 


Pp 
Hee 
—~~ @° = ge (x/H) (4.17) 
2— ‘d 5 
pU 
oO é 
and further: 
P H 
= ang wg (1 .O)) =i constant, = .K (4.18) 
Zeid 5 Pp 
OU 
fe) 
where Ps is the stagnation pressure. Combining Equations 4.17 and 
4.18: 
Be 4g Fae 
Thies Sierra (x/H) (4.19) 
s P 


The experimental observations shown in Figure 4.12 support the functional 


form of Equation 4.19. 


Regarding the constant nike the present experimental results 
give an rsa value of 24.7. Previous works indicate values of 24.2, 
50.2, 24.0 and 26.1 (Refs. 10, 11, 22 and 23 respectively). Taking an 
overall average experimental value of x as 26, the stagnation pres- 


sure is given by: 


i (4.20) 


s H 
= 26/ (> 


Further, the length scale for pressure could be written as: 


. 5B. 
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The experimental results plotted in Figure 4.13 support this formulation. 
The value of 8 at the wall is apo = 0.078 H. Since the wall pressure 
profiles agreed with Equation 4.3, one could now transform this equation 


to the more convenient form: 


p/P, = exp (- Teen (4.22) 


where ) = r/H. This is compared with experimental results in Figure 
4.14 where it is seen that P,, approaches zero at r = 0.22 H, which 
could be conveniently taken as the end of the impingement region. 
Equating the wall pressure integral to the initial momentum of the jet, 
the theoretical value of a will be 28.5 which agrees with the average 


experimental value of 26. 


In this section a simple method is developed to evaluate the 
functional form of Equation 4.19. In the impingement region, assuming 


that the total axial pressure remains constant and equal to Pe 


Pr 


(see the inset of Figure 4.13), one could write: 


Di. = BY . (4.23) 


Using Equations 4.15 and 4.20, Equation 4.23 is reduced to the form: 


p 
ae = 10% -9 . (4.24) 


where the coefficients have been evaluated using the present experimental 
results. Figure 4.12 shows that Equation 4.24 is satisfactory for x/H 
greater than about 0.92. The experimental points of Bradshaw and Love 


(read from their contour plots) lie somewhat below the plot of Equation 
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FIGURE 4.13 LENGTH SCALE FOR PRESSURE PROFILES 
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FIGURE 4.14 DIMENSIONLESS WALL PRESSURE 
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4,24, 


The variation of the wall shear stress in the impingement 
region is predicted as follows. Consider the Reynolds equation of 


motion in the x-direction: 


(4.25) 


where os and tere the total, normal and tangential stresses re- 


spectively, i.e.: 


— a 1a we du 
Oo, = (pi pu ot Zu 
(4.26) 
fee es (laminar) + T (turbulent) 
rx rx rx 


where the prime denotes fluctuation and the overbar denotes time-averaged 
value. Assuming that near the wall the measured pressure p is approxi- 
mately equal to the sum of time-average static pressure and the normal 


stress, then: 


du 
EE aah + 2u ax (4.27) 


An order of magnitude analysis shows (see Appendix D) that the viscous 
term is negligible in comparison to p, hence Equation 4.25 can be 
finally simplified to: 
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At x =H, the inertial terms vanish, hence: 
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Buty = T and 
rx xr 
dnd oO 19 bid, 
ear (r,)] Es a) Sede ot 
x=H x=H 
Hence: 
is 
Ss Op 
ee | r (22) dr (4.30) 
x=H 
fe) 


Using the previous findings on pressure distribution, it can be shown 


(see Appendix D) that the wall shear stress is given by: 


2 
oa = 0.18 ee 9.43 X exp (- 114 > 
om 
(4.31) 
where the maximum shear stress, Lie is given by the expression: 
Alay a web 
FF ane 0.16 pu, @? (4.32) 


and itroccursvat throtcd. 14. 


Figure 4.10 shows that Equation 4.31 agrees with the available 
data so long as A < 0.22. Returning to the neglect of the viscous term 
from Equation 4.27, it may be argued that its vertical gradient may not 
be negligible in comparison to the vertical gradient of -p. If this 
were the case then the coefficient 0.16 in Equation 4.32 would not be 
constant, but it would depend upon Ro (as can also be verified by di- 
mensional analysis). However, the experimental results show that no 
dependence upon Ro exists, at least within the range 60,000 < R S 
180,000. (This includes the present data and those of Bradshaw and Love 


(10)). ‘he reason for this result can be traced to the strong turbulence 
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near the stagnation point, as was outlined in Chapter II. 


Considering the boundary layer at the wall, and the maximum 
value, Va» of the radial velocity, in any vertical, which occurs at 
the edge of the boundary layer, (z = 6), it is of interest to study its 


variation with r. Dimensional analysis will show that: 


Seo 
3 


a 
; a & f. (i) (45.33) 


The available data are shown plotted in this form in Figure 4.15, which 
confirms Equation 4.33. Assuming the Bernoulli equation to hold, it 


can be shown that: 


Vin H 
rie om (8 


ce) 


(4.34) 


Equation 3.34 is seen in Figure 4.15 to agree with the data only for the 
lower values of i (i.e., X4 < 0.03). For estimating vs the empiri- 
cal curve drawn through the data (solid line) is more reliable. For 
values of A less than about 0.03, Equation 4.34 applies, and recalling 
Equation 4.22 it is possible to show that in this range Vn can be 
approximated by: 


Ud 
(v,,) = (rz, “";) r (4.35) 
X-<.0.09 


This equation will be compatible with axisymmetric stagnation flow con- 


siderations if the constant a (Equation 2.1, Section 2.1) is taken as: 


a= V228K, =y (4.36) 
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Substituting this into Equation 2.5 (Chapter II) the boundary layer 


thickness near the stagnation point will be given by: 


&) = Bee (4.37) 
A < 0.03 RO 


The data of Hrycak et al (22) are re-plotted in a different form in 
Figure 4.16. It is seen that initially = remains approximately con- 
stant and equal to about 1.92 x ia. Using Ro = 54,000, Equation 

myore Will give 0197 x ieee i.e., about one-half of the measured value. 
This is not surprixing in view of previous considerations which point 

to a turbulent boundary layer. The work of Wolfshtein (37) on turbu- 
lent stagnation flow shows that if the boundary layer is turbulent, its 
thickness will be larger than the corresponding laminar value by a factor 


which depends on the intensity of turbulence. 


4.3.2 Wall Jet Region 


. No data were taken in the wall jet region. However, theoreti- 


cal considerations (see Chapter II) showed that: 


6, = C, r 
u Cc 
ee ee (4.38) 
TT r/d 
aS 
CoMiiecmpe “Sp Pop 2S 
pu, /2 


where r is to be measured from a virtual origin, and Ce should be in- 


dependent of r. Consider first the data of Bradshaw and Love (10). To 
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obtain C the data on 6, are re-plotted in the form 6,/H vs 


Za 2 
¥tH in Figure 4.17. Poreh, Tsuei and Cermak (17) gave the relation 
6,/H = 0.098 (r/He”, which describes their data well. The value of 
r/H was between 0.5 and 2.5. This relation is also plotted in Figure 
4,17 and is seen to give somewhat higher values than those of Bradshaw 
and Love. Finally, Hrycak et al (22) gave the relation 6,/d = K(r/a)o'?? 
for values of r/d ranging from about r to 40. The constant K was 
070945 for, H/d = 10 and 0.1025 for!uH/d = 20. | This relation is also 
plotted for the values H/d = 10 and H/d = 20 in Figure 4.17. The 
data show that they can be described well by the straight line: 

55 r 

a 0.0792 Gt 0.17) average (4.39) 


ave. C, = 0.0792 and the virtual origin is located 0.17 H behind the 


stagnation point. 


For the velocity scale, if the reciprocal of ue is plotted 
against r, then a straight line would result. The data of Bradshaw and 
Love are replotted in the form U6 / nm vs r/d in Figure 4.18, where the 


resulting straight line is described by the equation: 


u 
m L052 
v, SE (Bradshaw and Love) 


ie. s oe = 1.052 and the virtual origin for velocities is located 1.8 d 


after the stagnation point. Hrycak et al gave the equation u/0, = 


=1,12 


i464 (cfd) which is also plotted in Figure 4.18. Again, a straight 


line can fit the data reasonably well. This results in: 


un 0.94 


vu = Gla = it (Hrycak et ai) 
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Finally, Poreh, Tsuei and Cermak gave the relation u A/V = bre 
Pitan tes where K is the kinematic momentum of the jet at the nozzle 
(equal to dou 24) Figure 4.19 shows this relation re-plotted in 


the form 7K/u vs r/H. The resulting straight line has the equa- 


ELon: 
5 oe Lei5 
VK (= - 0.14) 
Which Gan also be written as: 
bn othe 02 
Us (x/d) - 0.14 (H/d) 


(Poreh, Tsuei and Cermak) 


The above re-analysis is summarized in Table 4.2 below. 


C Locus of 
eee u Virtual Origin 
72,000 b020 + 0.14 H 
to to 
325,000 31 


Hrycak et al 


TABLE 4.2 VELOCITY SCALE — SUMMARY 


Bradshaw and 
Love 


Poreh, Tsuei 
and Cermak 


Thus, the velocity scale can be finally written as: 
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Considering now the friction factor Ces Poreh, Tsuei and Cermak give 


the relation (which includes the data of Bradshaw and Love): 


Z 
TH Hat 3 


pK (ORS! (e/ay7*°? 


FE iC is independent of r, then T «= ae and a plot of: 
£ fe) : 


2 1. 
TH VK Or 5 aa 0) 
ok = vs r/H 


would result in a straight line. This is shown in Figure 4.20 where the 


drawn straight line has the equation: 


2 
TH a , 0.209 
ae (c/H - 0.14)° 


Recalling the relation found previously for us from the same authors, 


the above can also be written as: 


1, ees Bak . 0.209 
Gg Claas. (ai5)r 
mM 
OF eeinatiy : ome Ghote - (4.41) 
1F fe) 


which was derived for data lying in bee 72,000 sons < 325,000. It 
remains now to consider the question of where does the wall-jet region 
begin. The answer lies in finding the value of r beyond which the laws 
governing the development of the wall jet (Equation 4.38) become valid. 
Careful examination of the available data showed that the beginning of 
the wall-jet region (r = r) can be best defined from considerations of 


the velocity scale us According to previous findings (Equations 4.33 
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and 4.38) u, can be written as: 


oe (a 
B 


H 6 
; * ie fa @ (4.42) 


for all values of r/B: It follows that ty must be a fixed portion of 
By ite. , x i/H = const. Penation of the available data gives an 
average value of r/ = 0.36, even though some variation exists among 
different experiments. Recalling that the impingement region ends at 
r/H = 0.22, it is obvious that the range 0.22 < r/H < 0.36 represents a 


transition between regions II and III. 


4.4 Summary 


An analytical and experimental study of circular turbulent jets, 
impinging normally, has been presented in this chapter. The beginning of 
the impingement region has been located at x/H = 0.86 from axial 
velocity considerations, whereas the end of this region has been located 
at r/H = 0.22 from wall pressure considerations. A semi-empirical 
method has ees developed to predict axial velocity and pressure in the 
impingement Terion: Similarity of velocity and pressure profiles in this 
region has been explored and a satisfactory analytical method to predict 
wall shear stress has been developed. Existing data in the wall jet 
region have been re-analyzed and were found to confirm the simple-exponent 
theoretical considerations of Chapter II. Average values of the constants 
describing the variations of the length scale, the velocity scale and the 


wall shear stress, have been obtained and summarized. 
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CHAPTER V - AXISYMMETRIC IMPINGEMENT OF DEVELOPING JETS 
et ee lntroduction 


The impingement of fully developed jets was considered in the 
previous chapter. However, when the nozzle is placed at a small dis- 
tance from the wall, the physics of the problem is likely to change 
because in this case the jet is not fully developed when approaching the 
wall. The study of circular, not fully developed jets impinging normal- 
ly, is the objective of this chapter. The main quantitative difference 
between impingement at "small" and "large" heights is that since the 
latter is fully developed when it enters the impingement region, it can 
be assumed to originate at a point source of momentum, whereas the 
former consists of an inner core of uniform velocity equal to uy and 
an annular shear layer surrounding the core. This type of geometry 
depends on both Us and d and it cannot be described by the previous 
assumption. This study, therefore, was carried out at impingement 
heights less than the length of the potential core, i.e., H/d < 6. 
Quantitative flow characteristics are shown in the definition sketch, 
Figure 5.1. The flow issuing from the nozzle behaves initially as a 
free jet. At a certain distance from the nozzle, Xo the static 
pressure rises above atmospheric, i.e., wall effects begin to be felt. 
THUS. WeOTy -x> Xo pressure gradients exist in both directions, x and 
r. The gradient in the r-direction is set up so that the pressure falls 
from a maximum value on the centerline to ambient values when r_ be- 
comes sufficiently large. In the range 0 < x Sens the flow can be 


identified to be an axisymmetric shear layer, i.e., for r < Tyo us Ut 
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form dc> ty> ae 0, u approaching zero for r = Lye This region 
can, therefore, be appropriately called the free jet region. In the 


range xX, <x < H, u increases initially from a value of u. at the 


axis to a maximum of u ate ee Aer and then decreases again. In this 


m 1 
region the difference between the static pressure and the ambient pres- 
sure is generally greater than zero and is simply denoted by p. 
Furthermore, the radial component of velocity, v, becomes comparable to 
the axial component u, in contrast with the free jet region where 
Viis<ou. It is appropriate, therefore, to call this the impingement 
region, 


2 
Ou 


: : e 
lie 3 Pp. is the centerline pressure then Po ne 5 = const = 


Oo . : 2 a 
» 1.e., the total axial pressure remains constant throughout, which 


indicates inviscid conditions near the centerline. The wall pressure, 
Ee has a maximum of Ps at the axis and decreases to zero a few nozzle 
diameters from the stagnation point. Beyond this distance the flow soon 


begins to behave as a radial wall jet. 


5.2 Experiments 


Details of the experimental set-up and measuring techniques 
have been described in Chapter III. A nozzle with inside diameter of 


0.923 in. was used throughout. It was initially intended to carry out 


< 6. However, it was found that 


eu] oo 


experimentation in the range 0 < 
when H/d was less than about 1.1, some static pressure was transmitted 
from the wall to the nozzle, its relative magnitude increasing with de- 


creasing H/d. ibe Pr is the total pressure at the nozzle and Po is 
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Z 
the static pressure, then Pr = Py as pu, fi2:. The relative static pres- 
pu. /2 
sure is Po/ Pr =] - aa . This quantity decreases from 1.0 at 
It 


H/d = 0 to about 0.01 at H/d =1.1. To avoid unnecessary complica- 
tions only values of H/d greater than 1.1 were considered. Further- 
more, it was found that the stagnation pressure Po» was equal to pu,” /2 
so long as H/d was less than about 5.5. For H/d > 5.5, ee I 

pu,’ /2, indicating that diffusion effects have penetrated to the center- 
line of the jet. Thus, only the range 1.1 < H/d < 5.5 has been consi- 
dered in this study. For H/d > 6 to 10 the flow is that of large 
impingement height as described in Chapter IV. A transition between 
large and small impingement heights exists in the range 5.5 < H/d < &) 
where Ey is an, unspecified as yet, upper limit. Determination of 


aT can be achieved from considerations of stagnation pressure. This is 


given by: 


—— = (51) 
ou §/2 
‘i Eg 5 . e1 <4 (Chapter IV) 
(H/d) 
ae ou.” ie 
Let § = H/d and pecs P,/ 5 An equation for p, can be found 


in the range 5.5 ee < = using a second degree polynomial approxi- 
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: 5 ta, &- 5.5) +a, E - ess (5.3) 


it is possible to determine aye a and aT from the conditions 


Rae 


eee This procedure gives af ie a, = dn hee 030325, f BO. 3s 


With these, Equation 5.3 becomes: 


1 -— 0.0325 (&- Sosy (5.4) 


~ | 
i] 


and it is shown plotted in Figure 5.2, in the range 5.5 < H/d caO-o; 
which constitutes the transition between small and large impingement 
heights. Detailed measurements of the flow field were made in the first 


five runs whose particulars are shown below. 
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TABLE 5.1 RANGE OF EXPERIMENTS 


Particulars of supplementary experiments are indicated on figures where 


these data are plotted. 
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5.3 Experimental Results and Analysis 


3-3<L Free Jet Region 


Typical velocity distributions in this region are shown in 


Paeure o.3. Lt is geen that for +f Sti u = Us: From continuity 


eee 0. The v- 


requirements it is easy to show that for r<r 
component is also shown plotted in this figure. It is seen that the 
above reasoning is supported by the data. For r > Tyo u decreases 


to zero, while v takes on non-zero values and in general v << u. 


These profiles are replotted in Figure 5.4 in the form =o vs mn. where 


r- ry = 1 
ne = ar andi) b= as thelvakue of  xyi— ty for which u = 5 U ma. che 


conventional function exp(- .693 ne) applies as is shown in Figure 5.4. 


The range r<r is called the potential core, whereas the range 


1 


r 24) is the shear layer. Annular shear layers have been studied 
earlier (see (28),(29) for comprehensive accounts). Using the integral 
momentum equation and the entrainment data of Hill (41), a simple analy- 
sis (see Appendix C) showed that ry and b vary almost linearly 
with x, a result that is supported by available data. The analysis of 
Squire and Trouncer (42) yields complex expressions for vy) and b 


which, however, can again be shown to be approximately linear. 


The present data on ry and b are plotted in Figures 5.5 
and 5.6 where they follow straight lines up to a certain value Fe x/ds 
their variation Stedetes beyond this value, owing to the presence of the 
wall. Theoretically, when x =0, b=0 and r,/d = 0.5. However, 
the data show that this is not the case and only if the x co-ordinate 


is adjusted by what is commonly known as the virtual origin correction, 
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would the above conditions be satisfied. 


Thus in the free jet region, vy) and b are given by: 


ys x 
7 Bee 0.39 - 0.065 Ze G5; 5) 
b ao 
= he Sua. dOnt 0.087 (5.6) 
d d 
ane virtual ‘origin is located at x =~ 1.32d so that, if epi 


ao (31320) st 1e32d5" then: 


. Oe nes, (5. 6A) 


re 
According to Equation 5.5, the length of the potential core will be 6d, 
a value that is commonly accepted for circular jets. Theoretical values 


for the slopes of rj/d and b/d are 0.077 and 0.100 respectively (see 


Appendix C), which compare favourably with the experimental values. 


The outer limit r of the jet could be defined as that value 


2 
of r where u =o. Using Figure 5.4, one could deduce that rH = 4) aa 
2.50. The end of the free jet region could be defined as the value of 


x where r and b begin to deviate from the linear variation. 


1 
Figures 5.5 and 5.6 indicate that these limits coincide for each run as 


would be expected. However, this limit is much better defined for r 


than for b. 


5.3.2 Impingement Region 


Beyond the value x = Xo» the flow is affected by the wall. 


This is generally accompanied by an increase in static pressure above 
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ambient, increase in the v-component of velocity and increase in ty: 
At the same time the distribution of a in the core (r aE) is no 
longer uniform. Typical velocity distributions are shown in Figure 5.7. 
In the shear layer (r 2 4,) the similarity = = f{(m) is preserved as 
it is verified in Figure 5.8. A complete Sat of the flow field 
requires knowledge of the quantities, Us» UL» Tys OH; Pp, and wu and 


v in the core. Consider first the static pressure, p. One could 


generally write: 
p = S&é(p, Ux Geos. 2h ot) 


or, by dimensional analysis: 


_ ae ZX) He cr 
fo) 
further: 

2879. 2 (0, %, 8) = 2, 5 (5.8) 
ou 219 ne oan kA | 2d da } 

fe) 

P. H H 
—37 £, (0, 0,4) £, (0, a (5. 8A) 
ou #2 

fe) 


and the relative wall pressure is: 


Pa i 
po Gea = Ge 


3 ) (S579) 


A. | pm 


For large impingement heights, the function f is independent of H/d, 


; rants a nia aA 
as was shown earlier (Chapter IV). Figure 5.9 shows existing data plotted 
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in the form P,,/P, vs r/H for different values of H/d. For H/d = 8, 
P,,/P. is close to the average curve for large heights, whereas data of 
different authors (15, 22) for H/d = 4 yield practically identical 


curves, while they are considerably removed from the large height curve. 


Consider now Equation 5.8. In the range of interest, 1.1 < 
4 we 
0U 
H/d See oe 1, which suggests that unless the functions fi 
and f, are of special form, they should be independent of H/d. This 


Z 


expectation is also supported by momentum considerations. Equating the 


integral of the wall pressure to the initial momentum of the jet: 


i] 
1) 


2 
an fs p, dr m ite (5.10) 


which, using Equation 5.7 reduces to: 


Cc fee) 
RSs SOBs See 3 r FH wm or, 41 
fe) fe) 


which again suggests that the functions fy and f. should be indepen- 
dent of H/d unless they are of special form. The available data ieee 
the range of small impingement height are plotted in the form aye 
vs r/d in Figure 5.10. In spite of scatter, the data appear to define 


a single curve, thus verifying the above reasoning. The Gaussian varia- 


tion would give here only a fair approximation. 


On the basis of the above, it is reasonable to assume that the 


function’ £ defined in Equation 5.7, is independent of H/d. It fol- 


a ies 


lows that the static pressure p at any point in the flow field will be 
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given by: 
Po yt) (5.12) 
where an overbar has been used to denote dimensionless quantities, i.e., 


2 
pu. 
2 


Dp p/ , z=2/d, r-=r/d. The centerline pressure, Poe will 


then be: 


Mo) 
{ 


f(z, 0) = £, (z) (5503) 


This formulation is verified in Figure 5.11, where data for H/d varying 


between 2.21 and 5.17 are seen to define a single curve. 


Finally, the complete pressure field can be constructed, using 
the present data and Equation 5.12. This is shown in Figure 5.12, 
where average contours of equal pressure are presented. The physics of 
this problem is now clear. The effect of the wall (expressed in terms 
of increased static pressure) extends a certain distance above the wall. 
This distance does not depend on the height of impingement H but it is 
proportional to the nozzle diameter d. Thus, with reference to Figure 
5.11 it could be stated that the boundary between free jet and impinge- 
ment regions on the centerline is given by z/d =1.2 - 1.4. This is 
also supported further by Figures 5.5 and 5.6 where vy) and, b “are 
seen to first deviate from the linear variation at approximately the 
same value of z/d for all three runs (an average is ~ 1.25). The 
lateral extent of this region at the wall can be obtained Erom Eiger 
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Since the total pressure on the centerline is constant, one 


could write: 


e 
Q 

HI 

be 
! 
"S| 

Q 

i 


f(z) 


Pais #6 Verified in Figure 5713, for H/d faneine froms2.21 to 5.17. To 
derive an equation for u/U5> an empirical observation is utilized which 
states that the product uv varies very little with x near the center- 
tine,\t.e., for small values of r. This is shown in Figure 5.14 where 
uv is plotted against r for different values of x/d. It is seen 


that for r< 0.1 in. uv is almost independent of x (or z). Thus 


uv = g(r) then eee (ruv) = 8, (x) and L a (ruv) | g, (0) const. 


IGE) 


The above can be easily transformed (see Chapter IV) to: 


c i 2 
or: 
ve z "a 
Sheth = wee dt 
Ue 32) d Ko 


The constants k and k are evaluated to fit the data and the above 


becomes: 


lo 


P . 
Z = 
; = 1.07 ‘7 0.07 (5224) 


Equation 5.14 is shown in Figure 5.13 to describe the data adequately for 


0.1 < z/d < 0.8. For 2z/d > 0.8, there exists a smooth transition between 
u 
Equation 5.14 and the line Se. For very small values of z/d a 
fe) 
Square-root type of variation would give infinite gradient of u. This 
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is physically unacceptable. If it is assumed that the flow very near 


the wall is of the stagnation type, then uu. Ce. Matching this with 


du 
Equation 5.14, so that u. and ee are continuous at the matching 
point one finds that: 
Zz 
=S < ° = 
Ae 0.14 : u/U, 2(z/d) (5315) 


Consequently, the variation of Po can be found using the Bernoulli 


equation which applies on the axis. 


p 
For gid “= 0214 BS SC a Yee y, (ke 2 
7 vu 7/2 d 
as 
(5.16) 
Po 1 
For z/d > 0.14 08 = 
ws d 
peu, /2 


These are plotted in Figure 5.11 and are seen to describe the data ade- 


quately. The variation of u, can also be described by the simpler 
A Me : : > . 
relation Ta = 1.0 ¥z/d »with the understanding that this is not to be 
fe) 
applied very near the wall. 


Considering the maximum velocity uw. it is assumed that 


eles ae Vz/d. The constant of proportionality was evaluated as 


WeL5 (u, = 1.15 ud) and this is shown plotted in Figure 5.15 together 


with the data. A simpler relation can again be taken as u/U, = 


oo) ved. In order now to predict u and v in the core, it will be 
assumed that in the range 0 <r < ty. the flow is inviscid. This is 
supported by the constancy of the total pressure on the axis. At the 


Same tite, irrotationality cannot be assumed because the total pressure 
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does not remain constant in the radial direction. Under these conditions 
the solution of Tani and Komatsu (15) will be applicable. Neglecting 


terms which involve powers of r greater than 4, u is given by (Chap. 


EE) 
a RUN a AOE? + oe) (5.17) 
¢ 1 2 | ‘ 
. Melby. -& 1 (iv) 
with o, (2) = i (ue + C,) and , (2) = (Cou, +3 ue ) 
where Ci> C, are constants. 
ou ob 5 
eae 2g) + Adar 
; “ne du 
Using the condition (=) = 0, 
or r=ry) 
¢ 
2 2 al 
ey = oy é 65.18) 
2 
With this and some further manipulation it can be shown that: 
ue u 
= 4° (RB) (5.19) 
Wow 
m c 


with j= r/r,- 


The data are eee in the form suggested by Equation 5.19 
in Figure 5.16. It is seen immediately that Equation 5.19 is, strictly 
speaking, a poor approximation. However, if one is interested in finding 
u rather than u - uo Equation 5.19 becomes adequate because u - u. 


is about an order of magnitude less than u. To demonstrate this, let 


the subscripts p and a denote predicted and actual values respectively, 
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Then the error will be in absolute value: 


It, - %, 


le] = 0.15 4 0.159, 


The maximum value of le| is plotted against beneath the main graph 


of Figure 5.16 and it is seen to remain less than about 5%. To predict 


Ty) consider again Equation 5.18: 
Ww 
ae) ee ee On Be es 
1 “a2 1 (iv) 
2 Cou. + BUS 
With u_ <2 vz/d as was seen earlier, ue and hat can be obtained. 


If one tries to evaluate Cc and C, so that the above relation fits 


the variation of r it will be necessary to use very large values of 
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written as: 


r/d = Cera. (5.20) 


where now C depends on H/d, and will be discussed later. This rela- 
tion is verified in Figure 5.5 with C evaluated empirically. As it 
could have been expected, Equation 5.20 intersects the free-jet linear 
variation at a distance zy from the wall such that zi/d is indepen- 
dent of H/d. For runs 1,°2, 5 respectively, zo/d et. 2), 1, ho and 
1.34. Taking an average value of z/d = 1.25¢,the value of C can be 


predicted as follows. 


Ty H-z 
cee 0.39 =90.,065. ¢ = ) (free jet region) 
AED gee 
fo) 
is C 
igen ae (impingement region) 
"vz Jd 


Equating and simplifying with z/d chal WAS 


C = 0.50 - 0.069 t (5.21) 


This equation is plotted in Figure 5.17 together with the empirical 
values of C and it is seen to be adequate. The radius of the core, 


3 was determined by means of the relation: 


alee 
Go. esa 
Ss Maat eN) 
1V 
ohare ri 8 “c 


1 ; ; 
and by observing that this ratio behaves almost as a This conclusion 
Cc 


can be obtained formally by requiring ue = const. This will yield: 


' J : | i 5 i ; ¥ ; a 
* Vl | : ; | i ¥ oy f a i r ye Sy vba tt 


pereee 2i¢T roam bayonlt ot | et bas ha ne hep 3 vo0 4 
roe nchisatakaae hsaawsevia a “Hat hye. pumgt® ak baranaoy ak 
xeant I satomndh ‘ait adonmensed, of: 4 nak taupe dhossoqan anes ype 

‘~asigphnt et’ Bh, ry aad “oe “Ei et x “ eonmaatd, & 4 fh 

wre oe. Y cm * BA,S Weteubenegans 2 4% af saz tot NY 
adinso 0 ae aul wie 125. i= \¢" ihe ent peat wa staat 
ca in - easel ae t 


ir at te 


4 «at eh | 2 aa 
oat sot soy ge ae ae 9 AR) 


. V a 7 7 
a Vi 
; - re) 


(nolge4 ssdmegnkget) 


; abi © WN a Wie @ateia tanta, bas gakaae 


L s 
= a 3 ‘ | 
=| ry 


ae fo ae a L emalbecoumialice . 


enn 


ane gor stipes igs Anish wirsagae St.8 it ry poate? andere LT ~ 
li 
eRe, oi | we apthe? ait isttnuypobs ad 0+ wepu ek 3 3h bos 3 to souluy an 
| Pd sehen as dts. Ya eneeut..ed bontarieasb esw) sagt" ; ' 
[ aa ror, oe - @ as ey ar, 2 on . - ine te. 

HW i ieee > | ; F » £ via 


ae aadl ~ git fw, 
wr ee 
nA Tg ee 


tre tavtonet: SLAP > ab sone angel ‘ okdew aids na nidbieaio “e has 
- 
hae iby aL » yay “ si. aioe wi ‘atsicsot pontalso. od ng 


, 
“s} 


2.0 


FIGURE 5.17 


3.0 4.0 5.0 
H/d 


VARIATION OF C WITH H/d 


118. 


Gl 


eeu 


j ¢ ire he, aig > uy . 7 ; 1 } F 
: 2 ale A i : Sey [tia S ‘ e : 
‘ph feta 9 Fo WoftAraay -TL.¢ Savore: 
fost - 3 ; > ae 
a is J 
j ot Lt , 
+ : a a 
a ‘ - 4 
m7 - : , 
: 
vy i ? 
A a ) gig t 
i 
; _ asf Lf 
a 1 ; 7 os a iy 
tha _ x , a 
aa “ 
. oe 2 Ja , 7 Ae 
ne ¢ 7 : 


ESO: 


Za, 


2 (5.22) 


4 2 
as + a, +a 


The coefficients Oo Qs a can be determined from the following boun- 
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dary conditions: 


2 fe 
ay Tie 
fo) 
e233) 
u d(u /U_) 
Zi pes : ECr Car | 
aL ee : v, 1.0 and meer cm 0 
These will give: 
u 
e. ) zjd z/d 
yepect Tyi6 42) ae (52) 


This relation is obviously superior to previous equations for predicting 
Boe. as shown in Figure 5.13. However, for purposes of predicting ry 
and for later analytical calculations, it is much simpler to work with 
the square-root relation. It remains now to predict the radial component 
of vetocity - “vi. bet Au Hetie Be and (A) = 7 (2 ~ vee By 


Equation 5.19: 
u = us 4 Au oO) (3425) 


Using this relation and integrating the equation of continuity, it can 


be shown that (Appendix D, Section D.2): 
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where f£'(z) = u/U, and F(A) is a function defined by: 
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F(A) = A + BOF, (A) + 28, ())A7 


With 6 =sAq°/ar2e*0cTS’ and: 
Coce 


r r 
FO) - fora aE OO - [re 
Oo ) 
EE Vi is the value of v at r= Th> then: 
v,/0, AL — — 
=-=>f'(z) F(1) = function of z only (S.27)) 
r,/d 2 


This conclusion is verified in Figure 5.18, where the data define a 
single curve while H/d is varied from 2.21 to 5.17. If our previous 
findings with respect to £'(z) and F(1) are substituted in Equation 5.2/, 


then: 


ea Ge mana (5,28) 
Vz - 0.07 

This is seen in Figure 5.18 to be in fair agreement with the data only 
for “2 < Ole. The discrepancy is believed to be due, mainly, to the 
fact that differentiation of f(z) is involved. Even though our approxi- 
mate prediction of He is reasonably good when applied to ul as such, 
it could be in considerable error when prediction of the gradient is 
attempted. Neperene eee the analysis has shown the correct dimension- 
less grouping of parameters, so that v, can be predicted using the 


empirical curve drawn through the data in Figure 5.18. 


Furthermore, dividing Equation 5.26 with Equation 5.27, the 


Vabiation. of Vv) is obtained as’ 
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Nea = ray = function of i only (S220) 


This formulation is verified in Figure 5.19 where the data are seen to 
define a single curve. Using our previous definition of F(A), it is 
possible to show that (Appendix D, Section D.2) Equation 5.29 can also 
be written as: 
y) 
OE er eT NO TS aa) (5.30) 


Vy ah 


This is plotted in Figure 5.19 and agrees generally with the data, even 
though it yields somewhat higher values. Considering the approximations 
involved in this derivation, the agreement is remarkable. If the flow 
were not only inviscid but also irrotational, then v would be propor- 
tioned sto. re, ige., v/v, = v/r,. This is also shown in Figure 5.19 


for comparison. 


With these results, it is possible to predict the velocity 
field within the core. The assumption of inviscid flow and the use of 
the truncated power series expansion for u has been fruitful, even 
though this flow model may not be realistic throughout the core. For 
example, careful examination of the full equations of motion will reveal 
that some shear stress is likely to develop within the core in the 
impingement region. At the same time, this flow model is not entirely 
free of inconsistency. Applying Equation 5.17 at r = T1> this would 
lead to Oe al ee L/u, which contradicts the relation Oe gs 


that was established before. Yet the former expression yields reasonable 


answers for not very small vaiues of z/d. It may be stated in conclu- 
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sion that the inviscid flow model should be regarded as a working hypo- 
thesis rather than as an accurate physical representation. It is 
satisfactory in so far as it has made it possible to analyse the experi- 


mental results in a systematic manner. 


Considering the shear layer in the impingement region, it is 
evident that u can be predicted using earlier similarity considerations. 
For this, it is necessary to know the scales us and b. The velocity 
scale us is given in Figure 5.15. In Figure 5.6 it is seen that the 
length scale b deviates from the linear law in the impingement region. 
It was shown earlier than rj/d = ce), i.e., the two variables 4H/d 
and z/d upon which rj/d depends, are separable. Assuming that a 


similar relation holds for b, then: 


2 = function of z/d only 
Oo 


with bo being the value of b at gz =s= Zo (i.e., the beginning of the 
impingement region). Figure 5.20 verifies this conclusion. The value 
of b, can be found from Equation 5.6 by setting x/d = (H - zid/d, 


which gives: 


b/d = 0.087(H/d) + 0.006 (5.31) 


According to this, the predicted values of b/d f£Or-runs J12)2,. 0° Les= 
pectively will be 0.198, 0.358, 0.456, whereas the corresponding experi- 


mental values are 0.200, 0.363, 0.445. 


After impingement, the flow forms a boundary layer on the wall, 


beginning at the stagnation point. Denoting the maximum value of v, at 
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the edge of the boundary layer, by v3 and applying the Bernoulli equa- 


tion: 
ve} D 
u, Als (5.32) 
pu, We 


Present data on Va are plotted in Figure 5.21 together with data of 
other authors (15),(43). It is seen that Equation 5.32 is adequate for 
r/d <1.0. Shortly after, the wall jet region begins. It was found 
@ariier, (Equation 5.15) that for z/d’< 0.14, uw. varies linearly with 
z with a constant of proportionality (see Section 2.1, Chapter II), 

2a = 2|u,|/d. That is to say, if ideal fluid stagnation flow occurs in 


this range, the constant a is equal to lu|/4, and therefore: 


(5.33) 


This is also plotted in Figure 5.21 and is seen to be adequate only for 
td <0. ta, Hence, stagnation flow may occur within a small neighbour- 
hood of the stagnation point, with radius 0.14d. Furthermore, using 


a= Ju, |/a and Equation 2.5 (Chapter II), the boundary layer thickness 


will be: 
fo = 8 (5.34) 
vR | 


Hrycak's et al (22) data for H/d =) 3. and Ro = 54,000, indicate 
& ~ .007.-°AE this yalue of R.. | Equation 5.34 will yield 
d fo) 
r/d-< 1.0 
6 


Be 0.00835 which is again (see also Chapter IV) about one-half of the 


measured value. 


re ae ie Tyee 
. j i P| iM 


A 
y AW 
Whe 


I 7 
¢ + : Bp © + 
«O5E \ oO H : edt To) 
Jia s 


-aupe Etheantse ale 2 wee - 4 o 


(i0..2) 
: bares re : 


ab 4 sc a ot 
suotiasads 2 


(Etre) i ag 5 


~ 


. | r’ 
+0 tai essay és fa fame on mt 
~~ on hl w! onto, Anoop deans chy 


antau- mgasioes “shat cuban inky (aa oeininansh 83 » > 
waphAotitt say udinwos oot” , (LK ‘suaigai) es Loman bas »\y Na os 
a Ale a ae 7 


VT pa he - 88 e 
ya , re Nii 4 “* ; 
(6£.0) | Ol ia m s 


sinless re ee M las ee WN "Gah sa (58) te a8 aor 


bint, Lite | 6.2 obdnops fi ) 7. ents 3A +010, | wa st 


12. 


Consider now the wall shear stress in the impingement region. 


Applying the same reasoning as in Chapter IV, it is possible to relate 


Ue to the local pressure gradient, i.e.: 


16 


meee op 
TG - ee 


O 


Since p does not depend upon H/d, TY 


this parameter, i.e.: 


18 
oO 


ee ee) 
2 d 
pu, 2" 


It follows that the maximum vaiue of we 


i = const = C 
ou 219 fm 
fo) 
Dividing: 
T 
fe) iG 
=< = h,@ 
om 


This is verified in Figure 5.22, whereas experimental values of C 


are tabulated below. 


G35) 


must also be independent of 


(6.36) 


Wiloeeper 


(5.37) 


(5.38) 
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It is seen that Cen is indeed independent of H/d, however, a slight 


dependence upon R seems to be present. 


From stagnation flow considerations (Chapter II), it was shown 
that initially gi es This was found to be the case for large impinge- 
ment heights (Chapter IV). However, this is not so in the present case, 
as Figure 5.22 shows, even though a stronger expectation of linearity 
exists here because there are better chances of the boundary layer being 
laminar at small impingement heights. This leads to the conclusion that 
shear stress measurements near the stagnation point are suspect. One 
source of error is the pressure gradient, which was estimated from 
Patel” s' (40) criteria to cause errors larger than 64 for r/d <.l.5.5 A 
more significant source of error near the stagnation point is the small 
thickness of the boundary layer. This was found earlier to be of the 
Orcer 0.0l7d, 1.6.5 since €@-= .9239in., 6°= 0.016 in. Considering 
that the outside diameter of the Preston tube was 0.049 in., its center 
was located outside the boundary layer and, therefore, the present mea- 
surements near the stagnation point indicated dynamic heads rather than 
shear stresses. This is supported by the shape of the shear stress 
curve at very small r/d, which is roughly of the form is oe Since 
for small values Giaar, ie ey (Figure 5.21) the dynamic head will be 
pv /2 « ee The data of Hrycak et al (22) indicate that the center 
of the Preston tube will be inside the boundary layer for r/d > 2.0. 
It follows that the present data on wall shear stress near the stagna- 
tion point must be regarded as meaningless. Yet, the data on maximum 


shear stress Tom are believed to be reasonably accurate. This expec- 
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tation is supported by the following consideration. The skin friction 
coefficient, Cow has an average value of 0.0069 for the present case, 


whereas for large impingement heights Cen = 0.32/(4)*. These two 
variations will become equal at H/d ~ 6.8 which is close to the corres- 
ponding value of 7.2 where stagnation pressures for large and small 
heights become equal. If the theoretical considerations of Appendix D 
are applied in this case assuming the pressure to be gaussian, the 
following relation will be obtained: 
TO = 1.47 eae 5 al. B(r/d)*, 7 eee 
(5389) 
This relation is plotted in Figure 5.22 and, in view of the above consi- 


derations, it is believed to be more reliable than the data for r/d < 0.8. 


5.3.3 Wall Jet Region 

In this region the flow assumes the pattern commonly known as 
the wall jet. In view of previous considerations (Chapter II) the 
analysis will be merely a demonstration of the fact that this region is 
independent of the parameter H/d. This can be shown formally (Chapter 
IV) by assuming that it is influenced by a single kinematic momentum 


term, mdU Ph Thus: 


6, = Cor 

u oO 

pln (5.40) 
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Consider first the velocity scale un The re-analyzed data of Hrycak 
et al (22), have shown (Chapter IV) iat uf, = 0.94/(5 - 130909) for 
Values om H/d equal to 4, 7; 10 and 20. This includes both small 

and large impingement heights and supports the independent development 
of the wall jet. For the case H/d = 4, the wall jet region begun at 
ri/d = 2.0. The results of Yakovlevskii and Krasheninnikov (27) for 
H/d = 3.5 and 10, also showed no influence of H/d (see also Chapter 


VEL). 


Considering the length scale Sy5 Hrycak et al (22) proposed 


the relation (directly from their data): 


5 


g/d = K(r/aye” (5.41) 
where K was given as: 
K = 0.0864 + 0.00081 (H/d) (Se42') 


These were obtained from different experiments using H/d = 2, 4, 10 and 
20. Within the range 1.1 < H/d < 5.5, K_ changes from 0.0873 to 
0.0908, hence it is permissible to accept an average value of 0.089 and 


assume K = const, i.e.: 


6/4 = 0.089 ees (5.43) 


This relation is shown plotted in Figure 5.23 for r/d = 4 to 40 which 
represents the range of experimentation in (22). It is seen that Equa- 


tion 5.43 can be represented by the linear relation 5, =, 0/6555" dani 


Cc, = 0.0765 (5.44) 
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The data of (27) show the same trend, but their value of C is slightly 


2 
lower (Cc, = 5075) These values of C, are also very close to the 
average value obtained for large heights, C, = .0792 (Chapter IV). 


2 


Present data on shear stress are shown in Figure 5.24, where 


u,,/U, is plotted against r/d with H/d varying between 1.22 and 4.55 


and for Ro equal to 13,800 and 55,300. It is seen that the data 
define a single curve. Dependence upon Ro cannot be ignored, however, 


the range of Ro considered is rather small and this dependence cannot 


be detected in these data. By Equations 5.40, u, should vary as r , 
and the following relation was found adequate: 


ae te 0.0794 


US 1 ede OLS 


(5.45) 


Assuming an average value of GC. = 1.0 (Chapter IV) and neglecting the 
virtual origin correction, the local skin friction factor Cy will be 
2(.0794/1.0)7, i.e.: 


CG. = 0.0126"; 31,800 <R_ < 55,300 (5.46) 


Summarizing this section, it has heen shown that the wall jet region 
begins at ri/d = 2.0, and it develops in the same fashion, as in the 
case of large impingement heights. The only difference lies in the 
values of to and loci of virtual origins which are now related to d 
rather than to H. A detailed analysis in this region will be presen- 
ted later (Chapter VII) where the general problem of oblique impingement 


is dealt with. 
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5.4 Summary 


The impingement of developing jets has been considered in this 
chapter. It was first shown that a given value of H/d can be consi- 
dered "small" in the interval [1.1, 5.5] whereas it can be considered 
large if it is greater than about 8.3. The range 5.5 to 8.3 is transi- 
tional. Further, it was established that the static pressure, including 
the wall pressure, is not affected by H. Wall effects extend about 


1.25d above the wall and about 1.4d along the wall, regardless of H/d. 


It was possible to predict the velocity field in the impinge- 
ment region, using the inviscid-rotational model as a working hypothesis. 
Knowledge of the free-jet region is necessary in order to provide infor- 
mation regarding the boundary conditions for the impingement region. 

An approximate analysis was developed to predict the annular shear layer 


which yields very simple expressions. 


The wall jet region was confirmed to develop independently of 
H/d, its properties being almost identical to those of large impinge- 


ment heights. 
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CHAPTER VI - OBLIQUE IMPINGEMENT OF PLANE TURBULENT JETS 
6.1 Introduction 


The case of a plane turbulent jet impinging normally has been 
studied by the author in detail (25),(26). A more general version is 
obtained if the angle at impingement, 6, defined to be the angle between 
the jet centerline and the wall, is allowed to vary. At present, there 
is relatively little experimental and theoretical material regarding 
oblique impingement. Schauer and Eustis (14) have studied mainly the 
wall jet region. Kamoi and Tanaka (24) were mostly concerned with tur- 
bulence characteristics but some wall pressure measurements were 
included. Mathieu (12) deals with the wall jet region for small angles 
of impingement (¢ < 30°). The present study covers the range 20° < $6 
< 90°, which is the most important from the practical point of view. In 
the impingement region, the objectives were to study wall pressure and 
shear stress, maximum velocity at the edge of the wall boundary layer 
and to define the limits of this region. No measurements were performed 
in the wall jet region, however, a simple theoretical treatment is 
presented and is compared with experimental data of previous investiga- 


tors. 


6.2 Experiments 


The experimental set-up has been described in Chapter III. 
Wall shear stress was obtained using a Preston tube with outside dia- 
meter of 0.046 in. and utilizing the calibration curves of Patel (40), 


as has been explained in Chapter III. A total of 15 experiments were 
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performed, of which the particulars are summarized in Table 6.1. The 


symbols appearing in this Table are defined in Figure 6.1. 


Free Jet 
Region 


Region 
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TABLE 6.1 RANGE OF EXPERIMENTS 
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6.3 Results and Analysis 


6.3.1 Free Jet Region 


Quantities Measured 


Impingement 


For normal impingement (¢ = 90°), the boundary between free- 


jet and impingement regions is located at & ~ 0.7H (26). 
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"large" impingement height, the centerline velocity, u,» can be written 
2 


asi ou, = f(o, pu, Gt AP lee) 2 It follows that: 
Se & 


For a free-jet, the corresponding value of axial velocity, Woe? is given 


by: 
cf wee 
ves veld 
or: 
Yok H a 
| aS (6.2) 
fe) v&/H 
where Cc. is a constant. Dividing Equation 6.1 with 6.2: 
u 
RC vee & 5 
2 = f., 7? >) (6.3) 
cf 


ig EO is the value of €&€ giving the boundary between free-jet and 
impingement regions, then one could define 5 such that u/uo¢ 


equals, say, 0.98. Hence: 


. BS 
ieee. f,G> >) 
which implies that: 
Oo — 
= = £4(6) (6.4) 


u 
All data on u, are plotted in the form poe ve J6/H. ay wigure 6,2. 
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It is seen, that for €/H <0.70, the data follow the free-jet variation 
(Equation 6.2) where the constant is is taken as 2.50. It follows, 


that O/H is independent of 9, i.e.: 


3 
oO 
am 80.20 (6.5) 


Schauer and Eustis (14) proposed that the distance from the stagnation 
point to the boundary between these regions is a fixed portion of H. 
This was verified only for the case @$ = 90° and it obviously contra- 
dicts Equation 6.5. The discrepancy is believed to be due to the strong 


assumptions introduced, as was outlined in Chapter II. 


6.3.2 Impingement Region 


6.3.2.1 Wall Pressure 

Distributions of wall pressure for different angles and heights 
of impingement are sherd in Figure 6.3. The stagnation point, s, which 
corresponds to the location of maximum wall pressure, is in general re- 
moved from the intersection of wall and jet centerline, by a distance s. 
This is not hard to explain if one considers first the flow of an ideal 
fluid. In this case, it is known (44) that a stagnation streamline (see 
Figure 6.1) must intersect a solid boundary at right angles. Assuming 
that stagnation-type flow occurs at least in a small neighbourhood of 
the stagnation point for flow Spices fluids, as pointed out in (36), 
then the stagnation streamline will intersect the wall at right angles, 
which causes the stagnation point to become offset. Figure 6.3 shows 
that the wall-pressure profiles are skew, the skewness increasing as 


@ decreases. Furthermore, for the angle ¢ = 30°, the pressure becomes 
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negative toward the side of negative x. This is believed to be due 

to the formation of a vortex. When the angle becomes sufficiently 
small, the boundaries of the free jet and the wall jet get so close at 
the side of negative x, that some of the fluid in the upper layers of 
the wall jet, which has very low momentum, is entrained into the free 
jet, thus completing a circulatory motion and forming a vortex. In view 
of this explanation, it is not surprising that the solution of the cor- 
responding potential flow problem (45) does not predict negative 


pressures, because it assumes the flow to be irrotational everywhere. 


It is now convenient to adopt the following convention. Re- 
ferring to Figure 6.1, let xy be distance along the wall, measured 
from the stagnation point, so that the positive axis x) points towards 
the right-hand side. In order to avoid use of negative values of xy 
it is observed that the negative axis x will be the same with the 
positive axis if the impingement angle is made 180° - >. Hence, by 
referring flow properties in the range x <0, to the case 180 - 9, 


it is possible to use always positive values of xy 


For pressure profiles that do not exhibit negative values, 


let b_ be the value of xy for which Pa 5 P (Pp. = stagnation 


s 
pressure). Then, Lf P,,/P. is plotted against Nn = x,/b a single 
curve results, as shown in Figure 6.4(a). Hence, P,/P. = g(n) and 


the function g(n) will be the same as that for ¢= 90°, which is given 


by (26): 
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Equation 6.6 is also plotted in Figure 6.4(a) and is seen to describe 
the data adequately for 1. <,1,6, being somewhat low for n2>1.6. 

This type of similarity has also been reported by Kamoi and Tanaka (24). 
However, it cannot apply when the profile takes on negative values. If 
p is the minimum value of the profile (p < 0), then Figure 6.4(b) shows 
that similarity is preserved if (P,, - P)/(, -p) is plotted against 


n= x,/b. Here b is defined to be the value of xy where Pi p= 


5 (Pp, - Pp). 


Complete prediction of the wall pressure requires, in addition, 
knowledge of the quantities Ps> b, > and s. With the same assump- 
tions as in the previous section, these quantities can be shown by 


dimensional analysis to be: 


Ep = = g,() (6.7) 
r= gy (4) (6.8) 
fig Cd | (6.9) 
= 2, (6) (6.10) 


The stagnation pressure, Ps? is shown to agree with the formulation of 


Equation 6.7 in Figure 6.5. Schauer and Eustis (14) gave the relation: 
P 
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which was derived by equating the integral of wall pressure to the normal 
component of momentum. However, the derivation was based on the assump- 
tion of symmetric wall pressure distribution, which does not allow for 
skewness or negative pressures for the lower values of @ and, therefore, 
it is only a crude approximation of reality as shown in Figure 6.5. 

Figure 6.6 shows the variation of p/P. with 9. It is seen that 
negative pressures do not occur in the range @ < 130°, whereas at 


wee 150>, Dp becomes a significant fraction of Ps: 


The variation of the relative length scale, b/H, is shown in 
Figure 6.7. (The present data were obtained for values of H/d equal 
to 45.5 and 68.2, whereas the data of Schauer and Eustis were obtained 
for the values of H/d equal to 20, 30 and 40. The data of Kamoi and 
Tanaka were taken at H/d = 12 and 16.) Gmethe zange 30° < ¢&< 90°, 
b/H is seen to change very slightly. However, for $ > 90° the varia- 
tion of b/H with @$ is more significant, and the following empirical 


relation was found to be in reasonable agreement with data: 


1 


. PS * (oe by" > os 90° (6.12) 
180-9 

Figure 6.8 shows the relative eccentricity, s/H, of the stagnation point 

as a function of >. Schach (45) has given the ideal flow solution for 

an impinging jet of finite width. The resulting expression for the ec- 

centricity s is of a very complex nature and is shown graphically in 

Figure 6.9, where Co denotes the width of the jet at an infinite dis- 


tance from the wall. In the range 30° < > < 90° this relation is 


almost linear, i.e.: 
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err whe 2 
q 2B15%4 90 , (6213) 


where $ is in degrees. 


In order to use this equation for the present results, it is 
required that the corresponding potential jet has the same initial mo- 
mentum and produces the same stagnation pressure as the real jet. 
Letting U_ be the velocity at infinity of the potential jet, then the 


following relations must be satisfied: 


Ca = pt ds = pv7/2 
OY p p > Be 0 
Using these conditions and the approximate expression for Pi (Equation 


6.11) it.can.be shown that: 


Bae 
. 90 
sin 


m|o 


The experimental data indicate that this relation is satisfactory, the 
constant of proportionality being equal to 0.194, i.e.: 
= 0.194(1 - 2)/sin 6 (6.14) 


= aa 
H 90 


With the above results, it is possible to predict the wall pressure up 
to the point where the minimum pressure, Ps occurs. An approximate pre- 
diction of pressure beyond this point can be derived by assuming that 
this pressure is essentially due to a point vortex located a distance a 
above the wall. If JT is the strength of the vortex, the wall pressure 


he at a distance 2 from the projection of the vortex on the wall, is 
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given by (see Appendix D, Section D.3): 


P 


2,-2 
ambos Hiblisd co) ] (6.15) 
Pp 
Py ih 
Tf Lo is the valueiof & where .— = 7 
Pp 
Py la eee 
—j = [le 0.410) ] (6.16) 
p 2 


The available data have shown that Equation 6.16 is a reasonable approxi- 
mation. If nN; No are the values of n where the minimum pressure, 
P, and its one-half value occur, the present data indicate that n 3.1 


and Ny = 5.1. Since 2=(M- n)b, Equation 6.16 can also be written 


as: 


Pa 


— = [1 + 0.1025(n - oi ik (6.17) 
P 


It is now possible to derive an Equation for the stagnation pressure, 
Ps? by equating the integral of the wall pressure to the normal component 


of the initial momentum of the jet: 


i pdx, = pu,“ sin 9 (6.18) 


Using the previous expressions, it is possible to show (see Appendix D, 


Section D.4) that the above equation reduces to: 


pa AS : Re ea (6.19) 
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where 8 is the ratio p/P, and b> by are the values of b_ for the 


positive and negative sides of x, respectively. Note that this nota- 


1 
tion was introduced in order to restrict $ in the range of values less 
than 90°, since p is the same for supplementary angles. In this 

s 


notation, the empirical relation expressed by Equation 6.12 can also be 


written as: 


f= Gey" 5 gs 90° (6.20) 


Introducing this into Equation 6.19, the stagnation pressure becomes: 


= (the een NaDPi ie and (6.21) 


: ra (r 44st) Yee 


pu, (2 1 
With the aid of the graphs of Figures 6.6 and 6.7, the above is plotted 


a. |r 


ci 


in Figure 6.5 and is seen to agree with the data reasonably well. This 


serves as a check of the presented results on pressure distribution. 


It remains to determine the extent of the impingement region 


along the wall. Letting x, denote the end of the impingement region, 


x could be defined as that value of x where lp| = 0.02 Pat When 
negative pressures occur - will not be single valued. Hence, an 


additional restriction is imposed on so i.e., it should be the largest 


value of x) satisfying the above condition. The results show that 


x,/H depends upon 9, as indicated in the following table. 


b ae i 


; y s > a) A 
‘ ' an i U Ex 7 7 7 7 
. Het ’ ri | isy. = h Sine ae 
y A , cy i 
th i ‘ oe) = im: 
av { a’ 
uy Mf a 


- 9d oats +H. ns te tr 


: rr, - TT » j ae : 
ney i ve f > . 4 cA ian tc a oa B 


ofy te2 ¢ to embed pre 
-eton atid ade oii xia 2eqeor fa te! wabte -” 
vest woukev lo apaan ai ot » leocamatins usennlic 


pkda at -anhgan 


ts 


(05.8) 


= 


(1g.8) > 


pose eb avads ods. <%.% + cent wines ewes os Yo Bt ofa 
batt few eideacanet capo sds doze song 00 se oat bs 2.9 02 | 
swosaaahinis aaa oe saison hoseonox aay toto as pe 


Ve nec 


| ‘wokgex somapyatgest as to sot ‘aerials tee to 


_ wont snsomaniqad add Yo See i gant JY gabtaoa’ tow ofa gad 7 
er ee ee _ ora ORES 


nis yesaidt Jbauvlev oigate od’ You Chbw £ ace vega: 


ssbguad oa ap, bhaorta Hh 4 sO8 ce fo hesogat ipa 

crioe sot eatupen om Lavtsthees’ oveds tilt gali@alses - ta pee 
re Mats » 

patent ay is eee 


Gaol 


ere ‘ariwotte “os ‘a 


Twat glen re 


3 » aga) 


i] 5 . ' ia 
Ve A Le 4 an 7 @ ste. i a 


5. a a! 7 
; af aa id ie 
Y «= a ee 


ESD. 


TABLE 6.2 EXTENT OF IMPINGEMENT REGION 


The value of x, /H increases for decreasing 4 in the range (o < 130°), 
owing to the "flattening" of the pressure profile and to increasing 
values of b. This trend is reversed in the range @ > 130°, owing to 
the Gecutrence of negative pressures which cause the impingement region 
to os considerably longer than it would be in their absence. It is 
noteworthy that in this latter range ( > 130°) values of x,/H pre- 
dicted from Equation 6.17 are in good agreement with tabulated values 


which were determined directly from the data. 


6.3.2.2 Wall Shear Stress 

Typical wall shear stress variations are shown in Figure 6.10. 
When @ < 90° measurements with a Preston tube present no difficulty. 
However, when @ > 90°, it is not possible to use an ordinary probe 
because the inclined nozzle leaves little room for equipment. A special 
probe was built to accommodate this configuration as shown in the 
following sketch. However , the length of the suspended part of the 


probe becomes enormous and measurement is not as accurate as with an 
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ordinary probe. This consideration is demonstrated clearly in Figure 
6.10. Considering the shear stress ps to be positive for x > U. 


then it will be negative for x, < 0, and naturally ue must vary with 


iE 


x4 in a continuous manner. Therefore, the following condition must be 


satisfied: 


where x) = of means that xy approaches zero from positive values. 
Hence, in the graphs of Figure 6.10 the slopes of [e140 for supple- 
mentary angles of impingement should be identical. Inspection shows 
that this is not the case, the deviations increasing with decreasing 9. 


Thus, shear stress measurements for @$ > 90° must be regarded as approxi- 


mate. 


In order to predict the wall shear stress the same reasoning 
will be applied as that used in Chapter IV. It is firstly assumed that 
the stagnation streamline intersects the wall at right angles. This is 
not unreasonable in view of previous considerations (see Section 6.3.2.1). 
If the equations of motion are written in a curvilinear system (t,> ty) 


where ty is distance along the stagnating streamline and t, is distance 
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perpendicular to it, then some curvature terms will arise. However, 
since the stagnation streamline intersects the wall at right angles, the 
curvature terms will vanish very near the wall and the coordinate system 
(t> t.) will coincide with the system (y, x,) shown in Figure 6.1. 


Using similar reasoning as in Chapter IV, one could write: 


dt 9 
Suds ( = (6422) 
OV fae 

y=o 


(See also Appendix D, Section D.1.2) 
For the case of normal impingement, it was found (26) that pressure 
distributions were geometrically similar with scales Bs and ee where 


Pp. was the centerline value and De was the value of x, where p/P. 


A! 
=D. Ds Assuming this to hold when 6 # 90°, with Po being the value 


of pressure on the stagnating streamline, one could write: 


-kn 
= : =x Di 3 = ¢ > k= 0.693 (6.23) 
Pp P.(n,) 3 eee g p ( 


Using Equations 6.22 and 6.23 it is possible to show that (see Appendix 


De Section D.L. 2): 


t/t, = erf(Vk n) - c,ng(n) (6.24) 
45 

SF = 1G) (6.25) 
pu, fz 


where Taw is the maximum value of o and Cc, is a coefficient which 


depends on }. The functions Cy and h are left to be determined 


empirically. Obviously, Equation 6.24 cannot be expected to hold when 
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negative pressures occur. Figure 6.11 shows the data plotted according 
to the formulation of Equation 6.24. It appears that a single curve can 


describe these data well, choosing C, = 0.2. The variation of C) with 


@ can be obtained directly from the data, using Equation 6.24 and work- 


ing backwards. This would show that Cy does vary somewhat with 6. 


However, the term associated with C in Equation 6.24 represents a 


1 


small correction and for prediction purposes, an average value of Cc, 


is sufficient. 


When negative pressures occur (¢ > 130°) the above analysis 
is not valid. This is demonstrated in Figure 6.12, where data for 
d = 150° are only crudely approximated by Equation 6.24. Considering 


the maximum shear stress, isi dimensional analysis will show that: 


vom H oe 
ae a h(o, RB.) (6.26) 
2 d H 
pu, f2 
H : ‘ 
where Ry = R, yi. However, the effect of viscosity on ee 


appears to be negligible, as was also found in Chapter IV, so that Equa- 
tion 6.26 agrees in essence with theoretical predictions (Equation 6.25). 
This is verified in Figure 6.13 where data, covering a range of Ry 
from about 20,000 to about 270,000, define a single curve. In the range 
® > 90°, lan (re oe error should be expected to be inherent in the 
data of Figure 6.13, according to previous considerations. A correction 
can be applied if it is assumed that the error is proportional to the 
measured value of shear stress. For the case @$ = 120°, this correction 


amounts to about 20% of the measured value as shown in Figure 6.13. It 
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can be shown (Appendix D, Section D.1.2, Equation D.28) that the maxi- 


mum shear stress at the side of negative x is related to the corres- 


ponding value at the side of positive xy by the equation: 


om), b, 
= eae (69827) 
Com Py 


This is shown plotted in Figure 6.13, where it appears to agree with the 


corrected data. 


Finally, owing to the presence of strong pressure gradients, 
the present measurements of wall shear stress are likely to involve 
errors greater than 6% (see Chapter III) when n is less than a limit- 


ing value, Nee This depends on 9 and is tabulated below. 


TABLE 6.3 VALUES OF n 


6.3.2.3 Maximum velocity, un 
Begining at the stagnation point, a boundary layer grows at 


the wall. If the value of longitudinal velocity at the edge of the 


boundary layer is denoted by ui Figure 6.14 shows typical variations 
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of this quantity. It is seen that u increases from zero at the 
stagnation point, to a maximum value and then decreases continuously. 
Regarding the accuracy of the data in the range ¢ > 90°, considera- 
tions similar to these with respect to wall shear are valid. If Wee 


is the maximum vaiue of u. then Figure 6.15 shows that profiles of 


u, are geometrically similar in the sense: 


= £m) (6.28) 
in the range 9 < 130°. 


Beltaos (25) and Beltaos and Rajaratnam (26) have shown that 
in the case of normal impingement, u, can be predicted by application 


of Bernoulli's equation, i.e.: 


pu.” 
= : = a 6.29 
SR = Py 3 b= 90 (6.29) 
This can also be written as 
un = 1. =) me 
Umax Ps 
u . 
2 2 
ne | See (ie apes (6.30) 
Unax 


Equation 6.30 is known to be valid for d = 90°. However, since Equa- 
tion 6.28 holds for all © less than 130°, it follows that £m has 
the form given by Equation 6.30. This is plotted in Figure OvbS wand 6 


seen to describe the data reasonably well. Yet, in the range dy > 130%, 
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this equation is only a rough approximation to the data, as seen in 


Figure 6.15a. From dimensional considerations, it can be shown that: 


Umax H i = 
alee, | Eb) (6.31) 


and this is verified in Figure 6.16. It is noteworthy that for 9 < 90°, 


£(>) is approximately a constant equal to 2.8. Since u is a continu- 


du 
m 


dx, 


ous function of x the value of ( must be the same for 
= 0 


6 
1 
supplementary angles of impingement. Using Equations 6.30 and 6.31, 


ihe 


this condition can be expressed as: 


Cul) 
max (b) 
180-9 - _ 180-9 (6.32) 


) cb) 


This together with Equation 6.12, leads to the following summary for 


u : 
max 


2.85 for =< 90° 


Umax H Bs 
se th 7 (6.33) 


1 
2.8 (Bee) sefou) aq. > 90" 


This is shown plotted in Figure 6.16.. The discrepancy between predicted 
values and data in the range ¢ > 90°, is likely to be due to the inac- 


curacy of measurements in this range. 


6.3.3 Wall Jet Region 


A short distance after the end of the impingement region, the 
flow assumes the familiar pattern of the wall jet. As no measurements 


were performed in this region, this section is mainly an attempt to re- 
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analyze existing data on the basis of the general theoretical considera- 
tions presented in Chapter II. The similarity of velocity profiles 
verified for $ # 90° by Schauer and Eustis (14), leads to the following 


(see Chapter II): 


6» = Cyx 

awe (6.34) 
u = Ca 
m m 


where the overbar denotes that x is measured from a virtual origin 
and not necessarily from the intersection of jet centerline and wall. 
Furthermore, the exponent a is given by: 


C 
Re -$ +54) (6.35) 


2FC., 


where F is a constant representing the quantity: 


co 
hy 2 
Ce aC.) 
u 
™m 2 
fo) 9 
ou 
equal to about 0.7, and C, is the local skin friction factor tol sie 


Consider first the skin friction factor Ce. According to 
previous considerations (Chapter II) this should be independent of x. 
The data of Schauer and Eustis (14) were re-analyzed and are plotted in 
Figure 6.17, in the form u,/u,, Ve" x/H: These data were obtained 
for values of Ro between 39,000 and 43,000. For ¢ varying between 
30° and 130°, and H/d between 20 and 40, it is seen that u,/u and, 


therefore, C is indeed independent of x/H, H/d and 9. Thus, for 


£ 
an average value of Ry = 41,000, Ce = 0.0054. Since C, isa 
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function of only Ro» earlier results, reported in (25) and (26) for 

@ = 90°, as well as data for the classical wall jet (¢ = 0°), should 
apply. The available experimental results are plotted in Figure 6.18. 
For the case of a classical wall jet, Sigalla (32) proposed the follow- 


ing Blasius-type relation: 


u a OJ 


Cc, = 0.0565 (“ 
£ v 


0.25 _ » 0.25 


Using his results on us and 6, it was found that (u6/v) 5 


F(x/d), where F(x/d) is a mild function, varying from 0.67 to 0.75 in 
the range of measurements, i.e., from x/d = 15 to 65. Considering 
that Sigalla's measurement of To involved an error of about 15%, it is 
permissible to assume that F(x/d) is constant with an average value of 


a With this, the above expression becomes: 


1/4 


G = 0.0790 Bs (Sigalla, re-analyzed) (6.360) 


£ 


5 
This relation was obtained for an estimated range of Ro of ~10 to 


mS = Lr, and it is plotted in Figure 6.18, where it is seen to be some- 


what steeper than the general trend of the data points. 


Also for the case of a classical wall jet, Myers et al (33) 
have reported measurements of wall shear stress for values of = in the 
range 7,100 to 56,500. An average relationship for Cy can be derived 


from their results as: 


Cx ‘= 0, 01g oa a (Myers et al, (G37) 


. re-analyzed) 


This is also plotted in Figure 6.18, together with three data points 
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computed directly from their tabulated measurements. The solid line 


drawn in Figure 6.18 is an average of all the data and has the equation: 


Ce = 0.0474 Bane (average) (6.38) 
Consider now the length scale Soe Equations 6.34 state that 6, 
varies linearly with x. The coefficient Co» describing the gradient 


dé, /dx is not expected to depend upon @$ or H/d, if account is taken 
of the derivations leading to Equation 6.34 (see Chapter II, Section 2.2). 
The data of Schauer and Eustis are shown in Figure 6.19. It is seen 

that a single straight line fits the data well, even though 9 was 

Ce Ae 30° to 130° and H/d was varied from 20 to 40. The equa- 


tionloLathis line isi; 


Pages x 
eens 0.089 G +41. 2) (6.39) 
i.e. Cc, = 0.089 and the virtual origin is located at a distance of 
1.2H before the origin of x. For the-case of normal impingement, 


(6 = 90°), Cartwright and Russel (18) reported C, = 0.056 for, Ry in 
the order of 0°. The virtual origin was located at x = -2.5H. For 
the classical wall jet, Sigalla (32) gives Cy = 0.065 for Ro in the 
order of Ties with the virtual origin at x = -0.5d. The results of 
Schwarz and Cosart (31) for $ = 0°, show C, to decrease from 0.085 
at R, = 13,500 £0 0,067 at Ro = 41,600, with the virtual origin being 
between x = - 7.2d and -15.6d. These considerations indicate that 
there exists a certain relationship between C, and Ro» and that 


generally C, decreases for increasing Ro However, it is not possible 
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at present to establish the relation between Cc, and Ro because Cy 
is also dependent on conditions at the nozzle, which vary among differ- 
ent investigations and cannot be accounted for at present. The location 
of the virtual origin appears to depend only upon nozzle conditions, as 


no consistent variation with Ro could be detected. 


An overall average value of C. is 0.0765, whereas the virtual 
origin seems to be located at an average distance of 1.35H before the 


Grapgin of x, t.e.: 


= 0.0765 = 1.35) (average) (6.40) 


For the case > = 0°, the quantity H becomes undefined, however, here 
the virtual origin is located at an average distance of 10d before the 
nozzle. The theoretical prediction for un is given by Equations 6.34 
and 6.35. To obtain an estimate of the exponent, a, consider the term 
C 

on 
2FC., 
this becomes C,/0.107. If this is less than, say, 0.06, it could be 


appearing in Equation 6.35. Taking F = 0.7 and C. = 0.0765 


neglected in comparison to unity and the exponent, a, will be simply 
equal to -0.5. This condition will be satisfied so long as Cy < 9.0064 
which in turn yields the condition Ri 2 20,000. In most practical 
problems, Ro is likely to be larger than this value, hence in what 


~ 


remains it will be assumed that a = -0.5, and: 
eae (6.41) 


u = C 
m 


From dimensional considerations, it is easy to show that: 
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which, in view of Equation 6.41, becomes: 


“sf — &@) 
Ls Vx7H 


ie) 


u C(O) 
or oS (6.42) 
fe) Vx/d 
where Co) is:a function pf >. Thus, the preceding approximation 


with respect to the exponent a has made it possible to show that H 
does not affect un The data of Schauer and Eustis were re-analyzed 
and are shown plotted in Figure 6.20, in the form (uj/u_)? visi x/d. 
Eaetion 6.42 states that when the data are plotted in this form, 
straight lines would result. This is shown to be true in Figure 6.20, 
and the slopes of the resulting straight lines are related to Co For 
the case of normal impingement, Cartwright and Russel (18) gave the 


-0. : ; 
x ae However, when this was re-analyzed it was 


relation u = 
m 
shown that in the range of x/d in which their data were obtained, 


Equation 6.42 was valid with er = 2641. 


Values of Ch are plotted against $ in Figure 6.21 where it 
is seen that the data points define a single curve, a fact that supports 
the conclusion Cc, = CC). This function is now predicted as follows. 
It has been shown earlier that in Equation 6.35, the term C./2EC, 


could be neglected for R. > 20,000. Recalling that Equation 6.35 is a 


TRE 


‘ ‘ The 
rea 
i wos 
1 
: 7 Ea F 
a ae Pees fy aan pis in 
1 7 . ‘ q "i 
teh eu thi ie in i et i 
S u f * : 
rene ry) eas piu e aartas hein 
: ‘ Hi 4 Ny 5 ee bert t Tis wal +) ti ' ey ae yy 
Tn ‘ hae RE ANN Ar een. 2 i a x igs . i pct 
Piney Dea ity D as Vey 6 hop 4 i 


: wea i mn, aw ite boa 


- 


Ms y Ain ? TTA IT eI 
rE nS if Pu 
CS : 
‘ can 2. » oe, pene). 


aol EON oe basaaong ‘vou ax ina ria 
ey trot oat. ce entent * cee 


ies 


Mg a ae 


Ls: 


O09 OvVl OZ OOl 08 09 OV 02 @) 


NOTOGY LAL TIVM SAIVOS ALIOOTIA O¢°9 HuNdTA 


P/X 


Ol 


O02 


(SI1SN3 8 Y3NVHDS 4O Viva) 

O€ O€l © 

O€ Oll B 

O€ 02 ey 
O€ O€ vw 

OV OS V 

0z 06 @ Os 
OV 06 O 

P/H (o)? 


09 


gewrs* avIT 36% BECIO# 
¥XG > 


f 


2trh 


“ent” iN fa, ine ‘ 


af 


\)\ 
7 \ Y 
if 
) ape 
: Pagan 1. bee “pail ie 
é oe ; v, rer 


+ 


H/d 


@ 30 
O co baer 14 


4 20 


GE RER LS 


OP Mag ep 120 150 180 
@ (degrees) 


FIGURE 6.21 THE FUNCTION C9) 


GREE 


a 


180. 


result of the integral momentum equation (Chapter II, Section 2.2) and 
noting that this term represents frictional effects, it is reasonable 
to assume that frictional forces could be neglected from the momentum 
equation. If M, is the momentum of the wall jet in the direction of 


positive x and M, is the momentum in the direction of negative x, 


then: 
M, -M = pu, d cos (6.43) 


co 
2 ; ae ; ; 
But in general M = f ou dy, or taking similarity into account: 
[o) 


epee Aree nr sa (6.44) 
tio 


Letting Cul and Cie correspond to the positive and negative x 


respectively, Equation 6.43 can be re-written as: 


2 2 2 Z 

- = 4 

pdu CaF (Cy C12 ) eu, d cos > (6,45) 
or the equivalent: 
2 Z Z 

r a 6.46 
Cu C2 m cos ¢ ( ) 

where i = =: Note that if o is taken less than 90°, then C4 

Z 


coincides with CD) whereas Cio will coincide with c (180° - >). 
One more relation between Cal and C2 can be obtained if it is 
assumed that the momenta in the two wall jets (which are preserved 
quantities if friction is neglected) are in the same ratio, as in the 
corresponding case.of a potential jet. If it is recognized that this 


ratio of momenta is determined by the flow conditions in the neighbour- 
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hood of the stagnation point, where the flow is locally of stagnation 
type, then this assumption does not appear unreasonable. This assump- 
tion will yield (44): 

ul Lit cos o 


ete ee 1e= .cos o) (6.47) 


Solving Equations 6.46 and 6.47: 


= m/(1 + cos 6)/2 
Oe On 0004; (6.48) 


(@) 
| 


= — 
C2 mY(1 - cos $)/2 
To evaluate m consider the case $= 0°. This represents the classical 


wall jet for which Cul = 3.32 (29), and naturally Cio = 0. It fol- 


lows that m = 3.32. Furthermore, since C2 could also be taken as 


the value of Cr for the supplementary of 9, Equations 6.48 could be 


1 


reduced to a single equation, i.e.: 


CD) /5.5(. + cos 6) (6.49) 


This is shown to describe the data adequately in Figure 6.21. Note 
that the chosen value of 3.32 for m would give a value of Coa 

Ow? ello Od 3. This is not in agreement with the average value of 
0.0765 ehesined earlier (Equation 6.40) and the discrepancy is attributed 
to virtual origin effects. Substituting Equation 6.49 into Equation 6.42, 
the velocity scale is finally given by: 
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ee | | (6.50) 
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6.4 Summary 


The problem of oblique impingement of plane turbulent jets 
has been studied in this chapter. It was established that the free jet 
region extends up to 70% of the inclined impingement height, H, regard- 
less of the angle of impingement. The wall pressure and wall shear 
stress, as well as the maximum velocity near the wall were studied in 
detail in the impingement region. A semi-empirical method, based on 
similarity properties, was developed to predict the wall pressure and 
velocity. The wall shear stress was predicted by the same method, as in 
the previous chapters, properly modified to account for the effects of 
obliqueness. The extent of this region along the wall was established 
from wall pressure considerations and it was found to be a portion of 


H that depends on the angle of impingement. 


The wall jet region was studied utilizing previous experimental 
results. The analysis was carried out in accordance with the theoretical 
considerations of Chapter II. It was shown that the local skin friction 
factor depends only upon the Reynolds number Ro» and an average empiri- 
cal equation was proposed which covers the entire range of experimental 
results (5600 eo <= 300,000). “The length scale was shown to vary lin- 


early with x, however, the coefficient C, describing its gradient 


was found to depend upon the Reynolds number. At present, it is not 
possible to establish this relation because of the presence of signifi- 
cant nozzle effects which vary among different investigators. A satis- 


factory theoretical method was developed to predict the relative velo- 
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city scale u/U5> which was shown to be proportional to (x/d) os 


the coefficient of proportionality oe being a function of only 6. 
This function was determined with the aid of the momentum integral 
equation and an analogy drawn from the corresponding case of a potential 


impinging jet. 
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CHAPTER VII - OBLIQUE IMPINGEMENT OF CIRCULAR 
TURBULENT JETS 


7.1 Introduction 


The oblique impingement of plane turbulent jets was studied 
in the previous chapter. Another configuration that is also of consi- 
derable practical importance, is the oblique impingement of circular 
turbulent jets. At first sight, it may appear that this problem is a 
generalization of the normal impingement of circular jets involving the 
same escalation of difficulty as when passing from the normal to oblique 
impingement of plane jets. This, however, is only a crude representa- 
tion of reality. When dealing with plane jets, the oblique impingement 
involves only one parameter in addition to those describing the normal 
impingement, i.e., the angle of impingement 6. tf a eireular set 
impinges at an angle different than 90° the flow remains axisymmetric 
only in the free-jet region. Therefore, in order to describe the 
impingement and wall-jet regions, it is necessary to introduce an extra 
Space co-ordinate. This could be distance measured perpendicularly to 
the plane of symmetry (i.e., the plane containing the jet centerline 
and normal to the wall), or it could be a polar angle measured from the 
trace of the symmetry plane on the wall. Thus, when passing from normal 
to oblique ere of circular jets, two new variables must be 
introduced, and this makes the problem more difficult than any other 
studied previously. Yet, the understanding gained in previous chapters 


will be helpful in this case. 
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ment region, whereas the wall jet region has been studied earlier by 
Yakovlevskii and Krasheninnikov (27). Their work includes experimental 
and analytical results regarding the velocity field in the wall-jet 

region. The analysis is of the integral type and makes use of a number 
of rather arbitrary assumptions, such as conservation of flow rate and 
kinetic energy in the impingement region. Furthermore, an arbitrary 
physical model is introduced to represent the wall jet. This consists 
of the assumption that the wall jet is produced by a cylindrical nozzle 
of variable thickness. The analysis yields complex equations that are 


difficult to solve. 


The objectives in this chapter are, therefore, two-fold. 
One is to study wall pressure and shear stress in the impingement region, 
and the second to obtain additional experimental data in the wall-jet 
region so as to complement those of (27). Furthermore, a relatively 
simple analytical approach is presented, based on the general theoretical 
considerations of Chapter II, which results in simple analytical expres- 


sions for the velocity field and the wall shear stress. 


Consider now a circular turbulent jet impinging obliquely as 
shown in Figure 7.1. The wall pressure in excess of ambient will be 
denoted simply by Ps? the stagnation pressure is Be and it occurs 
at the stagnation point, S, which is removed from the intersection of 
jet centerline and wall (point 0) by a distance s. For convenience, 
two co-ordinate systems will be used. In the impingement region, it 
was generally found more convenient to use the Cartesian co-ordinate 


system (x, y, z), whereas the wall-jet region was found easier to study 
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(b) Plan view 


FIGURE 7.1 DEFINITION SKETCH 
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by using the cylindrical co-ordinates (r, 6, z). Note that both systems 


originate at the stagnation point S. 


7.2 Experiments 


The experimental set-up has been described in Chapter III. A 
total of nine experiments were performed, whose particulars are listed 
in Table 7.1. The nozzle Reynolds number, Ro» was varied from 35,000 


to 100,000, whereas the relative impingement height, H/d, was between 


Poi beed atlasks 


34,700 
89 ,600 
89,600 
34,700 
34,700 
89,600 
89,600 
100,000 
34,700 
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TABLE 7.1 RANGE OF EXPERIMENTS 


Wall shear stress was measured using a Preston tube (Chapter III) of 
0.046 in. outside diameter. As the line of action of shear stress is 
not generally known before-hand, measurements were performed only in 


instances where its line of action could be accurately estimated. 
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has 3 Analysis and Results 


sore Impingement Region 


The objective in this section is to study the wall pressure 
and shear stress. Measurement of wall pressure is straight-forward 
(Chapter III) even though it is fairly laborious. However, regarding 
the wall shear stress, a major difficulty is present, i.e., the direc- 
tion of shear stress is not known at an arbitrary point on the wall. 
Considering the "limiting" or "wall" streamlines, the shear stress will 
act along the tangent of the wall streamline passing through the given 
point. Intuitively, it appears reasonable to assume that the wall 
streamlines are a family of curves orthogonal to the contours of equal 
wall pressure. Before this can be assumed, however, a better justifi- 
cation must be furnished. For a three-dimensional boundary layer, the 


equations of motion at the wall (z = 0) reduce to (46): 


Cro) 


ne 
dg 
< 
See 
N 
i 
° 
| 
isi 
Ce a 
Q QO 
N a 
i) 
Ne 
N 
I 
e) 


(7.2) 
vs th 
ry * »(&),., 


In Equations 7.1 and 7.2, p is the static pressure above ambient, 
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u, v are the velocity components in the x and y-directions respec- 


tively and u is the viscosity of the fluid. 


If the velocity profile is assumed to be collateral in the 
immediate vicinity of the wall, as is usually the case (46), then for a 
range of very small, but finite, values of z, the ratio v/u remains 
constant on the vertical through any given point (x, y, 0). Quantita- 


tively this is expressed as: 


Peeves a Pe Cesc) acon ONS 2 meee 


with Zo being a small but finite length. This immediately implies 


that: 


2 4 
ov i Aroe av 8 Apo 
(3) y 6(32) ave 7 » | | 
Z=0 Z=O dz /z=0 dz /z=0 


Substituting into Equations 7.1 and 7.2, it is easy to show that: 


(dp/dy) dort es Say 
(dp/dx) |, Tox 


which shows that the total wall shear stress 1 is in the same direc- 
tion as the vector Cees which in ar is orthogonal to the equal- 
pressure contour passing rhough the point (x, y, 0). This then shows 
that the wall streamlines and the wall shear stresses are orthogonal to 
the equal-pressure contours. This conslusion will prove useful in the 


remaining part of this chapter. 
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7.3.1..°Wall Pressure 

Typical wall pressure contours are shown in Figure 7.2, for 
d = 60° and 30°. It is seen that the contour shapes deviate increas- 
ingly from the circular pattern of normal impingement, as @ decreases. 
This, in turn, shows that wall streamlines will no longer be linear 
with the only exception of the streamline on the symmetry plane which 
is the x-axis. An interesting feature of these plots is the absence of 
negative values of pressure at the side of negative x. This is in ap- 
parent contrast with the observations for the plane jet (Chapter VI) 
where negative pressures did occur. Yet, this does not constitute a 
contradiction. It is well known that if an isolated vortex tube forms 
within a mass of fluid, it must end either on flow boundaries (solid or 
free surfaces), or it must close upon itself, i.e., forming a vortex 
ring. In the case of plane jets the occurrence of negative pressures 
was attributed to the formation of a straight-line vortex which ended 
on the sidewalls of the experimental set-up (see Chapter III). In the 
case of the circular jet, a vortex can close only upon itself, thus 
forming a vortex ring surrounding the jet. However, from the side of 
positive x there is no reason why concentrated vorticity should occur 
and, therefore, it is impossible for the ring to form, which explains 


the absence of negative pressure in this case. 


In order to predict the wall pressure, an inductive, semi- 
empirical approach will be employed. Consider firstly the case of 
normal impingement, ¢ = 90°. At a given point (x, y) or (r, 0) the 


wall pressure is P(x, y) = P, (x) because, due to symmetry, dp/d0 = 0. 
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It was found (Chapter IV) that the radial pressure profiles are similar, 


aE e 
P. 
—~ = g(n) 
Ps 

where n= r/b and b is 


7.3 implies that: 


pc, y) = 


The pressure on the x-axis 


P(x» 0) 


It follows that: 


P(x» y) 
P(X Oo) 


profiles along lines parallel to 


value of pressure is used (which 


ate > 


for all values of 96 


figures the following convention 


the value of £5 


is: 


where 


oo |- 42) 


Equation 7.6 states that, in the case of normal impingement, pressure 


is used as length scale. 


as Figures 
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Equation 


(7.4) 


Cs) 


(7.6) 


the y-axis are similar if the maximum 


occurs at y = ©) as pressure scale and 


This property of similarity is retained 
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Hence: 


= g(f) = e (77%) 


P12 


where k = 0.693. Here, however, the length scale b is not constant, 
but it varies with x as shown in Figures 7.4 [(a) to (c)j. Assuming 
large heights of impingement, it can be shown by dimensional analysis 


that: 


x 
a i G >) (773) 
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The data are seen (Figure 7.4) to agree with the formulation of Equa- 
b x 
(=) ef =) 
H 90 b°*H 
straight line parallel to the x-axis. The data indicate that the 


90S) -="const, i.e. ,. 4£ isa 


tion 7.8. When $ = 90°, b 


linearity of b/H is preserved for angles other than 90°, however, the 


slope of the line depends upon 9. Ef b3 is the value of b at 


x = 0 (associated with the profile on the y-axis), the above considera- 


tions can be expressed as: 


(7.9) 


b 
where - and a are both functions of only @$ and, therefore, they 


can easily be evaluated empirically. Regarding the principal pressures 


Py and Po» these are also similarly distributed, much like in the 


case of plane jets, as shown in Figures 7.5 [(a) to (c)]. Defining n 
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as |x| /b b being the value of |x| where P) > “VoIP, the 
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principal pressures are given by: 
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By dimensional analysis it is easy to show that: 
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It is now possible to derive a general expression for the pressure at 


any point (x, y) on the wall as follows: 


p P P 
1.2 
See = = 8(0) 8 (n) 
Ps Pild Ps 
Py ws 2 
or <a UES ORDA pucet te pay (7.12) 
Pp 2 2 
s by 9 b 


In Equation 7.12, by should be used when x > 0 and b, when x < 0. 
On the basis of the above, the wall pressure can be predicted entirely 
if the following quantities are known: The eccentricity s, the 


stagnation pressure Ps» the length scales b,> bo> b3 and the slope 


Q. By dimensional analysis, it is easy to show that: 
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The first two formulations are verified in Figures 7.6 and 7.7 where 
the data for different values of R and H/d are shown to define 
Single curves. Similarly, the expressions involving the length scales 
are verified in Figure 7.8, whereas values of a are plotted against 


g@ in Figure 7.9. 


Qualitatively, the variations of s and P, are similar to 
those obtained for plane jets (Chapter VI). However, here the relative 
eccentricity s/H is ae one-half that of the plane case. The rela- 
tive length scale b,/H increases with decreasing values of 6, like 
in the plane case, but its rate of increase is considerably larger. In 
contrast with the plane case, the relative scale b,/H increases with 
decreasing 9. This is attributed to the inability of the jet to 
develop concentrated vorticity at the side of negative x. Yet, as 


b is always greater than b the pressure profiles on the symmetry 


1 


plane are generally skew, the skewness increasing as @ decreases, 


Diz 


much like in the plane case. This is clearly demonstrated in Figure 7.10, 
where the ratio b/b, is plotted against 9. This ratio is indica- 
tive of the degree of skewness of the pressure distribution and is 

termed the "skewness index". It is now possible to derive some useful 
expressions related to properties of the wall pressure distribution. 

The equal-pressure contour of value p/P. = 1/m can be obtained using 


Equation 7.12: 
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which is symmetric about the x-axis. Dropping the minus sign for con- 


venience, and using Equation 7.9, this becomes: 


b 2 
re 3 nm _ pki 
as Gar - S/H TD (7.24) 


where j = 1 when x>0O and j =2 when x<0._ If ay is the 


area of the contour P,,/P, = 1/m, then an is given by: 


1 
ad fr ydx C7745) 
m 
%9 
where x are the values of x where this contour intersects the 


rip! 


x-axis, i.e.: 


(7.16) 


Using the available results and making some obvious approximations, 
Equation 7.15 can be reduced to (see Appendix D, Section D.5): 
T 2£nm 
2s -—— 757 
a a4) b3 (b, + by) ( ) 
The area of the contour P,/ Ps = 1/2. 4s then: 


25,25 3 b3 (b, +0 >) (7.18) 
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since k = %n2 = .693. Hence, the ratio a_/ao equals &nm/k. When 


Pel? = 1/m, this can be written as: 


(7.19) 


Equation 7.19 states that the aiiineasnre is similarly distributed if 
it is referred to the area a of the contour which a certain value of 
p determines. The scale of pressure is Ps and the scale of area is 
ao» i.e., the value of a when Poe > P.. This could be easily 
checked for the case @ = 90°. It follows that the property of simi- 
larity with respect to contour areas is preserved for all 9. Contour 
areas were obtained directly from the data (with a planimeter) and 
Equation 7.19 is verified in Figure 7.11. This result makes integration 


of wall pressure very simple. The total pressure force acting on the 


wall, Pp. will be: 
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fe) fe) 
anyP 
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Equating ae to the normal component of the initial momentum of the jet: 


tz 4% oe 
Te) Ce ae ee sin > 
Using Equation 7.18, this becomes, after some manipulation: 


aig ay” Be Ee (7.20) 
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where b, b,/H : i, = 12.2393 4 #Whened , = 90°60 bite ba = bs, = b 


and sin 


1. The present data (Chapter IV) gave b,. = 0.078. A 


if 


summary of available data is tabulated below. 


Poreh & Cermak 
Tani & Komatsu 


Hrycak et al 


Bradbury 


Present data 


TABLE 7.2. VALUE OF hye 


2 


Using an average value of b,.= 0.0814, the value of — () at 
90 ou 219 d 
) 


d = 90° is predicted to be 52.2 which is in excellent agreement with the 


average experimental value of 52 obtained in Chapter IV. The empirical 
curves of Figure 7.8 were drawn so that they pass through b = 0.0814 

at o = 90°. Using these curves, Equation 7.20 can be evaluated and is 
plotted in Figure 7.7 where it is seen.to agree well with the experi- 


mental results. 


Considering the eccentricity s of the stagnation point, it 
is not possible to derive an expression from potential flow considera- 
tions, owing to the fact that no potential solutions can be found for 
flows that are not axisymmetric. However, s can be predicted if it is 


assumed to be proportional to a characteristic linear dimension of the 
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wall pressure distribution. Such a dimension is the distance x _ of 
c 


the center of gravity of the wall pressure distribution from S. This 


Cc lo) 
xP, = fh (lates Ci. 21) 
—_ © —- &O 


This equation can be reduced to (see Appendix D, Section D.6): 


is given by: 


Chis aph 2 fe] ib = 2 (7.22) 


Taking a factor of proportionality of 2.9, indicated by the experimental 


results, the eccentricity s becomes: 
Xo 
=wal.J9 i 6i:223) 


This equation is shown plotted in Figure 7.6 and is seen to be an ade- 


quate description of the data. 


With the information already available it is possible, in 
principle, to derive an equation for the family of curves that are ortho- 
gonal to equal-pressure contours. These curves will represent the wall 
streamlines, as was pointed out earlier, and inspection of the available 
relations shows that a derivation of their equations would be rather 
laborious. The end of the impingement region could be taken as the con- 
tour P,,/P. = 1/m with m being a large number so cher P,,/P. mo ee 
For example, m could be taken equal to 50, so that P,,/P. = 0.02. Tt 


is reasonable to assume that beyond the impingement region, the wall 
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streamlines are linear, retaining the direction they had when crossing 
the limiting pressure contour. Inspection of the pressure contours in- 
dicates that these straight lines intersect the x-axis at different 
points,i.e., there is no single focus of streamlines. These considera- 
tions are in agreement with photographic evidence presented in (45) for 
the case of a potential impinging jet. Yakovlevskii and Krasheninnikov 
(27) state that in the wall jet region the streamlines are straight and 
originate at the point O (Figure 7.1). This implies a single focus 
of streamlines. This apparent disagreement with the present conclusion 
can be explained as follows. Even though the streamlines originate at 
different points on the x-axis, the distance of these points from 0 

or S remains finite and is of the order of the distance OS(= ss). 
This distance, however, is about an order of magnitude less than H. 
Measurements in the wall-jet region are made for distances r at least 
of the order of H, hence it makes little difference whether a measure- 
ment of velocity is made by orienting a probe towards S or O or any 


other point in-between. 


7.3.1.2 Wall Shear Stress 

Measurements of wall shear stress in the impingement region 
were performed, as a rule, only along lines where the direction of shear 
stress could be estimated reasonably well. The x-axis constitutes such 
a line and typical shear stress measurements on this locus are shown in 
Figure 7.12. The quantity os is a normalized value of the shear 
stress 1 defined as a e = td, /4pv" where d, is the outside 


diameter of the Preston tube (.046 in.). In Figure 7.12, the data 
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for the case x = 0, y > 0 indicate measurements along the y-axis. For 
the case ¢ = 60°, the y-axis was almost perpendicular to the pressure 
contours, hence such measurement was permissible. Owing to probe inter- 
ference with the flow, the data for x <0 are mostly of qualitative 
value. This is demonstrated clearly in Figure 7.12 where the 
continuity condition (dt,/dx) grat ler (dt /dx) deco is not sat- 
isfied. It may be noted that values of shear stress along the x-axis 
are the most important from the practical point of view, because they 


indicate maximum (x > 0) and minimum (x < 0) variations. 


In order to predict to consider again the equation of motion 
in the z-direction. Strictly speaking, a curvilinear frame of reference 
is distance along the stagnating 


should be used (x z,) where 2z 


he Tg 1 


streamline, xy is distance perpendicular to it on the symmetry plane 
and Yi is distance permendicular to the symmetry plane. However, if 
the stagnating streamline is assumed to intersect the wall at right 
angles, it is permissible to use the Cartesian system (x, y, z) or the 
system (r, 0, z) if z is kept very small. The Reynolds equation of 


motion in the z-direction is (after some manipulation of the equation 


given in (47), and using cylindrical co-ordinates): 


ov, Vv, ov Ov 00 (xt) l dT py 


Vy Dea al ve = ‘ eae " a2 aaah a or ele) a2} 
where Ca ye and Ty, are total (laminar plus turbulent) normal and 
shearing stresses respectively. In particular, 0, is given by: 

_ 5 ov, 
ga Bi pt Pa (7.25) 
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where the overbar denotes time average and the prime denotes fluctuating 


value. At the wall (z = 0) the inertial terms vanish and Equation 7.24 
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Much like in the previous chapters, here again the viscous component of 


becomes: 


ee will be assumed to be negligible in comparison to the other terms. 


Furthermore, it is assumed that the measured pressure p is approxi- 


mately equal to p+ ov! Hence, Equation 7.26 is finally reduced to: 
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where T and T are the components of T inthe: r° iand: 6 
or 08 fe) 
directions respectively. This equation will now be solved for the cases 
@ = 0° and 180°, i.e., along the x-axis. In this case r=x and 
since the flow approaches axisymmetric properties near the symmetry 
or fe) 


plane, terms involving a Wilt Vanish. Denoting’ (7 2). as (Ts 
ole, 0=0,7 


Equation 7.27 becomes: 
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This can be reduced to (see Appendix D): 
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of Nh the value of mn where the maximum shear stress, T m? CCcurs. 
Oo 


These functions are shown graphically in Figure 7.13. Furthermore, n 


m 
is a function of 9 and: 
vom H 2 
Cy, "Sega ) = h() (7530) 
pu, 


The functions a Nn 6? and h(d) are left to be determined empiri- 
cally. In Appendix D, Section D.1.1, Equation D.17 and Equation D.18, 


it is further shown that ae has the same value for both semi-axes 


(x)= 0 <and x >-0) and that if C1). aes 


are the values of T 
om om 


for the semi-axes x > 0 and x <0 respectively, then, in order to 


Satisfy continuity off"1. at x= 0, the-ratio tT (2) p., = must 
Oo om om 
equal the skewness index: 
a 2) by 
aries) = (733) 
i 1 
om 


The variation of ia with $ is evaluated empirically in Figure 7.14. 
Using this graph, and the auxiliary plots of Figure 7.13, Equation 7.29 
was applied to the cases 9 = 60°, 45° and 30°, and the results are 
compared with experimental data in Figure 7.15(a) to (c). It is seen 
that theory is in good Bor eanent with experiments in the main portion of 
the impingement region. It was pointed out earlier that when x < 0, 
significant errors are involved in measurements of ae Litt ts 
assumed that this error is proportional to the magnitude of To? then 
the curves Le sa will be reliable in the range x < 0, as they involve 


only relative values. Regarding the friction coefficient C, defined 
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FIGURE 7.13 AUXILIARY FUNCTIONS Oe qa 
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SIMILARITY OF WALL SHEAR, SYMMETRY PLANE 
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Zaku. 


in Equation 5.30, only the data for x >O0O are reliable. Dimensional 


analysis will indicate that: 


T ie. - 
RS a Ce (7.32) 
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TABLE 7.3 MAXIMUM SHEAR STRESS (x > 0) 


As was the case at normal impingement, there seems to be no consistent 
> 2 : : a 2 
variation of C, with Ro: At the same time, c, & ) appears to be in- 


dependent of the angle 9, as was the tendency in the case of plane jets 


(Chapter VI), and. therefore an average value will be adopted, i.e.: 


c.) = 0.166 (7.33) 


(2) 


Regarding C, , the value of C, when x <0, Equation /.31 provides 


a useful guide. Figure 7.16 shows this plotted (using Figure 7.10) 


o% 5 Re 
[thy 
-iss disne 
i =s ifs 
" a i. : 
= aoe = ®) a a9 tit i Gane 


inobounadl tia | os ‘4 y x vr sau ote uno cob a 
| ‘adit ssnotbat ite 


(S£..%) 


s-thh wd ‘oi orougas 


ih, outs oman ad3 9A | of Asie i. Se esha 7H 
ait seis We gnaw a mt yonabaga sel ew ae 4 sabe ods te ta 


iguk ee wv enh eulev eoert me ernngeens bas’ .(ZV¥ tedqadd) 


ss =" 


f' cs : ae ’ G Py (és) 
CRO i Bn ess rat ian, cy 


fi Te a aT Pat 


| "®, 
to it. ‘ enna .0 > x nodly e to oulav ods, P) getbisge 


” 


» 


cor .¥ o-xiigh pais porsodq, an evo ats ‘ wees obhug, leheon: ; 


4 é 
iar a 
7 7 ‘ 
ee od - 
ye ‘ i - 7 - ub one ip his 
-_ = 7 as. a 5 ‘a aan ie ie (i Pos Lea) faa _! Ln nar 


O MEASURED 
4 CORRECTED 


30 


p 


FIGURE 7.16 


| 60 
(degrees) 


EVALUATION OF C, 


(2) 


90 


Pe 


A = 
7 HK ‘ A ; as i 
I D y ; : 
iM : : 
sss r 7 ! re : h, n 
, 7 Lt ‘ ee 
y 1 1 i ’ 
i he te 
c 
i i 
i = F : i cals 
y om 
a i a a 
f a" oy) = ¥ 3.1 av L 
6 ‘ive =| rer) ' 
: re | ie ‘ 
1 a 
\ 4 j — ig 4 : 
4 i 
~~ 
a a 
‘ ae 
a } ‘ f Y te 


heady (Ose sa 


é ( ¥ 
: , 7 an | mer Be Ee = wee “i 
7 ei ‘i a. i r 5 © egal wetter 
Sei) mye 2 


le ei? . ' 
t ‘ / fe ie 
Fre tL ae 
‘ * ! Ja 
in 
ce | 
i 
* 
i 
j 
” fe 
i 
> I 
© 
¥ + 
j 
1 
4 
j 
m= | 
~ 


, / ) ae , 
FA : { on ‘ i i on 
‘ : : oa ¢ Pie ’ : : in . 
> ity Pee 1.@ a jess tad a Tass, ee | ere 


223 


together with the measured data. It is seen that there is no agreement 
but this does not necessarily disprove the analysis, as the data are 
unreliable in this range. If it is assumed that the error is proportion- 
al to the magnitude of the wall shear stress, it is possible to apply a 
correction based on the requirement of continuous shear stress variation 
at x = 0. The corrected values are also shown in Figure 7.16 and they 


are much closer to the theoretical curve. 


7.3.2 Wall Jet Region 


A short distance after the end of the impingement region, the 
flow assumes a pattern similar in many respects to the type of flow 


known as the wall jet. 


On the basis of considerations stated at the end of Section 
7.3.1.1, it is permissible to assume that for all practical purposes, 
the flow takes place along radial lines originating at the stagnation 
point S. Yakovlevskii and Krasheninnikov (27) assumed the origin to 
be at the intersection of the jet centerline and the wall. Considering 
that measurements in this region are made at radial distances that are 
at least an order of magnitude larger than the eccentricity, s, this 
difference is immaterial. Furthermore, it will be assumed that no sig- 
nificant skewness occurs in the boundary layer. This assumption is 
implicit in (27) even though it is not stated explicitly. It is evident 
that skewness is likely to increase for decreasing values of the angle 
of impingement, so that this assumption can be adopted with the quali- 


fication that it may be unrealistic for the smaller values of }. 
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Yakovlevskii and Krasheninikov (27) reported that profiles of the radial 
velocity component are geometrically similar in the usual manner of wall 
jets. However, their experimentation was limited in the range 45° < 
oO". In order to further explore similarity, measurements were per- 
formed for the rather extreme case @$ = 20°, at polar angles 06 = 0° 
and 45°. If u is the radial component of velocity and u is its 
maximum value in the vertical, plots of u/u vs 2/6, are shown in 
Figure 7.17. The results confirm similarity of profiles at this low 
value of 9. A summary of average data is presented graphically in 
Figure 7.18 together with the classical wall jet profile. Experimental 
results are in good agreement among themselves, even though agreement 
with theory is only fair. It is reasonable to expect that skewness 
would cause the profiles to deviate from similarity, so that it may be 
concluded that at least in the range 20° < > < 90° skewness effects are 


not significant. 


Consider now the following analysis. Neglecting viscous 


stresses, the Reynolds equation of motion in the radial direction is (47): 
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this equation becomes: 
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Introducing the usual boundary-layer approximations, it is assumed that: 
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ae 2 =e and iw 58 ve (7.36) 


Furthermore, since in the wall jet region the radial direction is almost 
the same with the direction of a streamline, it could be further assumed 
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In addition to these, since v << ves it could be assumed that the 


8 


stress V6 ee which acts between radial planes is much smaller than 
the stress ee = ae which is produced mainly by the dominant velocity 
gradient (dv /dz). Considering also the first of relations 7.37, it 


follows that: 
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dv ov. l eae 
‘ror * ‘ete = p 22 ne 
Consider now the time-average equation of continuity (47): 
Ov ov 
197 nets ae 
et Eee a 2?) 
Since v, < v and pe i. this reduces to: 
ine r00 OZ ; 
d(rv_) d (rv) 
oes SUE cy OnE = 0 (7.40) 


These equations (7.38 and 7.40) appear to be the same with those govern- 
ing axisymmetric flows. Yet, the flow is not assumed to be axisymmetric 
because angular variations cannot be neglected, even though they do not 
appear in the governing differential equation. In view of these, the 
flow in this region could be labeled as "quasi~axisymmetric", because 
even though it is not axisymmetric, it does retain some properties of 
the axisymmetric radial wall jet. It is now convenient to employ a 
more familiar notation, i.e., y ekeec designated as u and w re- 


spectively. Then, Equations 7.38 and 7.40 become: 
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These relations are identical with the starting equations of the conven- 
tional axisymmetric wall jet. Assuming similarity of shear stresses 

2 ' aO.. , , 
T= pu, g (6), and proceeding as in Chapter II, it is possible to derive 


the following: 


6, du 6 a6 
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This relation implies that: 
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where the constants C,> Cy» C3 must be universal since the functions 
f, F and g are the same regardless of Ro H/d wenn. | Equations 


7.47 and 7.48 require that: 


6 = C,r (7.49) 


and Equations 7.46 and 7.49 lead to: 
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where a is an exponent and k(0) is a function of O' The integral 


momentum equation, together with Equations 7.49 and 7.50 will yield: 
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a f= = (1 + ) | (7.51) 


where Ce is the local skin-friction factor ge 1s is the wall 


shear stress in the r-direction, and Fy is a constant equal to 
co 


2 
Hee de. The skin-friction factor. C can be shown to depend only 


f 
fe) 
upon Ro? using the same arguments as in Chapter II, i.e.: 


qa* Cc, (R,) (FiR52) 
The quantity u,/u = VC ./2 is plotted against r/d for different 


values of ¢$, 6 and H/d in Figure 7.19. It is seen that u,/u, 
does not depend on r/d. Values of Ce obtained from this graph are 
tabulated below. Inspection shows that the only dependence of Ce is 


indeed on the Reynolds number, Ro: 
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TABLE 7.4 SKIN-FRICTION FACTOR 


Data from Table 7.4 are plotted in Figure 7.20, together with the empiri- 
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cal equation for Ce arrived at in Chapter IV. It appears that the 
exponent —-0.3 of Ro gives a variation somewhat steeper than the general 


trend of the data. Hence, the following relation is proposed, using an 


exponent of -1/5 as in the case of plane jets, i.e.: 


nile 


C. = 0.0980 Ro (ig53) 


f 


Figure 7.20 shows that Equation 7.53 is a reasonable average of the 


available results. 


Regarding the length scale So» it was shown earlier that it 
grows linearly with r, at a rate that should be a universal constant, 
apart from some Reynolds number dependence. Yakovlevskii and Krashenin- 
ikov (27) reported experimental data for © = 45°, 60° and 90°, H/d = 
3.5 and 10 and @ = 0° and 90°, at a value of Reynolds number, R_ = 


9 


69,000. The length scale 6, was found to grow linearly with r at 


2 
a slope of 0.075, which is close to 0.0792 and 0.0765 that were estab- 
lished in Chapters IV and V respectively. The present results for the 


case @ = 20° are shown in Figure 7.21, and they are described adequately 


by a straight line of slope 0.0755. 


It remains to evaluate the velocity scale us Consider 
first the exponent a in Equation 7.50. Strictly speaking, this is 
always less than -1 as Equation 7.51 indicates. Using C, = 0.075 and 


F, = 0.7, the term C./4F,C becomes C./0.21. This will be less than, 


2 
say 0.06, so long as Cr is less than 0.0126, a condition that will be 


satisfied so long as R_ > 20,000. In this range, it could therefore 
fo) 
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be assumed that a ~*~ -l, i.e., Lr The velocity scale, u> is 
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also given by: 


22 
u, = £(p, pu, d. H, $9, r, 8) 


which by dimensional analysis becomes: 


_m . =\ fib. 0, (7.54) 


, 1 , 
Since us a rune it follows that: 


an H _ b@, 0) 
Ud r/H 
fe) 
u 
_m _ h@, 8) 
oF U =F (7.55) 
fo) 
which implies that un is independent of H. Furthermore, if U6 / Un 


is plotted against r/d for fixed $» and 98 , straight lines should 
result. This is shown in Figures 7.22 (a) and (b) where the data define 
straight lines after a certain distance from the stagnation point. Note 
also that data with same $ and 6 but different values of H/d, fall 
along the same line. The virtual origin seems to be located at the 
stagnation point, with the exception of the case > = 20°, 6 = 45°. 

The slopes of these straight lines equal the reciprocal of h. The 
data of (27) were re-analyzed in the same manner and it was again found 
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Experimental values of h($, 6) are tabulated on the following 


page, together with values of r/Hs to being the value of r where 
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us begins to obey Equation 7.55, and it can be thought of as the begin- 


ning of the wall jet region (see also Chapter ry"). 


to a 
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34,700 
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Yakovlevskii & 
Krasheninnikov 


(27) 


TABLE 7.5 EXPERIMENTAL VALUES OF h(>, 0) 


This table shows that h depends only upon 9 and 90 according to 
theory. The quantity r/H can also be shown to depend only upon @¢ 


and 6. This appears to be the trend, however, the data on r, are 
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not as consistent. It remains now to evaluate the function h(d, 6). 
Consider first the particular cases @ = 0° and §@= 180°. # These cases 


are of particular interest because they represent maximum and minimum 


values of h for a given . Let: 
h, = hd, 0°) 
(7.56) 
h, = h(d, 180°) 


A relation between hy and h, can be obtained from momentum considera- 
tions, as follows. Consider an elementary control volume determined by 
two planes intersecting at the centerline of the jet and whose traces on 


the wall form a small angle AQ about the x-axis, as shown in Figure 


In23 


AO/2 
WALL PLANE 


FIGURE 7.23 CONTROL VOLUME 
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The planes are OPM and OPM whereas the plane MP,P. is the plane 


defined by the outlet section of the nozzle. The angles OMQ, OQP, 
and OQP, are all right angles. If the angle is not very small, 


then the elementary momentum in either direction of the x-axis will be: 
lee} 
am = | pu*rAQaz = pAdu 78.rF 
p p Pm Oca 
fo) 


or, using Equations 7.49 and 7.55: 


Pie 5tED 
aM, 5 = 0F,C,A6h, U, d (7.57) 


which shows that dM is preserved. Regarding the term hy 9 this 
> 


should be understood to be hy Whenm exer Ommrancd h, when x < 0. The 


same applies to dM, 2° The difference: 
> 


Z 


D ) (7.58) 


22, 2 
dM, a dM, = pF Cou, d A8 (h, = h 


will equal the x-projection of dM. the momentum at the nozzle contain- 


ed within the angle ag* (see Figure 7.23). Strictly speaking, in 
order to evaluate dM, >? distance r should be measured from the 

. 
stagnation point S. However, when $ is not very small, s/H < 0.1, 


and if r/H is taken sufficiently large, no significant error will be 
introduced. Another objection may be raised because some momenta may be 
entering the control volume near the point 0, owing to the bending of 
streamlines there. These are assumed to be negligible, which is not 
unreasonable if $ is not very small. Finally, frictional forces are 


neglected, and this is reasonable for R, 2 20,000 (see also Chapter IV). 
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The elementary momentum, dM | is: 


1 2 
A®@ 2 1d 
dM = 2 oT pu, he 
hence: 
2 2 2 2 iT Daeae. il 
pF, CU, d A8 (h, - hy ) = GZ emd US A8” cos > 
which reduces to: 
1 
2 2 Aé 
hy - hy = Ky Ge cos 
il 1 
where K, = ———. Ti AD oe AG are very small, A@ = tan AO, 
z 4F,C, ; 
(from Figure 7.23): 
Si: 
ao ‘ alr i S s ae Sei ee 
A8 AG 72 MQ” 
Thus, Equation 7.59 becomes: 
2 eval cos > 
DI oe thas stand 
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(7.59) 


hence, 


(7.60) 


(7.61) 


One more relation between h, and h, is obtained by considering the 
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elementary flow rate, dQ between two planes normal to the wall, located 
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at a small distance Ay/2 on either side of the x-axis (see sketch on 


previous page). 


co co foe) 


dQ ff vovas Ay fv tye,8, f € Age 
0 0 OM. 
Fo 
ise: dQ,» = AyFC)(U,d)hy > 
dQ h 
and a = = (7.62) 
1 1 


Note that dQ onthe strip Ay is preserved, i.e., it is independent 
of xr. Furthermore, the flow in this strip has the properties of the 
symmetry plane and approximates a two-dimensional condition. Hence, it 
is assumed that dQ,/dQ, is equal to the ratio of flow rates from a 


corresponding plane jet of ideal fluid, i.e., (44): 


ce? pee dey’ cues (7263) 
Qs h 1+ cos > : 


The experimental results, taken from Table 7.5 confirm Equation 7.63, as 


shown in Figure 7.24. Solving Equations 7.61 and 7.63, hy and h, are 


given by: 
VK, 
hy = — 1+ cos $ (7.64) 
Y sin > 
VK, 
h fe 1 1- cos $ (7.65) 
2 2 
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d (degrees) 


FIGURE 7.24 VERIFICATION OF EQUATION 7.63 


Data on h, and hy are plotted in Figure 7.25. When @¢ = 90°, 


ie On Behe The average value of h in Chapter IV was found to 


90 


be. 2.0; The data of (27) indicate hog = 1.19. Furthermore, it 


appears that the experimental data of Figure 7.25 will be best described 


at hog is taken as 1.10 which is also an average between 1.0 and 1.19. 
Moreover, the definition of KX = we will give hog = eb Ua 
te 2v4F{C, 


oe. oF is taken as 0.7 and C as 0.077 which is an average between the 


a Z 
value 0.079 found in Chapter IV and 0.075 found in this Chapter. Thus, 


Equations 7.64 and 7.65 become: 


h, = 1.10 Eee poe (7.66) 
Y sin 
he ce. Tea 1 - cos ¢ (7.67) 
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FIGURE 7.25 THE FUNCTIONS hy AND hy 
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These are plotted in Figure 7.25 and are seen to describe the data ade- 


quately at least for 20° <9 S90". Considering the limiting case 


@ = 0°, it is seen that lim h, = 0, which is in agreement with physical 
$0 
considerations. However, lim h, =o, which is physically unaccept- 
+0 


able. This is attributed to the fact that for very small 94, the 


formulation of Equation 7.59 is no longer valid. 


The function h(¢, 8) can now be predicted in general, as 
follows. Consider the locus of points on the plane z= 0, for which 
the maximum velocity in the vertical, uw. is constant. These loci will 
be called “isotachs" as they are lines of equal uy even though they 
are not, strictly speaking, such. Some contours of equal u, 2re 
presented in (27), for $ = 45°, and they appear to be almost elliptical 
in shape. It will be assumed, therefore, that for any angle the isotachs 
are elliptical. For the case @¢ = 90°, the isotachs are circles which 
indicates that if a, b are the long and short semi-axes of the ellipse, 
then the ratio b/a must be unity when @¢ = 90° and it should decrease 
for decreasing values of $. Consider then one such ellipse as shown 


in Figure 7.26. 
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FIGURE 7.26 ASSUMED GEOMETRY OF ISOTACHS 
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This ellipse will be completely defined if the semi-axes a, band the 
eccentricity e€ from the stagnation point are specified. Obviously, 

at. @ = 90°, an="b* and= ce= 0: For the following analysis it is only 
necessary to assume that at least one isotach is an ellipse. For any 
point (r, 8) on the periphery of the ellipse, uf Us is constant, 


ae em: 


ee const = sigs (7.68) 
£ Yr 


2 
Ce ie tae eee) aeae es! (7.69) 


=at = a- € ghee Sellen OFS a: 
Furthermore, T =a ce ry =a an = Agan (E/ad Hide 


Fe 


Recalling Equation 7.63, it follows that: 


3 
en cos 
Since (Equation 7.68) a = => Equation 7.69 can be solved to 
1 ik 


yield (see Appendix D, Section D.7): 


| 2 Pee 2 
cos @ cos 8 + Weos 6 + G) sin > sin 0 
a a - (7.70) 


he 
eee 2 
2 awe 
(1 + cos $)[cos 0 + © sin 6] 
This equation was derived on the assumption that at least one isotach is 
elliptical. However, using this relation and working backwards it can 


be shown that all isotachs will be elliptical with the same value of b/a. 


It is now necessary to establish the variation of b/a with }. 
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Consider the case 6 = 90°. Then: 
Om 90 2hoo 
Ti id r/d SS Cfer) 
fe) VA/A, 
u., 
where A is the area of the isotach ao const (1r*) and A, is the 
Oo 


area of the nozzle Pid wae It is reasonable to expect that this law, 
expressed in terms of isotach areas will be generally valid, i.e., for 
@ # 90°. (A similar property was discovered earlier with respect to 


wall pressures in the impingement region.) 


For elliptical isotachs: 


Using some obvious geometrical relations and Equation 7.68, it is 


possible to show that: 
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On this isotach: oA aides  /d ; using Equation 7.72 it follows that: 
h 
b 2 
u mv nr 
SE P| Se ee Ae (7.13) 
US A/A 


Comparing with Equation 7.71: 
(ch, +h.) v2 = 2h 
1 2 a 90 


or, by Equations 7.66 and 7.67: 


/ i = 
: i 
i> - 
co , 
(1%) | So 


erie oe A bee Gigp sino dose 
. vol alt “wad sp0gee od akdonounees on at piieie ‘hi 
rol get filly ehtainanen od nite saute foniont ee a Ts 
‘eld sagen dakw adi tveo bonavenath sdiew worsens: 3 : ; 

a  ecningor nieamgada ond ¢ 


= ; : i 
re 
st ot at noibteligat ome Rene 


Pe ek ae 
Att leer op 


RX 3 Sih yeahs ine ek 


sar: i : i 


| Aa - 
awa d mae dy Se v si ila | ‘teak “s ¥ 


cc 


(£8. %) 


| ik nation fake aston | 


rr vt years Liat, 


we 


: il 
‘ id if 
z + 2 F PS 

4 ¥ . ; a 

{ 


ve - pie fe ; 
_ + nT 


of 
~—.¥ 


7 aoe 


248. 


w |o 


= sin > C1 TA) 


Substituting Equations 7.74 and 7.66 into Equation 7.70, h is given 


PinalLly as: 
- 10 + 
nto, @) 1.10 , $= £OS 2608 8 (7.75) 
Ysin > cos-8 , Sin 68 
sin > 


Experimental results are available for 60 = 45°, 90°. and 135°. Con- 


Sider eirst. 6.=—90°. Then: 
At, $90) = 110 (ain bees ro) 


When 6 = 45° or 135°, Equation 7.75 becomes: 


Cs i ae P15) Lt as07, Cos @ 
h(d, 135°) ae ‘ ; =) Carn) 
Ysin > 1+ (sin 6) 


These are shown plotted together with available experimental data in 
Figure 7.27. It is seen that agreement is reasonably good with the 
exception of the smaller values of d. The variation of the function 
h(¢, 6) is illustrated more clearly in Figure 7.28 where h is plotted 


against 6 for fixed values of $6. 


Finally, consider the beginning of the wall jet region, O° 
It is easy to show that xo/H is a function of $9 and 96 only. 


Adequate data to define this function exist only for 0 = 0%. If ry 


denotes the end of the impingement region, this quantity can be found 
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from the pressure similarity curve and the values of b/H (Figures 7.5 
and 7.8). It was found that the ratio ro/85 is independent of $6 
when © = 0° and is equal to 1.9. If it is assumed that this property 


holds for other values of 6, then: 


O 
Lr. = £ (0) 


Applying this relation at 6 = 90°: 


r 
ES = £(6) = (2) = £(0) = 1.9 
* $=90° 


which implies that £(6) is actually a constant, equal to 1.9. 


7.4 Summary 


An advanced case of impingement on smooth walls has been 
considered in this chapter. This is the oblique impingement of circu- 
lar turbulent jets. In the impingement and wall-jet regions, the flow 
is three-dimensional, i.e., it is necessary to use three space co- 
ordinates for its description. The wall pressure has been predicted by 
means of an inductive, semi-empirical approach. Wall shear stress in 
the impingement region is generally very difficult to measure and predict 
owing to uncertainty regarding its line of action. However, the direc- 
tion of wall shear is known in a particular case, i.e., along the 


intersection of the symmetry plane with the wall. Measurements showed 
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general agreement with theoretical considerations based on the same 


approach as that used in previous chapters. 


For prediction of flow properties in the wall jet region, a 
theory has been developed, which is based on boundary-layer simplifi- 
cations of the Reynolds equations of motion. It was shown that the 
flow in the wall jet can be considered to be quasi-axisymmetric. This 
means that flow takes place practically along radial lines, with 
shearing stresses between radial planes being negligible. The length 
scale grows linearly on each radial, at a universal slope. The local 
skin friction factor Ce remains constant along each radial. Moreover, 
it does not depend on any other parameter of the problem (H/d, 6, 8) 
but upon the nozzle Reynolds number, Ro: Using all available data, 
it has been possible to derive an average empirical expression for Cy 
which indicates that Ce oc bay tee much like in the corresponding case 
of plane jets (Chapter VI). Thus, the wall shear stress can be calcu- 
lated if the velocity scale, u. is known. It was shown that u., ol /r. 
However, area exist significant variations of un, in the angular 
direction. These were predicted by assuming the contours of equal us 
to be elliptical in shape, based on data presented in (27) and on 
available knowledge regarding normal impingement (Chapters IV and We 


It was again verified that considerations of whether the height of 


impingement is "large" or "small" do not influence the wall jet region. 
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CHAPTER VIII - CIRCULAR TURBULENT JET IN AN 
OPPOSING, INFINITE STREAM 


Sab Lneroduction 


The impingement of turbulent jets upon solid boundaries has 
been studied in the preceding chapters. However, impingement is not 
necessarily associated with solid boundaries. For example, one could 
conceive of two streams impinging against each other at an infinite 
variety of geometrical configurations. In this chapter, the impingement 
of a circular turbulent jet against an infinite uniform stream will be 
considered. This constitutes : simple geometrical configuration, 
however, there is relatively little information in literature with 


respect to this problem. 


Jets in counterflow are used as aerodynamic stabilizers in 
afterburners or turbojet engines. For the hydraulic engineer, this 
problem is of interest with respect to efficient dilution of effluent 
discharges. Furthermore, this problem represents a limiting case of 
the more general problem of jets issuing in streams. Depending upon 
the angle between the original directions of the jet and the stream, 
many different configurations are obtained. For example, jets in co- 
flow and jets in cross-flow are the configurations associated with the 


angles 0° and 90° respectively. Jets in counterflow occur when 


this angle is 180°. 


Basic mean flow characteristics are shown in the definition 


sketch, Figure 8.1. The jet issues from a nozzle of diameter d with 
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STAGNATION 
Ve é 
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Ht Ucel USTAGNATION *— 
POINT 


FIGURE 8.1 DEFINITION SKETCH 


initial velocity Ue against an infinite, uniform, stream of velocity 


(u 


1 <0). A stagnation point forms on the centerline at a distance 
x from the nozzle. This distance is called the "penetration length" 
and could be considered to be a quantity analogous to the height, H, 
associated with impingement on solid boundaries. It was shown in the 
previous chapters that H was a "key" prarmeter for the quantitative 
treatment of the problem. Thus , it is not unreasonable to speculate 
that in the present problen, = will also be an important parameter. 
There is a physical difference, however. Whereas for impingement on 
solid boundaries, H is an imposed (or, independent) parameter, here 
the penetration length is a dependent eae. Treatment of problems 


in terms of dependent quantities is somewhat unorthodox, however, it can 


be perfectly logical under certain circumstances. Other things being 
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constant, x will be a single-valued function of the stream velocity 


u This implies that, mathematically, one could treat either para- 


L’ 
meter as independent. 

Returning to Figure 8.1, if u, v are the time-average com- 
ponents of velocity in the longitudinal (x) and radial (y) directions 
respectively, knowledge of the flow field amounts to knowing the 
function u = u(x, y), since v = v(x, y) can be obtained by integrat- 
ing the equation of continuity. As it is common in jet flows, longi- 
tudinal velocity profiles at different distances from the outlet can be 
reduced to a single similarity curve if proper choice of velocity and 
length scales is made (28). The problem then reduces to predicting the 


similarity functions and the variations of the scales. 


The first contribution to this problem appears to be a brief 
paper by Arendt et al, Rouse (48) in 1956. Experimental results were 
Reesenced in a dimensionless form suggested by dimensional analysis. 
These results were rather approximate because only the magnitude of 
velocity was measured, its direction being estimated by a trial-and- 
error procedure. The correct choice of determining parameters of the 
problem was pointed out, and it allowed derivation of a very simple 
expression for the penetration length. Sui and Ivanov in 1959 (49) 
and Sui in 1961 (50) proposed, after extensive experimentation, empiri- 
cal formulae and nomograms for predicting the velocity and length 
scales and the similarity function. Timma (51) dealt mainly with 
similarity of velocity profiles and the dependence of the penetration 


length upon initial conditions. Ilizarova and Ginevskii (52) reported 
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some experimental results and Ginevskii (53) gave a complex theoretical 
treatment, which however, is mostly of qualitative value. Sekundov (54) 
presented an integral analysis for predicting the penetration length 
and the velocity scale. The resulting expressions are rather complex 
and confirmation was only partial, since comparison with experiment 

was made only with respect to the penetration length. Measurements of 
the mean flow and turbulent correlations were reported by Rao (55). 

The experimental data with respect to turbulent correlations are, 
however, rather crude, as is pointed out by the author. For example, 
the turbulent shear stress was found to have a large negative value on 
the eas which is impossible as it must vanish at this location 
because of symmetry. In the following sections it will be shown that 
the time-average velocity field can be predicted by means of simple 
equations based on theoretical arguments combined with empirically 


established assumptions. 


352 Experiments 


Careful search of the Russian literature (Refs. 49 to 54) 
disclosed that a wealth of experimental data with respect to the mean 
flow was already available. Therefore, only one experiment was con- 
ducted to serve as an additional, independent check on the available 
results. Details of the experimental set-up and techniques have been 
described in Chapter III. The values of Us and uy for,.this,ex= 


periment were 248 and -23 fps respectively. 
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Seo Analysis and Results 


In what follows, it will be assumed that this problem is ana- 
logous to impingement at large heights. By this is meant that the 
penetration length should be sufficiently large so that the main portion 
of the flow in the jet can be assumed to originate from a point source 
of momentum. Arendt et al, Rouse (48) have shown that this assumption 
leads to the result x,/d «xX, if 2 is defined to be the speed ratio, 
u./|u,|. The constant of proportionality was reported [(48),(55)] to 
be about 2.4, which means that if =e is required to be, say, larger 
than 10, 2 will have to be larger than about 4. A more precise 
determination of this limit will be considered later, but these values 
give an idea of the range in which the present study will be valid. 
Assuming then a "large" penetration length, it is possible to distin- 
guish again three flow regions (Figure 8.1). Firstly, in a region 
beginning at the nozzle and extending up to some distance before the 
Stagnation point, one could expect the flow to have considerable simi- 
larities with the free jet. The region close to and upstream of the 
stagnation point is an impingement region where the v-component of 
velocity takes on large values while the static pressure rises above 
the ambient. Finally, the region downstream of the stagnation point 
can be assumed to be occupied by approximately potential flow. This is 
clearly demonstrated in the experiments of Rao (55) where it is shown 
that the turbulent correlations become very small as the stagnation 
point is approached. Transitions between these regions cannot be 


precluded. 
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Sree: Velocity Similarity Function 


A typical variation of u in the radial direction is shown 
schematically in Figure 8.1. From a maximum value of ul on the 
centerline, u decreases continuously, approaching the free-stream 
velocity u; as y increases. If u, is the centerline value of 


Cine Ue Sra then, similarity can be expected to be in the 


1 MW. ec 1. 
form (u - u,)/u, = f£(p)..wonere ol =sy/b..and b fisethat value of 'y 


=u 


where u - uy = $u. This is shown in Figure 8.2 where the present 
data and those of Sui and Ivanov (49) define a single curve. In this 
figure the quantity D/d which appears in the legend denotes the 
relative dimension of the opposing stream with respect to the nozzle 
diameter, d. The conventional free-jet similarity function 

(f = e7 0-693n*) is also plotted in Figure 8.2 and is seen to be ade- 
quate only in the range 0<7n < 1.2. For n> 1.2, the data give 
considerably higher values of f, which in turn shows that in this case 
the jet expands faster than a free jet. Inspection of the legend in- 
dicates that these profiles are in the range x/%,, =< 0.73. No data 

are available near or beyond the stagnation point, as measurements with 
conventional equipment are practically impossible close to the stagna- 
tion point. i order to obtain an understanding of the flow field 
beyond the stagnation point, it is reasonable to assume that it can be 
approximated by the potential flow solution for a three-dimensional 
point source of strength m in a uniform stream, Uy: The source is 


located on the centerline at a distance Xo from the stagnation point. 


This distance will be initially unspecified, but it must be such that 
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the shape of the stagnation stream surface is reproduced as closely as 


possible near the stagnation point. 


x) 


POINT” STAGNATION POINT - 


SOURCE 


ge to ah terete es 


FIGURE 8.3 THREE DIMENSIONAL POINT SOURCE IN A UNIFORM STREAM 


With reference to Figure 8.3, the potential flow solution yields (56): 


(8.1) 
Vv, = Ju, | sin 8 
where k=m/4tT. From these, it is possible to find an expression for 
u - Uy> oe Ae 


u-u. = k cos. o/er (8.2) 
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At y = 0 : u_=u - = k/x,° (8.3) 


u 
m c al} 


where x) is distance along the centerline measured from the point 


source. It follows that: 


u-u x 
er 1 
a ae oa 7 ee 
x y 

If b is the value of y where (u - u,)/u,, = 5: and n= y/b , 

the above can be reduced to: 


3/2 


= [1 170) Sone l= (8.5) 


with b = vY0.59 x (8.6) 


This ascertains that similarity holds even beyond the stagnation point. 
It is now reasonable to assume that similarity exists throughout the 
flow field with a universal similarity function given by Equation 8.5, 


alin, G see 


3/2 


= £(n) = [1+ Gson- le (8.7) 


This is verified in Figure 8.2 where it is seen that this function des- 
cribes the data adequately even when n becomes as large as 4.0. It 

is not hard to explain this agreement in the range = < is Por smait 
values of n, where the flow is of the turbulent jet nature, the 


gaussian law is adequate. For large values of n, the flow becomes 
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almost potential (as can be verified in (55) where turbulent intensities 
tend to zero for increasing values of y), and Equation 8.7 is more 
appropriate. As these two functions practically coincide when n < 1.2, 
Equation 8.7 can be used for all yn. It is of interest to note here that 
a conventional approach, i.e., assuming an eddy viscosity model for the 
shear stress, and then using similarity considerations to reduce the 
equation of motion to an ordinary differential equation, is not fruit- 
ful. The main reason for this appears to be the lack of similarity in 
the distribution of the radial component of velocity v and the turbu- 
lent shear stress T. This can be seen if v is evaluated from the 


equation of continuity and substituted into the equation of motion, 


2 , ' : Ey sails 
using T = pu g(n). Such substitution will indicate some restrictions 
upon the quantities b, uv Up which in turn can be shown to lead to 
absurdum. 


‘This limitation was not realized by Ginevskii (53), who im- 
plicitly assumed similarity of shear stress distributions and obtained 
the shear stress by a polynomial approximation satisfying some obvious 
boundary conditions. The function f(n) describing the velocity profile 
was then derived by means of an eddy viscosity model. This lack of 


similarity is the reason for the somewhat unconventional approach em- 


ployed in the present work. 


8.3.2 Velocity Scale 


For a turbulent jet propagating in an infinite, uniform strean, 


and for a large penetration length, it is reasonable to assume that (see 
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also Rouse (48): 


ee a ay) (8.8) 


2 De ni 
where M, (= pnd U.. /4) ais the initial momentum of the jet. With the 


same assumptions, us will be given by: 
Le Be £4); Pp, u,, x) 


which by means of Equation 8.8 reduces to: 


us as £4» Oyaees =) 
uo, : 
On? DO = LT ee (8.9) 
e) Pp 


everindat ai data are plotted in this form in Figure 8.4, where it is 
seen that, excepting the values i = 2 and 3, they define a single 
curve. Thus, it appears that the present formulation is correct for 

X > 4, which checks out with previous approximate considerations. Con- 
sidering now the quantitative form of Equation 8.4, it is known that 
for a free jet ig Lyx. Also, for jets in co-flow, Tani and Kobashi 


(57) found that ue 1/x. It can easily be verified that here u 


varies also as 1/x for zie ds odie Bink 
Werk 
ciel) Dy Aah Se SON eye Ee ON 70 (8.10) 
U d Dae 


This is shown plotted in Figure 8.4. Note that Equation 8.10 can also 
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be written as: 


ie. * foes. 
U ~  x/d yee, 
fe) 

which is the familiar law of axial velocity decay in free jets. For 


ar > 1.0, Equation 8.3 should apply. By plotting: 


ae 
@ 2) vs. x/x 


Ud 


Oo Pp 


it was found that a straight line, intersecting the axis at o =. 0.6, 


described the data reasonably well. Thus: 


G@) ix 
og: et $  x/x, > 1.0 (8.12) 
Oo (x, /x,) 
where x, =x - 0.8 Soe This relation is verified in Figure 8.4. In 


the range 0.7 < x/x < 1.0, a smooth transition joins the laws of 
Equations 8.10 and 8.12. Comparing Equation 8.12 to Equation 8.3, the 


required strength m of the point source was found to be (Appendix D, 


Section D.8): 


Sion 22) eek (8.13) 


where Q (= nd-U /4) is the initial discharge of the jet. A complete 
fe) fo) 
prediction of un requires knowledge of the penetration length, 5 


which is evaluated in the following section. 


8.3.3 Penetration Length 


Recalling Equation 8.6 and applying dimensional reasoning, a 
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is given by: 
prog d's ae Pog uae 
Pp p 


where of is a constant. An experimental value of os = 2.4 was 
given in (48). Applying Equation 8.9 at es = 1, where ui = 4,> 


it can be shown that ee = £,(1.0). From Figure 8.4, an average value 


of e skcy PAs(\e Hence: 


ir 


x /4 2.6 A . (8.14) 


When A _ becomes less than about 4, the above reasoning will not be 
valid and in this case it could be shown that x,/d = o,@); where 


$, A) is not necessarily a linear function. 


8.3.4 Length Scale 


Using dimensional arguments similar to those leading to 


Equation 8.9, it can be shown that: 
b = £ Gc/x (8.15) 
es 5 ( 7 


This formulation is substantiated in Figure 8.5 where the data define 
a single curve while \X is varied between 7 and 20. The corresponding 


variation of a free jet is also shown for comparison. 


To develop a quantitative expression for the function fo. 
consider the locus u = 0, i.e., y = y); Sekundov (54) arrived at an 
expression for the gradient | dy /dx, which was too complex to permit 


integration and it was assumed, with some indirect justification, that 
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FIGURE 8.5 LENGTH SCALE VARIATION 
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FIGURE 8.6 LINEARITY OF use 
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dy, /dx is a constant, The results available here permit a direct test 


£ thi (pal e = = =_ — 
) : assumption. If ie y,/b then £(n,) u,/u, and No 
~-l A 
pat eo am Further, Io!* = ngb/x,- Using average data from Figures 
m 


8.4 and 8.5, values of eae were calculated and they are shown 
plotted in Figure 8.6, which confirms the linearity of Nee in the 


range sic <— 0.7; Thus, the following empirical relation is adopted: 


y = 0.26 =x (8.16) 


When X + © , the flow becomes that of a free jet. In that case, 
Equation 8.16 will represent the locus of points u/u_ SOO Bee.s 
the locus where u is almost equal to zero. Using the definition of 


y_ and Equations 8.7, 8.10 and 8.16, it is easy to show that: 
fo) 


b Lele 
b= or ef (228)? an 


where € is substituted for 2 This relation ds .verified,.in 
Figure 8.5 at least in the range of available data, Biche sy Oe/3iiyiel OF 
x/%,, >1.0, .the»potential variation of s\b ,will\be ..b, =,)70.59 x1 (Sec. 


8.3.1), which can also be written as: 
b/x, 3 .Ou77 (6.-10.. 8) (8.18) 


In Figure 8.7, the relations expressed by Equations 8.17 and 8.18 are 
shown graphically. It is seen that the two functions have a point of 
tangency at & = 1.4, a property that justifies use of Equation 8.17 


for & < 1.4, and use of Equation 8.18 for € > 1.4. 
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FIGURE 8.7 LENGTH SCALE - GENERALIZATION 


8.3.5 *Geometry 


It is now possible to construct the generalized geometry of 
the flow field. As an example, two loci will be considered in this 
section, i.e., the locus y = y (u = 0) and the stagnation stream- 


surface, y = bp (Figure 8.1). 


The locus Ves cen be found using the definition of Le 
and the present findings regarding similarity of velocity profiles and 
the variation of the scales. Calculated values are shown plotted in. 
Figure 8.8. It is seen that a reaches a maximum value of about 
0.18 ss at zIx, 220777. Ginevskii's (53) theoretical development 


predicts this to be 0.11 see i.e., about 40% less. 


- ae 


ae 


“us ob pinoy it ae asa od Otome, 
siege maton 6b com imige ae ileal 
we ee eo C08 war ahah 


“| 
oa: 


se Sa nbhakath ods gatev bawod- ed: eed oo x pilicoee 
th mabnioad seal ‘fo vatvothmte gorbengor “agatbatt saorg wba 


ae 


ae “fad - zo} noe eri ou Loy bosesuatad caaileos is fective vad 


as oe al . ee 
a “es adeno sepbaereoe aa 2) yntevoni9 EOF e* th * Es , 
Rae Mere) Ud Mores ees 
} | ber Aw, 4 
oP a ab +e. ee a “= i gee +} a od ron preierees a: a 
di \ : Le 
he. ee. eo a eae 
i Nee ae "y | ; at i pl 


210. 


RETURN 
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FIGURE 8.8 FLOW GEOMETRY 


EE ye is the radius of the stagnation streamsurface, then Ys can 


be determined from the definition of this surface, i.e.: 


Vs 


f 2nmyudy = 0O 


0 


From this equation it can easily be shown that (see Appendix D, Section 
D.9) Yolen depends only upon 1 and calculated values are shown 


plotted in Figure 8.8. It is seen that ye reaches a maximum value of 


about) 0.2080 )-"at x/x. = 0.75. 
Pp Pp j 


It should be emphasized that these loci must be regarded as 
time-averaged. Instantaneous configurations are shown in the photo- 
graphs of Figure 8.9. (These photographs were taken using ordinary dye 
injected together with a water jet. The set-up was installed ina 
1.5 ft. wide laboratory flume and it was used for vizualization purposes 


only.) In theory, the stagnation streamsurface envelops the entire 
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flow originating from the jet. However, this is not true in the present 
case owing to turbulent exchange, so that the enveloping curve would be 
considerably higher than the stagnation streamsurface (broken line, 
Figure 8.8). This effect can also be seen in the photographs of 
Figure 8.9, and it is obviously beneficial from the practical point of 


view, where maximum diffusion is usually desirable. 


8.3.6 The Question of Infinitely Wide Streams 


In reality, infinitely wide streams do not exist. Therefore, 
it is necessary to consider the question as to what stream dimension 
can be considered to be infinite. If D is the lateral dimension of 
the opposing stream, and d. is the maximum diameter of the stagnation 
streamsurface, then the ratio D/d, gives an indication of the degree 
of contraction of the opposing stream, due to the developing counterjet. 
Assuming the opposing stream to be circular in cross section, the dis- 
charge Q will be Q= mu, /4. At the section where d. occurs, 
the velocity Ue (assumed to be uniform) will be larger than uy and 


it will be given by: 


4 1 
Se ds Pk SU Be (8.19) 
ey a2 
1- GS? 


If D/d, is greater than, say, 3.3, then us/u, will be always less 
than 1.10 which means that the development of the jet will hardly be 
affected by the size of the stream. An indirect confirmation of this 


limit is presented in Figure 8.10, where the relative penetration length 
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is plotted against \. If the stream can be considered to be infinite, 
then x,/d should obey Equation 8.14. Otherwise the size D will be 
an important parameter. It is seen by inspection that our prediction 
is reasonable, even though it is difficult to establish precisely the 
limiting value of D/d. in this manner. Since d. 42s 0. 29 ie 


this criterion can also be written as: 


D 


which is close to see = 2.0, given by Sekundov (54). The above 


relation can also be written as: 


Be ee i (8.21) 


which implies that the limiting value of D in terms of the jet size 


d is dependent upon the speed ratio i. 


8.4 Summary | 


The problem of a circular counterjet has been treated in this 
chapter from the point of view of impinging jets. It was observed 
firstly that the penetration length could be used as an independent 
variable instead of the stream velocity uy): Thus, dimensional analysis 
was carried out in analogy with previous chapters using the penetration 


length in the role of the impingement height. 


The velocity similarity function was predicted by combining 


potential flow theory with turbulent jet knowledge. The velocity scale 
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was found to vary as in a free jet in the range al < 0.7, whereas 
it conformed to potential flow considerations for — Pe SAN The 
length scale was found to grow considerably faster than in a free jet 
and it was predicted by utilizing the empirical observation that the 
locus of points where u = 0, is a straight line. Finally, the question 
of infinity was considered and it was shown that a stream size D can 


be treated as infinite, if D/x,, is Tareer than about 1.9. 
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CHAPTER IX — SUMMARY 


The flow characteristics of turbulent impinging jets have 
been studied both experimentally and analytically. Starting with the 
relatively simple geometry of a circular jet impinging normally on a 
smooth wall, the study was extended step-by-step so as to include 
relatively complex geometrical configurations such as the oblique im- 


pingement of circular jets. 


With a view to practical application to scouring problems, 
emphasis was placed in predicting wall pressures and shearing stresses 
in the impingement region and a satisfactory semi-empirical method of 
analysis has been proposed. Some preliminary considerations of the 
scouring problem (Appendix A) have shown promise with respect to pre- 


dicting scouring processes on the basis of the present findings. 


Flow characteristics of the wall jet region have also been 
considered. It was shown that a simple analytical approach can predict 
flow properties reasonably well so long as the nozzle Reynolds number 
is larger than about 20,000. onreein was provided by a re-analysis 


of experimental data of previous investigations and by data collected 


by the author. 


Further, the problem of a circular counterjet was studied as 


being essentially an impingement problem and a simple semi-empirical 


method has been proposed to predict the time-average velocity field. 
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Finally, in the course of this study it was shown that when 
the impinging jet is not fully developed when approaching the stagna- 
tion point, the physics of the problem changes drastically in the 
impingement region. A limited study in this direction (axisymmetric 
case) has been performed in order to furnish some fundamental under- 
standing of the difference from the fully developed impingement. This 
could provide the foundation for further work on more complex geometries, 
which was not pursued here, owing to rather limited practical interest 


of this type of flow. 


The present work was confined to study of time-average 
quantities. However, the importance of turbulent fluctuations cannot 
be underestimated. In view of the current trend of attacking sediment 
transport problems in terms of probabilistic arguments, it is evident 
Phar tecbulence measurements in the impingement region are desirable 


and expecially with respect to wall pressures. 


Experiments with rough walls are also desirable, since mobile 
beds are usually rough. However, the difficulty may arise of how 
would wall pressure and shear stress be defined in the case of rough 
boundaries. This is commonly overcome by averaging over an area that 
includes a large number of roughness elements. Is it then correct to 


express forces on roughness elements in the manner introduced in 
Appendix A? Research in this direction (see also Chao & Sandborn (58)), 
would be an excellent contribution to the field of sediment transporta- 


tion. The impinging jet provides a good means of doing so, because of 


the simultaneous existence of shear stresses and pressure gradients, 
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thus affording an opportunity to relate forces on roughness elements to 
the local (areal average) values of pressure and shear stress. It is 
known at present that local scour is caused by two kinds of singulari- 
ties. One is the impinging jet, the other is the vortex located near 
a boundary. It appears reasonable to speculate that an approach 
similar to the present could be fruitful with respect to the vortex 
starting, say, with a simple rectilinear vortex located near a smooth 


boundary. 
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APPENDIX A - A CRITERION FOR SCOUR 


Ae eine roduct ton 


In this section it is shown how the knowledge gained so far 
can be applied to the scour problem. When a turbulent jet impinges on 
an erodible boundary, scour may or may not occur. The first task, 
associated with scour from jets is, therefore, to determine the condi- 
tions under which scour will take place. This question can be put into 
several forms. For example, given the geometry and velocity of the jet 
and given the properties of the erodible bed material, determine 
whether scour will or will not take place. Or, given the geometry and 
velocity of the jet, determine the properties of the bed material so 
that scour will not occur. If scour is expected to take place, the 
next task is to determine its magnitude and geometry and their develop- 


ment with time. 


In this Appendix, which serves as an illustration of practical 
application of the present results, only the first question will be 
considered for the relatively simple case of circular turbulent jets 
impinging normally on a boundary consisting of cohesionless material of 


uniform size D. The answer to this question will form a criterion 


for scour. 


A.2 Theoretical Considerations 


A criterion for scour has been proposed earlier by Poreh and 
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momentum of the jet and H is the initial impingement height and D 
the size of the uniform bed material, then if everything else is kept 
constant and M is allowed to vary, there will be a value of M) 

below which no scour will occur. Poreh and Hefez wrote the following 


dimensional relation for the critical value of M,: 


My rs £(O,Y.> Ds Hy. V) 


where po, V are the density and kinematic viscosity of the fluid res- 
pectively, and xe is the submerged specific weight of the bed material, 
Tee ya og(s - 1) if s is the specific gravity of the bed 


material with respect to the particular fluid. 


By dimensional analysis, the above relation can be reduced to: 


wy VET 
M Vv sane 3 Z 
(@) 

nd2 2 


Since M =0 =a Un » this can also be written as: 
Oo 


3 
H e . ay For ames Dy (A.2) 
d°U 
(6 


Experiments with air jets impinging on sand beds, showed that the in- 


fluence of the Reynolds number, Ro? was negligible* and for the range of 


ee ee Wi el he OER ee ee 

* This is not surprising in view of our previous findings on smooth 
walls. It was established earlier that even for smooth walls, shear 
stress and pressure did not depend on R). Since erodible wails are 
rough, it may be expected that the scouring forces will not depend 


on Ro for a stronger reason. 
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experimentation (D/H = .00029 to .0022), the following empirical rela- 
tion was found adequate: 

3 

H mS ie 

EMAC gente 12.8 = (A. 3) 


de 
(@) 


which can be further reduced to: 


H/d 


aa Se a die 05. 5G (A.4) 
U,/vgD(s - 1) 


Tsuchiya (8) considered the same problem. His approach, however, was 
based on questionable theoretical derivations so that it was not pos- 


sible to obtain any consistent criterion, such as the above. 


Consider now the following simplified theoretical model, which 
makes use of the critical shear stress concept. The total tractive 
force on a bed particle will consist of two components, one due to the 


shear stress oe and another due to the pressure evadient, ap/diug aves : 


aa ot ta Ap ja» 
where T is the tractive force and ay is a horizontal area represent- 
ing the "area of influence" for the particular grain. Obviously, 
a, = cD” with the constant cy being a coefficient which depends 
upon the degree of packing of the bed material. The quantity Ap, 


represents the pressure difference between the front and back sides of 


the particles, whereas ay is a vertical area representing the effec- 


tive area of projection of the grain. Much like ais this area will be 


c p? 1 ae is an equivalent, or total, shear stress, which 


» Vs 


(fiwA) ; is nat ~~ a $.tf 
f sy 
(dA) RB 


aw Favenegd | saline ait now Laois er 
~3iiy 388 any ae ee 08 acco co 


: ei a, ne .oktost 


dodgy. sdstiont shits bottesqute : | 
avbtaee fedon o6T - .agnoeed, sense a Lasttin 20, 
ats od “whi one _apmeormo2 ss He sats ato be «a 


Sud ably soars  isiemab ‘etd of pub roridane ba oF emebli: a 


i 


Se 


~snivesa ion tarda i "ig sa ; | 
‘thveb sions wolustoria ani stall Fanaa to iat edd 2 


* bria dosti gnpt2t¥ta09 & gated ? Yanan ota, ait haa 


Pe weazinnup oi Apiieres had arts to stom, te, sana co 
40 ‘obvi dos bat sage odd onmuaip’ ostis xa 33h orecsnt eis esneea3qa7 . 


~soR as gat zanna nore ddctsaeyn 6b . gf vaoeacdy ceotabsana oils 


gel - “Ri Rew hal ag skit out: Sastong st to mphsantons o_anan avtd 
a Aq.s = 


dats renee: weadu ef995. +0 ‘sanpaetons 4 as at ot 
- ) Si } > i | ee 
A -~ , : ; Pal : ree 


~~ 
se 7 Py = zi - a 
é , 
a ~, ‘, . 4 =m 


288. 


incorporates pressure gradient effects, it will be equal to T/a,, 5 B= vad 


29 
1 a ld Ap. 
el 
dp. 
or x = a - aD ce (A.5) 


where Q = const = cy/c This derivation is ordinarily valid for 


L 
uniform bed materials with grains of spherical shape. Since the con- 
stant © depends mainly upon the shape of the individual grain and the 
degree of packing of the bed material, allowances can be made to extend 
this equation to non-uniform materials of shapes other than spherical. 
Furthermore, it is obvious that this analysis will be realistic so long 
as the grains are very small in comparison to the width of the impinge- 
ment-region, lfm r, denotes the radial extent of this region, it is 
reasonable to restrict particle sizes so that D is less than, say, 
ae@ss 


one-tenth of Pees 


Dice) Oe di xe (A.6) 


For a smooth wall and large impingement heights, r, ~ 0.2 H (Chapter 
IV). When the wall is rough this value is likely to change. Since 
roughness tends to intensify turbulence, it is reasonable to expect that 
r, will decrease with increasing D. When D becomes equal to or 
then the pressure difference Ap, will be simply 


greater than ry, 


equal to De: Obviously, most practical situations will be such that 


condition A.6 is satisfied. 


If t is the critical shear stress associated with the par- 
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ticular combination of fluid and sediment, then a condition for non- 


scouring will be: 
max  T as (A. 7) 


Considering large impingement heights, it was shown in Chapter IV that: 


ay ee a(x). 3 Beet Mes OO shi Ss 9 
with: 

4B 7 
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hes 

9) 2 
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2 d Pp 
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When the wall is rough, the functions g and f as well as the coef- 


ficients C, and Ss will generally depend upon the relative roughness 


D/H,. so that: 


finweG. = 0.16 5 IS i Ngus ly & ee: lim g(\) = 8, 0A) 
D/H>0 p/wo ?P D/H>0 


and lim of (A) = falA) 
ry 
D/H>0 


where 85? fy represent the forms occuring for smooth walls. 


Using these relations, it is possible to reduce Equation A.5 


to the form: 


J (A.8) 
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The quantity inside the brackets is a function of both \ and D/H. 
Denoting this by F(A, D/H), the above equation becomes: 


2 
2 
tT. = LFA, 2 (A.9) 


The maximum of aN will occur at i= w while da is determined 


from: 
oF 
2 ee. (A.10) 
A=) 
m 
which implies that i _ =x ys It follows that: 
m m “H 
Pt hag Si 
Ue pu, GP GG) (Avi) 


where G is a new function which depends only upon D/H. 


Considering now the critical shear stress, ee the well-known 
Shields’ diagram gives the normalized critical shear, T /pg(s - 1)D, 
as a function of the particle Reynolds number, u,D/v. However, the 
influence of viscosity is of little consequence in most practical sit- 


uations, i.e.: 
se = Bog(s - 1)D (A.12) 


where § is a constant. This condlusion is also supported by the 


results of Poreh and Hefez (6), where no viscous effects could be de- 


tected as was pointed out earlier. 


Substituting Equations A.11 and A.12 into Equation A./ the 
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condition for non-scouring can be expressed as: 


D 
G(=) 
H/d H 
> (A.13) 
U6! wy Oi) Ud © 
where w, has been substituted for the quantity V/g(s - 1)D , which 


has the dimensions of velocity, and is related to the fall velocity of 
a single particle in still fluid. Taking the limit when D/H > 0, 
then Le will not contain the pressure gradient term and at the same 


tame ~C, = 0.163 Hence, at the limit: 


H/d Me 0.16 


of eae 
Us! Wy. = B ‘ Gi o at 


This relation will hold for the smaller values of D/H, and it is 
obviously independent of this parameter. This explains the empirical 
finding of Poreh and Hefez (Equation A.4), whose experiments were for 


D/H always less than 0.002. 


Considering now the larger values of D/H, i.e., D/H > r,/H, 
obviously this will represent impingement on a relatively smooth, but 
curved wall. To evaluate forces on a particle it will now be necessary 
to integrate a and p on the surface of the grain. These forces 
will now be praptredinal 2 lls Fy and Ps which in turn are proportional 


to pU ga TRC. It follows that in this case the criterion of non-scour 
fe) 


will simply be: 


ana >paconst wey Dwar bets (A.15) 
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Since r, is likely to decrease with D/H, the condition D > r, will 
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be satisfied at a value of D/H which will be less than 0.2, i.e., of 
the order of 0.1. The criterion of Equation A.13 can also be obtained 
by dimensional analysis, as in Equation A.2, neglecting viscous effects. 
Analogous results can be obtained when the height of impingement is 
small. Using the results of Goanitter V, it is possible to show that in 


this case: 


pu, mi D dh 
sa = als Tag leon (r) = (gh d -t | (A.16) 
dr 
_— us — — Tom 
where h_ (r) hc aheeae: ha = P,,/P, Pui wid des Con eset 05007 
om pu, [2 


Thus, the criterion for non-scouring will be: 


D 
Wy G 1 @ 
U. os y 8 CAP 17) 


whereas when D/d > 0, this will reduce to: 


* fm D 


> const (D> PS) (A.19) 


A.3 Comparison with Experiments 


Experimental results are available from three different sources. 


(i) Poreh and Hefez (6). Data were obtained with air jets impinging on 
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sand beds in the range 0.00029 < D/H < 0.0022 (0.23 mm <p < 0.90 mm): 
(ii) Tsuchiya (8). Data were obtained with water jets impinging on 
sand beds in the following ranges of parameters: 


0.4435 m < D < 15.0 mm 
3,840 < Ry < 58,600 
LO) sop ds) 80 

0500222” < D/H <. 0.0375 
Critical velocities were obtained by keeping other parameters constant, 
while varying the velocity US and recording the relative amount of 
sediment that moved. The critical velocity was defined to be that value 
of US at which 0.5% of sediment in the impingement area of the bound- 
ary was dislodged. 
(iii) Chiu (59). Data were obtained for air jets impinging on a nearly 


uniform sand bed of median size D = 0.26 mn. Other parameters were 


in the following ranges: 


eee, §< 4B/de ah 50.28 
8,780 < R < 42,800 
0.000804 < D/H < 0.00243 


In these experiments, the bed was allowed to be scoured and the final 
maximum (or equilibrium) scour depth was measured. This was plotted 
against Uy in each case (keeping other parameters constant) and the 
value of US Garr poet ing to zero scour depth was obtained by extra- 
polation. This extrapolation was considered to be fairly reliable 


because the scour depth varied with qe. linearly (excepting the very 


high values of U,): 


Experimental data from the above sources have been re-analyzed 


the form suggested by Equation A.i3 in Figure A.1. 
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The data exhibit considerable scatter which, however, is not uncommon 
for scour experiments. Thus, our formulation is reasonably well veri- 
fied and for practical purposes the average curve drawn through the 


data could be used to obtain limiting values of the parameter: 


H/d 
U,/¥gd(s - 1) 


For the lower values of D/H (D/H < 0.0012) this parameter 
remains constant, as was predicted earlier, with an average value of 
3.5. Substituting this into Equation A.14, the corresponding value of 
8 will be 0.013. A conventional value from Shields' diagram would be 
about 0.04. There are two reasons for the discrepancy. Firstly, 
Shields' value was obtained for plane flow, whereas the present data 
are for axisymmetric flow. This difference, however, is not expected to 
account for the present variation of 300%. The following second reason 
appears to be more satisfactory and lies in the definition of the 
critical condition. Shields' diagram was obtained by plotting shear 
stress versus rate of transport and extrapolating to zero rate of 
transport. Later experimental and theoretical studies (see for example 
Paintal (60) and (61) and Yalin (62)), showed that some transport takes 
place even when on) 0. This situation can be explained if the tur- 
bulent nature of the flow is considered. Turbulence is responsible for 
the existence of a certain, however small, probability of movement, 
which in turn, causes some transport to occur over long periods of 


time. Using theoretical derivations of Paintal (61) for plane flows, 


the corresponding probabilities of movement for 8 = 0.013 and 0.04 
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respectively are 1073 and tne 


Another interesting feature in Figure A.1 is the tendency of 
the data to become independent of D/H as D/H increases, in agree- 
ment with previous considerations. Finally, it is noteworthy that no 
significant differences exist between air-sand and water-sand data, 
even though these represent an over thousand-fold variation in submer- 
ged specific weight of sediment. In the phenomenon of sediment trans- 
portation there exists a quantitative difference between transport by 
air and transport by water. This is due to the fact that in the former 
case, (s - 1) ~ 1,000 whereas in the latter case (s -1)~ 1.0. Thus, 
fall velocities are much greater in air. This implies that saltating 
particles, strike the bed with much greater momentum in air than in 
water, a fact that contributes to increase transport. When dealing with 
tiation of motion, however, very little transport occurs and move- 
ment is mostly by rolling, so that the above factor is absent. This 


explains why no quantitative difference exists between air-sand and 


water-sand data in Figure A.l. 


Using 8 = 0.013, the criterion for scour for small impinge- 


ment heights will be: 
U 
) 0.013 (A. 20) 


The quantity Ce varies slightly with nozzle Reynolds number, and for 
m : 


R_ between 31,000 and 54,000 (Chapter V) has an average value of 0.007, 
fe) 
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—— water (A. 21) 


In practice, almost always Us is much larger than w, So that scour 


will almost always occur for small impingement heights. 


A.4 Summary 


A simple application to the scour problem associated with 
impinging jets has been praseaned in this Appendix. It was shown that 
the question of whether or not scour occurs, can be answered by exam- 
ining a single dimensionless parameter which expresses the erosive 
capacity of the jet. Experiments confirmed theoretical predictions 
over a wide range of the specific weight of bed meter ieh: with respect 


to the surrounding fluid. 
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APPENDIX B - ERROR ANALYSIS OF PITOT CYLINDER 
TECHNIQUE 


The Non-Rotating Pitot Cylinder Technique is based on the 
assumption that the pressure readings taken from the three orifices of 
the probe (see Figure 3.7) are related to the local static and dynamic 
pressures p) and Py, respectively by the equations: 


2 lle wl al + K;Py S) REE 235 (Bil) 


where K, are three calibration factors which depend on the angle $¢ 
(see Figure 3.7). Calibration of the instrument was carried out in a 
uniform stream, i.e., Py and De have the same value in all three 
locations 1,2,3 of the orifices. However, if there exist gradients 
of PS and Py in the direction normal to the axis of the probe, then 


it is reasonable to assume that: 
Dp = pl, TK,P a 3a te eee (B. 2) 


where Pod and Pi are the values corresponding to the locations of 


the orifices. Assuming constant gradients between 1 and 3, the follow- 


ing relations apply: 


dp, dp, 
P63 a 5? +A dx A Geen Go ae A dx 
(B.3) 
dp dp 
4 na 4 Vv 
Bog =. Doo tS aoe a eee: 
dp, dp 


where the gradients FO Fm are evaluated at the central orifice, 2, 
x x 


x denotes distance normal to the probe axis, and A is the distance 
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between successive orifices. 


The first step after measurement of Py> Po> P3 is to compute 


the factor: 
pau, 6B 
K = ie aes : d=< 0 
P3truPs 
(B.4) 
Di iP 
Ors K = ee : ¢ > 0 
i Ware a) 


from which the value of can be obtained from the appropriate graph. 


Using Equations B.2 and B.3, K will be given by: 


dp dp 
A ) Vv 
Saat aa ahs E POP ax | 


i 05 OE ee (B.5) 
nse ws Gok a Bax 
v2 
dp. dp. 
Jee aoe and dn «are both zero, then K = (K, - K)/K, = Ky) > and 


this can be assumed to be the "true" value of K at the location of 
orifice 2 where the measurement applies. The extra terms in Equation 
B.5 will be responsible for the introduction of a certain error. As 
"true" distributions of Py» Py are not known before hand, the measured 
values (which are obtained on the assumption that gradient terms are 
negligible) will be used to obtain an estimate of the errors introduced 
because of gradient terms. In order to apply a correction to K, it is 
necessary to add terms equal in magnitude and opposite in sign to the 


gradient terms on both the nominator and the denominator of the fraction. 


Hence, a corrected value of K_ will be: 
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dp dp 
(@) V 
Put Po rar fe eettlede 
KS (B.6) 
dp dp 
= Ee ee 
Ps Po dx 3 dx 


The values of Pi» Po» Ps K,» Ky appearing in Equation B.6 should, 
strictly speaking, be "true" values. These, however, are not known 
before hand and, therefore, as a first approximation the measured values 
can be used. For the case of impingement at small heights, where 
gradients are most intense, error eeedeaton will be provided for two 
profiles. One is in the free-jet region where dp ,/dx = 0 and the 
other in the impingement region, where both dp /dx and dp /dx are 
different from zero. Letting h denote corresponding manometer height 
and considering that here A = 0.0141 in, Equation B.6 becomes: 


| dh we] 
h, - h, + 0.0141 Gee K) + ea 


Re (B.7) 
Vv ie} 
h, - hy - 0.0141 & Se: = 


The analysis of error is performed in Tables B.1 and B.2, where r de- 


notes radial distance and the subscript c denotes corrected values. 


It is seen by inspection that these errors are generally with- 
in acceptable limits, with the exception of large values of distance 
from the center, where the percent error is magnified owing to the very 
small magnitudes of hy and hos Thus, excepting regions where hy 
and ho become very small, the error in static pressure is less than 
3.5% whereas it is less than about 4.0% for the dynamic pressure. Since 


v « Yp , the error in velocity measurement will be less than 2%. 
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The error of the angle ¢ is generally less than about 1°, becoming as 


large as 2° in regions where Po? Pe tend to zero. 
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APPENDIX C - ANNULAR SHEAR LAYER 


Consider an annular turbulent shear layer, originating at 


x = 0, as shown in Figure C.1. 


No 


FIGURE C.1 DEFINITION SKETCH 


The Reynolds equation of motion in the x-direction is (neglecting laminar 


stresses in comparison to T = -pu'v'): 
du Ot (46 od dpe SLOT) oa 
on 0 ox e or or (C.1) 


and the longitudinal component of velocity, u is given by: 
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u = (C.2) 


where n = (r - r,)/b, b.= 72 2 £ being the value of r where 


2 Le 2 


u= 0, and the function f(mn) is closely approximated by (28): 


2 
3/2, (c.3) 


£(n) (1 -n 


The conventional approach in these problems is to pursue a similarity 
analysis (29). However, in this case it can easily be shown that con- 
ventional similarity assumptions [such as v/U, = £,(n) : t/pU,” = g(n)] 
will lead to absurdum. Therefore, an integral analysis is preferable. 
Multiplying Equation C.l by r and integrating from r=0 to r&= To> 


and neglecting the pressure gradient term: 


2 
2 a Ria ees 
if riiidnre €, 2 pide, (C.4) 
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Using Equation C.3, this becomes: 


Z2 
B=” 0.065 0 “baa c13 Wee (C.6) 
fe) oF gi. 
Thus Equation C.4 can be written as: 
—2 —2 ee ce ie 
r, + 0.13 b+ 0.626br,-7F = 0 (Ci7) 


where the overbar is used to denote non-dimensionalized quantities with 
respect to d_ (i.e., Bur rj/d etc.). It is then possible to express 


T in terms of b by solving the quadratic, i.e.: 


Tr, = 0.25 - 0.0325 b* - 0.313 b (C.8) 
This is plotted in Figure C.2 where it is seen that it can very accurate- 
ly be approximated by the linear relation: 


cn x -050.= 0.338 & (c.9) 


One more relation between and b can be obtained if the entrain- 


ment data of Hill (41) are used. Let Q _ be the total flow-rate through 


the jet. Then, Hill's entrainment coefficient C, is defined as: 


aiond : : 
C, = a 4 (incompressible flow) (C310) 
x 
oO 
r 
2 
But Q = OO ae , which after manipulation, can be written as: 
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30%. 


2 2 
§ n{U5y, + 2K,Ub> + 2KU,br |} (c.11) 
itt ik 
where Ky fone = 0.13 and K, Scbetdnv= t0.045 
Oo O 


It follows that Equation C.10 can also be written as: 


Ea. =o ee 1 
25 + 2K)b + 2Kyb r} = Z &9 (C.12) 
dx 
2 
; _ md 
(Since Q, ay ee Us): 


The data of Hill (41) indicate that for small values of x = x/d, {€ 


can be approximated by a straight line, i.e.: 


GC 2a+tBe 3 “= s< 5-6 (C21) 


with @° = 0.083 : 6B *=— 0,036; 
Substituting Equation C.13 into Equation C.12 and integrating, we obtain 


after some manipulation: 


47) +1040 +3605, -axt58x +1 (l14) 


Using Equation C.9, this can also be written as: 
be 4.1.57 Bb - [0.297 x + 0.064 x] =" "0 (C.15) 


which yields: 


po ‘ [o(x) - 1.57] (C.16) 
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ith OG) = 92.46 + 1.19 x + 0.256 x (c.17) 


The function (x) is shown graphically in Figure C.3, where it is seen 


that it can very accurately be approximated by the straight line: 


» | 


d(x) = -1.57 + 0.45 (C<18) 


Using this relation and Equations C.9 and C.16 the scales ry and b 


are given by: 


b Ys 0s er (c.19) 
By) Ph 0-50 0.079 x (C. 20) 
; Hee : 
Lt bi 72 is the value’ of r = at where u= 5 es then Equation C.3 


will show that b = 0.442 b which implies that: 


4/2 


b = Ost0) = (C521) 


Considering now the entrainment velocity, Vo? this can be predicted as 


follows: 
aoe 
2TL 5V, oie 
Vv C, 
Hence: v= il ae which by means of previous 
is o > f 
2 
findings can be written as: 
eee 0; 083 0, 0362x (c.22) 


= g(0.71 + 0.12 x) 
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This is shown plotted in Figure C.4, where it is seen that ve increases 
with x, approaching a constant value of about OL025). Thus, (for = = 


& 4:9, 1.e, for the transition and for the fully developed jet, one could 


assume vs = const = 025, 
0.03 
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0.01 
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FIGURE C.4 ENTRAINMENT VELOCITY 
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APPENDIX D - DERIVATIONS 


D.1 Prediction of Wall Shear Stress in Impingement Region 


Wall shear stress in the impingement region is predicted in 
general by applying at the wall the equation of motion in the direc- 
tion normal to the wall. Then the local pressure gradient normal to 
the wall is related to the shear gradient along the wall. The perti- 


nent equations are different for plane and circular jets. 


Deel Circular Jet 


The total normal stress is a =- pt 2u ou RarseLy, Lt 


is shown that the viscous term is negligible. 


ou. 2) OCEOOIeS > 2 ts) 2 gies 
2u cont lg 5. aa = Toa 2 (by continuity) 
Cay) Sha. a: \ 
oe or dr vee vie Ease , (De) 


where v is velocity paramell to er and ed is its maximum 
in a given vertical occuring at the edge of the boundary layer. When 
r> 0 it is reasonable to aeenne that Fe + (ou_,7/2) mp (as is also 
shown in Chapter IV). This relation leads to Equation 4.35 which can 


be used to evaluate IS Gv, and Equation D.1 becomes: 
r=0 


= < (rv) ~ 226 5 Se (D.2) 
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At the same time, p 42 Pigiod pies Hence it will suffice to show that: 


Und 
Qu 2 VIG k = seep 
p Ss 


or, after some obvious manipulations: 


228 
R oe te Ky = sh) Ws} (is) 


which is always satisfied for a turbulent jet. The above considerations 
apply to normal impingement. In the case of oblique impingement, simi- 
lar reasoning can be applied for wall shear stresses on the symmetry plane. 
In this case the angular gradient is zero and, therefore, the continu- 

ity equation will be identical with that of axisymmetric impingement. 


Using the results of Chapter VII, the condition to be satisfied will be: 


4 /2k/K 
Sb eae SD 
ropes (b/H) (D.4) 


where now A and b/H depend on the angle of impingement height. The 
right-hand side in (D.4) will be less than about 20 if the angle of 
impingement, > is greater than 10°, which implies that (D.4) will al- 


ways be satisfied in this range. 


A similar development can also be applied to the case of small 


2 
impingement heights (Chapter V). Here es pu, /2 and the necessary 


conditions (using results from Chapter V) will read: 


R >> 10 (D.5) 


which again is always satisfied. 
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The starting equation for prediction of wall shear stress will 


be: 
Pee - = Gp (D.6) 


where z is distance from the stagnation point measured along the nor- 
mal to the wall. This relation applies for normal impingement as well 
as for oblique along the symmetry plane (see also Chapters IV, V and VII). 


Assuming similarity of pressure profiles: 


2 
p= pF) soe EID. Mirra 1 a 603 
Differentiating: 
op Dp. ob 
SPive st. S pene es ae 
y pt eg a 
P 
op db 
and is 22) =b (==*) Hip (2) n7¢" (D.7) 
\ oz AO pw \dz 0 s \dz 0 


Using Equation D.7, and integrating Equation D.6, the wall shear stress 


will be given by: 
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For the case of large impingement heights, dimensional reasoning will 


show that: 


b 
ae. Z aes x4 
: hy Ge >) and = hy Gp ) 
OP. Pe db 
It follows that oe, = h, (9) i and Ge, = h, (9), where 
oh oh 
Ba ee ee A oe eee ON 
Ng Laci), ao (prs Lste7 |, ao 
With these, Equation D.8 becomes: 
Tol, &, | Spee of RE (n) - h, (4) g,(n) (D.9) 
fo PEL ps 4 he Digs ON 


If Nn is the value of mn where the maximum value of Th occurs, then 


this maximum, hers is given by: 


it 
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ge foe qv 
se - oe aes h,(o) + 2m, (0 | g,(n,).- bh, (o)-2,n,) (D.10) 


Letting a(¢) and 8(d) represent the coefficients of 81 and 85 in 


Equations D.9 and D.10: 
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Ss 
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ag)'M) - 8B e'() = 0 (D.12) 
g,'(n,) a 
Letting oO = ———=-— and ha [g, ,) - On g,(n_)] the above 


85'(n) 


can also be written as: 


oO = 
ie be [g,(n) - 5, 8,()] 
(D.13) 
ap. 
Ton ube Peatobey Mdesig. 8, Qe] = ied 


which implies that: 


ai H 
Cs sit, =ih@) (D.14) 
pu, 


The quantities qa and On depend on "hn and Figure 7.13 (Chapter V) 
shows these functions graphically. The variation of "mn with the angle 


of impingement 96, is left to be obtained empirically. 


In the case of normal impingement, Le ~ 1.76 which gives 


aes 2.295 and one 0.32. Using also the result cee = 0.078, 


Tle will be given as: 
- 9.43 Xe. (D.15) 


where \} = r/H. This is identical to Equation 4.31 in Chapter IV. 


Consider now the case of oblique impingement. It is reasonable 


to assume that the length scales defining the wall pressure (b, »b, and b4) 
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approach a constant ratio between themselves as the wall is approached, 


i.e., for small values of z, by/b, = (¢. and b,/b, =¢ so that c 


a 2 1 


and Cy) are independent of z. This implies that: 
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a relation that was verified experimentally (Chapter VII). From Equa- 


tions D.13, it now follows that: 
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The condition of continuity of wall shear variation at x = 0, requires 


that (from Equation D.13): 
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It can be easily checked that this is indeed satisfied by means of 


Equation D.18. 


Consider now the case of small impingement heights. The steps 
leading to Equation D.8 are identical, therefore, Equation D.8 will be 


the starting equation in this case. In Chapter V it was shown that: 


2 Ae 
p= pU “/2 and —P—=£,4, 5 
s fo) ou 219 ral (al 
re) 
db 
Therefore, a = he (2/d), hence Cone = const. 
b = OP. 
Further, the term Se om. = const , hence: 
Ps fe) 
% 
— = ag,(n) - B 8, (n) (D.19) 
Pe 1 


where @ and 8 are constants. Using a Gaussian approximation to 


| ee halen . 
P,,/P5> Lee = 0.62. Hence n= 5-65 = 1.85 which gives aes 2.345 


ALE 


/ i) Ny ae i 4 i 
~gpd mort (IL, siaaqegy 
be “ i a! Y Be 


a 


orkbpax <u> x as nabsnsaay sede Live Ye: 


; | to anapi xd sor taties osha af, ats 


ae ie Ig 
4 toe ie 


~ " Al | 


( 2 = ha % zea 
i ’ ; f He ; : : ‘ ¥ iy q I r i " ! ir A as am ® / eT b- +" ie 7 
agasa iad? emai nett ‘eae “Ser eit alka We 


sanity ‘i gH PaeegqhtO AE voune Baits 
ca rete | ; > "i mi ay i ; 
ie e 7 ’ =. : t 
| eas ine eae Ae 
init : ay $4 ‘i ae Shek. 
—_ . sands. wf)” eave 
Hs - ’ J e 1 
. a6 Pr Es 
PF \peotteg. , dado = Ga) . yD 
: I a | 
) ee ih eu 7! eo oe ae : or 7 
a! eae Pe ne Bis anes Wotsi'es aut 
? ; oe te Tah . be . \2 Le, tee ; <_ 


rt 1 «Baal f i Pa ines = re 


i | 93 jen eines sce aa ata ic 4 i - a 


* 


\ a) favia ia a site» By — 


318. 


and On = 0.385. Thus the wall shear stress can finally be written as: 


Z 
pecs , r/d ane 


which is identical to Equation 5.39 of Chapter V. The maximum wall 


shear stress will be simply: 


Dela Plane Jet 


The total normal stress is again e,. See po Zit e (see 


Chapter VI for definition of symbols). From continuity: 
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which will always be satisfied for a turbulent jet. 


The starting equation for the prediction of wall shear stress 


is: 
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This is identical with Equation 6.26 of Chapter VI, if B/a is set 


equal to C Using the subscripts 1 and 2 to denote the sides of 
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positive and negative xy respectively, and reasoning in a manner simi- 


lar to that of the previous section, it is possible to show that: 
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DeZeeDenivataon oF Equations 5.26 and 5.30 


The equation of continuity is: 
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r r 
oe Re ee soit ; io2 
rv = 7 UOT Aust) f dAdA + Au v5) if ¢ A dA 
0 0 
F, (A) FLO) 


Lhismeansalso be wihitten ase 


f ’ ' 13 
v/U : Au, rjd BF, (A) u'd j Aud Lele 
r,/d U ry A 2 Us US r 


(D.30) 


1 
But ui =U £(z) and PA do[f(z)] .2 where o is a function of H/d, 
Oo 


Cc 


related to C as defined in Equation 5.20. Since US ie a Bee Os) uy air 


follows that Au_ = pu f(z) > 6) SND. 15. Substituting these into 
c 


Equation D.30 and simplifying: 
Le 
= - 5 £'(@)FQ) (D.31) 


where F(\) = A+B (Fy + TEDL. This equation is identical to 


Be 3 


ac qd) ees 
~. reg 
> vi cs ka 
ak yt 
| “gatsbageai ins esd sotsoupl & : 
=. : ' “ey 
oe \ ’ i 4390 7 ot ae ~ ah ve ~ a\ 
a, 
e egy Q Dy leer ; : : ‘48 
CA) pee : 
SO a ine 
(NE.ay 9 = eee | eee a 


. | : ) \ i D “4 Mee eae he Raia a | i 
Sb\H (OX iis 3 at. Olexade -& (PSD © t ban (5) 20 su tua 


hie 


; “i i ‘ * aa | 
1 ow Cx if is ie e5are 08.2 sotsappe Serten re, a oJ beiatas . 


asa stort oribiwdtiadie, 28,0 = ay joa * he “gud awollod 


, ‘eae | Maen | oF ae | 
(18 <a) oe oe A poe Np . - ri 
a a “me : ru ms ihe, of. Te 


O9 ali se - ‘aot aap va Fecye +e = wh. = (At “oreude 


Equation 5.26. To evaluate F(\) it is necessary to find FLO) and 


F, (A), which are related to (A) = x7 (2 ~ ay 
r r 
FQ) - frre afore cones! 25° 
0 0 
r 
4 6 
. SS Oe 
PO) -f AX )dA = 5 G 
0 
4 6 
F, + 2F, S42 ae) 
Hence FO) = At B(22° ~ ) and F(L) = 1+68 therefore: 
PEO a | So) 2 2 
2 Seiya) Te well HEA 2D] (D.32) 


which is identical to Equation 5.30 with § = 0.15. 


D.3 Derivation of Equation 6.15 (Chapter VI) 
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With reference to Figure D.1, the velocity qy due to a vortex 


of strength IT |. is: 


qy = qo = T/2tr 


The resultant velocity will be in the direction of the wall and has a 


magnitude of: 


Applying the Bernoulli equation (and using te =a + 07): 
souls - 2 oes Te (D. 33) 


The minimum pressure p will be: 
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which is identical to Equation 6.15. 


D.4 Derivation of Equation 6.19 (Chapter VI) 
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Using the similarity property of Figures 4.5a and 4.5b: 
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a has to be taken equal to one which is in reasonable agreement with 
the theoretical value of 1.065. Using ae = 1.0, Equation D.42 


becomes: 
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which is Equation 6.19. 


D.5 Derivation of Equation 7.17 (Chapter VII) 
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Puprther; vif G5 = 2nm/k , then: 
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D.6 Derivation of Equation 7.22 (Chapter VII) 


Consider firstly the integral: 
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Evaluating the integrals this becomes: 
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D.7 Derivation of Equation WATEY 


The equation of the ellipse is: 
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In terms of polar co-ordinates, (r, 6) this can also be written as: 


(E cos 6 - BMG vs é sin a 


: = oa (D.47) 


Expanding and introducing the ratio r/r, which equals h/h,, this 


a 
becomes: 
WS Be 2 16 
r ali 7 € se Lie KE 
o eae cos 6 + ® - A a.) cos O.+ 
MARY a 2 (D.48) 
A i a De 
+ cae), oo GC) sin’@ = 1 
1 
But r sate, i.e, rj/a = 1+ (e/a). Hence: 


: 2 € 2 a : 2 
cos 0 + (1 the G) sin'6) - 


ig : € 
© ja+D 


2, 
S Zena ee) cos 6 eer -S = 0 
a 


Further, dividing by 1 + (¢/a): 


2 2D. 
oP ree = }oo870 . ©) sin’ol AM AS eS tos 6 (= =) = 0 


This relation can be solved for r/r, to give: 


2 2 He 2 
Z ©) cos 6 + (1 + =) (1 - £) (cos 6 + —= sin 6) 
a a a a b2 


—cos 8 
r 
ae tie 2 


Gr 4 Ses "6 + =. 
a b2 


oe aeuoet eos Ode (afb)- in-8 [1 - (e/a)? ] 
ett ee NS as te td (D.49) 


(1 + ia) Tones 6 2 tarny aac) 


I+ 


ginco) 


. Lee ERM ey i adh Cena 


jt } it 

: ne re 
'en nos i pt on = 

ie 
(viay | | a yor 
; ea : iM, 
eh ald alawpe olde 7\1 okie eth gahsuaRae Bae eee 
Fa ue ni abel : ee on 
. ‘ oa core ae 
> tae f a} 


| +.6 wes ODE - ‘e ihe ‘og i 


aie) 


qe; p ee ; iy, ms xe ag “ : I a 
‘si ee a ‘gh ry ree 

. =. ee ees A te , ae a 

0 '= G “)) v iar tale e nik es + Ff ¢ + a ‘s ; an, 


a oe tpvity oe ais wr vi wd tine ants 


(Ente he Ope ae - ngs D4 oto 2 


i { 
mt acl 8 


332. 


But ¢€/a = cos ¢ , hence: 
i |. Sos $.¢os' 0 2 Weoso + oo mer tang (D.50) 
oy (1 + cos atieoa’o + Care e179] 


Since x/r, is a positive quantity only the "plus" sign is meaningful. 
Recalling that r/r, = h/h,; Equation D.50 becomes identical to Equation 


7570. 


D.8 Derivation of Equation 8.13 


The distance Xo from the point source to the stagnation 


point is given by (56): 


m 
Xo = anu, | (D. 51) 


But this was found to be equal to 0.2 5° hence: 


m = am |u, | (0.2) 


Using Equation 8.14, this becomes: 


m = 4m(0.2)7(2.6)° dé ua / Iu, | 
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which is Equation 8.13. 
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D.9 Stagnation Streamsurface, y = Wits 


The starting equation is: 


350% 


V5 
| 2myudy = O (D.53) 
0 
But: u = uy + u fn) (Equation 8.7) 
Ys Ns 
2 
| yudy = b | nu f + u,)dn 
0 0 
where ns = y,/b. Hence: 
Ys Ns Qs 
| yudy = ue \e, | nidn + uy | nan 
0 0 0 
Evaluating the integrals, this becomes: 
cae Led 1 pnke 
0.59 : Pla Riu 
y1 + 0.59 n 
s 
-1 
’ any? 1 ; : . 
Letting F(n) =n 1 - and using Equations 8.9 and 
4059 rie 
8.14, this reduces to: 
x 
F(n,) heh eee A owe (D.54) 
Pp 
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= function of = only (D.55) 
Pp Pp 


This function can be evaluated easily using Equation D.54 and Equation 


S215. 


